Back to home page

MITgcm

 
 

    


File indexing completed on 2024-07-17 05:10:38 UTC

view on githubraw file Latest commit acacc28f on 2024-07-17 03:59:01 UTC
5ed655852f Jean*0001 #include "COST_OPTIONS.h"
47dd97aab7 Patr*0002 #ifdef ALLOW_AUTODIFF
                0003 # include "AUTODIFF_OPTIONS.h"
                0004 #endif
39154404ca Patr*0005 
acacc28f7f Jean*0006 CBOP
                0007 C     !ROUTINE: COST_FINAL_RESTORE
39154404ca Patr*0008 
acacc28f7f Jean*0009 C     !INTERFACE:
                0010       SUBROUTINE COST_FINAL_RESTORE( lastDiva, myThid )
39154404ca Patr*0011 
acacc28f7f Jean*0012 C     !DESCRIPTION:
                0013 C     Restore cost function value
                0014 C     This is required in conjunction with DIVA
                0015 C     heimbach@mit.edu 5-Mar-2003
39154404ca Patr*0016 
acacc28f7f Jean*0017 C     !USES:
                0018       IMPLICIT NONE
                0019 C     == Global variables ===
39154404ca Patr*0020 #include "EEPARAMS.h"
                0021 #include "SIZE.h"
                0022 #include "cost.h"
                0023 
acacc28f7f Jean*0024 C     !INPUT/OUTPUT PARAMETERS:
                0025 C     lastDiva :: True if last AD run from DIVA sequence, False otherwise
                0026 C     myThid   :: my Thread Id number
                0027       LOGICAL lastDiva
                0028       INTEGER myThid
                0029 
                0030 #ifdef ALLOW_DIVIDED_ADJOINT
                0031 C     !LOCAL VARIABLES:
                0032       LOGICAL exst
                0033       INTEGER idivbeg, idivend
                0034 CEOP
                0035 
                0036       _BEGIN_MASTER( myThid )
                0037 
                0038         INQUIRE(FILE='divided.ctrl',EXIST=exst)
                0039         IF (exst) THEN
                0040           OPEN(UNIT=76,FILE='divided.ctrl',FORM='formatted')
                0041           READ(UNIT=76,FMT=*) idivbeg,idivend
                0042           CLOSE(UNIT=76)
                0043         ELSE
                0044           idivbeg = -2
                0045         ENDIF
                0046 
                0047         IF ( idivbeg .EQ. 0 ) THEN
                0048           lastDiva = .TRUE.
                0049           OPEN(UNIT=76,FILE='costfinal')
                0050           READ(76,*) fc
                0051           CLOSE(76)
                0052         ELSE
                0053           lastDiva = .FALSE.
                0054         ENDIF
                0055 
                0056       _END_MASTER( myThid )
39154404ca Patr*0057       _BARRIER
                0058 
acacc28f7f Jean*0059 #endif /* ALLOW_DIVIDED_ADJOINT */
39154404ca Patr*0060 
acacc28f7f Jean*0061       RETURN
                0062       END