/[MITgcm]/MITgcm/pkg/dic/dic_read_pickup.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/dic_read_pickup.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.4 by dfer, Mon Apr 7 20:31:16 2008 UTC revision 1.5 by jmc, Thu Apr 10 00:49:25 2008 UTC
# Line 16  C     === Global variables === Line 16  C     === Global variables ===
16    
17  C     == Routine arguments ==  C     == Routine arguments ==
18  C     myThid  :: my Thread Id number  C     myThid  :: my Thread Id number
19        LOGICAL ph_isLoaded        LOGICAL pH_isLoaded
20        INTEGER myIter        INTEGER myIter
21        INTEGER myThid        INTEGER myThid
22    
# Line 32  C     == Local variables == Line 32  C     == Local variables ==
32        CHARACTER*(MAX_LEN_FNAM) fn, dFileName        CHARACTER*(MAX_LEN_FNAM) fn, dFileName
33        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
34        LOGICAL fileExist, shareExist        LOGICAL fileExist, shareExist
35          INTEGER iG, jG
36        INTEGER fp, iL, i, ioUnit        INTEGER fp, iL, i, ioUnit
37    
38        COMMON / LOCAL_DIC_PICKUP / shareExist        COMMON / LOCAL_DIC_PICKUP / shareExist
39  CEOP  CEOP
40    
41        ph_isLoaded =.FALSE.        pH_isLoaded =.FALSE.
42        ioUnit = errorMessageUnit        ioUnit = errorMessageUnit
43    
44  C--   pickup file name :  C--   pickup file name :
# Line 67  C-    look for file = {fn}'.data' Line 68  C-    look for file = {fn}'.data'
68          INQUIRE( FILE=dFileName, EXIST=fileExist )          INQUIRE( FILE=dFileName, EXIST=fileExist )
69        ENDIF        ENDIF
70        IF ( .NOT.fileExist ) THEN        IF ( .NOT.fileExist ) THEN
71    C-    look for file = {fn}'.{iG}.{jG}.meta'
72            iG = 1+(myXGlobalLo-1)/sNx
73            jG = 1+(myYGlobalLo-1)/sNy
74            WRITE(dFileName,'(2A,I3.3,A,I3.3,A)')
75         &                fn(1:iL), '.', iG, '.', jG, '.data'
76            i = iL + 5 + 8
77            INQUIRE( FILE=dFileName, EXIST=fileExist )
78          ENDIF
79          IF ( .NOT.fileExist ) THEN
80  C-    look for file = {fn}'.001.001.data'  C-    look for file = {fn}'.001.001.data'
81          WRITE(dFileName,'(2A)') fn(1:iL), '.001.001.data'          WRITE(dFileName,'(2A)') fn(1:iL), '.001.001.data'
82          i = iL + 5 + 8          i = iL + 5 + 8
# Line 101  C--   Read pickup file Line 111  C--   Read pickup file
111            CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )            CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
112            STOP 'ABNORMAL END: S/R DIC_READ_PICKUP'            STOP 'ABNORMAL END: S/R DIC_READ_PICKUP'
113          ELSE          ELSE
114            WRITE(msgBuf,'(4A)') 'DIC_READ_PICKUP: ',            WRITE(msgBuf,'(2A)') 'WARNING >> DIC_READ_PICKUP: ',
115       &      'will restart from approximated pH'       &      'will restart from approximated pH'
116            CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )            CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
117          ENDIF          ENDIF
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22