/[MITgcm]/MITgcm_contrib/rpa_layers/layers/layers_readparms.F
ViewVC logotype

Diff of /MITgcm_contrib/rpa_layers/layers/layers_readparms.F

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

revision 1.1 by rpa, Tue Sep 15 19:16:53 2009 UTC revision 1.2 by jmc, Wed Sep 16 18:04:49 2009 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "LAYERS_OPTIONS.h"  #include "LAYERS_OPTIONS.h"
5    
6  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
# Line 31  C     k           - index Line 34  C     k           - index
34        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
35    
36  C--   Default values for LAYERS  C--   Default values for LAYERS
37          
38  C     The MNC stuff is not working yet        C     The MNC stuff is not working yet
39        layers_MNC = .FALSE.        layers_MNC = .FALSE.
40        layers_MDSIO = .TRUE.        layers_MDSIO = .TRUE.
41          
42        DO k=1,Nlayers+1        DO k=1,Nlayers+1
43          layers_G(k) = UNSET_RL          layers_G(k) = UNSET_RL
44        ENDDO        ENDDO
45          
46        WRITE(msgBuf,'(A)') 'LAYERS_READPARMS: opening data.layers'        WRITE(msgBuf,'(A)') 'LAYERS_READPARMS: opening data.layers'
47        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
48       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
# Line 59  C     Close the open data file Line 62  C     Close the open data file
62    
63  C--   Make sure the layers_G we just read is big enough  C--   Make sure the layers_G we just read is big enough
64        DO k=1,Nlayers+1        DO k=1,Nlayers+1
65          IF ( layers_G(k) .EQ. UNSET_RL ) THEN                  IF ( layers_G(k) .EQ. UNSET_RL ) THEN
66            WRITE(msgBuf,'(A,I4)')            WRITE(msgBuf,'(A,I4)')
67       &        'S/R LAYERS_READPARMS: No value for layers_G at k =', k       &        'S/R LAYERS_READPARMS: No value for layers_G at k =', k
68            CALL PRINT_ERROR( msgBuf, myThid )            CALL PRINT_ERROR( msgBuf, myThid )
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22