/[MITgcm]/MITgcm_contrib/submesoscale/code/gmredi_readparms.F
ViewVC logotype

Diff of /MITgcm_contrib/submesoscale/code/gmredi_readparms.F

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

revision 1.3 by dimitri, Fri May 30 22:24:25 2008 UTC revision 1.4 by zhc, Fri Mar 12 18:31:00 2010 UTC
# Line 45  C     GM_Small_Number :: epsilon used in Line 45  C     GM_Small_Number :: epsilon used in
45  C     GM_slopeSqCutoff :: slope^2 cut-off value  C     GM_slopeSqCutoff :: slope^2 cut-off value
46        NAMELIST /GM_PARM01/        NAMELIST /GM_PARM01/
47       &          GM_AdvForm, GM_AdvSeparate,       &          GM_AdvForm, GM_AdvSeparate,
      &          GM_ExtraDiag,  
48       &          GM_InMomAsStress,       &          GM_InMomAsStress,
49       &          GM_isopycK,       &          GM_isopycK,
50       &          GM_background_K,       &          GM_background_K,
# Line 53  C     GM_slopeSqCutoff :: slope^2 cut-of Line 52  C     GM_slopeSqCutoff :: slope^2 cut-of
52       &          GM_maxSlope,       &          GM_maxSlope,
53       &          GM_Kmin_horiz,       &          GM_Kmin_horiz,
54       &          GM_Small_Number, GM_slopeSqCutoff,       &          GM_Small_Number, GM_slopeSqCutoff,
55       &          GM_Visbeck_alpha,       &          GM_Visbeck_alpha, GM_Visbeck_length,
      &          GM_Visbeck_length,  
56       &          GM_Visbeck_depth,       &          GM_Visbeck_depth,
57       &          GM_Visbeck_maxval_K,       &          GM_Visbeck_minDepth, GM_Visbeck_maxSlope,
58       &          GM_Scrit,       &          GM_Visbeck_minVal_K, GM_Visbeck_maxVal_K,
59       &          GM_Sd,       &          GM_facTrL2dz, GM_facTrL2ML, GM_maxTransLay,
60       &          GM_MNC,       &          GM_Scrit, GM_Sd,
61       &          GM_SM_Ce,       &          GM_MNC
62    C SM (3)
63         &          ,GM_SM_Ce,
64       &          GM_SM_Lf,       &          GM_SM_Lf,
65       &          GM_SM_tau       &          GM_SM_tau
66  c     &          GM_dumpFreq,                
 c     &          GM_taveFreq  
67    
68        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
69    
70  C--   GMREDI_READPARMS has been called so we know that  C--   GMREDI_READPARMS has been called so we know that the package is active.
71  C     the package is active.  C note: cannot test for GMRediIsOn outside GM-redi S/R (since is in GMREDI.h)
72        GMRediIsOn = .TRUE.  C       + is not initialized (if useGMRedi=F) => useless => removed
73    c     GMRediIsOn = .TRUE.
74    
75  C--   Default values GM/Redi  C--   Default values GM/Redi
76        GM_AdvForm          = .FALSE.        GM_AdvForm          = .FALSE.
       GM_ExtraDiag        = .FALSE.  
77        GM_AdvSeparate      = .FALSE.        GM_AdvSeparate      = .FALSE.
78        GM_InMomAsStress    = .FALSE.        GM_InMomAsStress    = .FALSE.
79        GM_isopycK          = -999.        GM_isopycK          = -999.
80        GM_background_K     = 0. _d 0        GM_background_K     = 0. _d 0
81        GM_maxslope         = 1. _d -2        GM_maxSlope         = 1. _d -2
82        GM_Kmin_horiz       = 0. _d 0        GM_Kmin_horiz       = 0. _d 0
83        GM_Small_Number     = 1. _d -20        GM_Small_Number     = 1. _d -20
84        GM_slopeSqCutoff    = 1. _d +48        GM_slopeSqCutoff    = 1. _d +48
85        GM_taper_scheme     = ' '        GM_taper_scheme     = ' '
86          GM_facTrL2dz        = 1.
87          GM_facTrL2ML        = 5.
88          GM_maxTransLay      = 500.
89        GM_Scrit            = 0.004 _d 0        GM_Scrit            = 0.004 _d 0
90        GM_Sd               = 0.001 _d 0        GM_Sd               = 0.001 _d 0
91        GM_MNC              = useMNC        GM_MNC              = useMNC
# Line 94  C     Set to zero to kill Submeso at run Line 96  C     Set to zero to kill Submeso at run
96        GM_SM_Ce            = 0.06  _d 0        GM_SM_Ce            = 0.06  _d 0
97  C     GM_SM_Lf is the minimum Front Width scale.  Typically 1-10km  C     GM_SM_Lf is the minimum Front Width scale.  Typically 1-10km
98        GM_SM_Lf            = 5.00  _d 3        GM_SM_Lf            = 5.00  _d 3
99  C     GM_SM_Lmax is the maximum upscale.  Typically 100km  C     GM_SM_Lmax is the maximum upscale.  Typically 100km
100        GM_SM_Lmax            = 1.00  _d 5        GM_SM_Lmax            = 1.00  _d 5
101  C     GM_SM_tau is the momentum mixing timescale: 1/d to 1/week  C     GM_SM_tau is the momentum mixing timescale: 1/d to 1/week
102        GM_SM_tau           = 8.64  _d 4        GM_SM_tau           = 8.64  _d 4
103  C     Will be converted to an inverse squared timescale later  C     Will be converted to an inverse squared timescale later
104                            
105    
106  C--   Default values GM/Redi I/O control  C--   Default values GM/Redi I/O control
107  c     GM_dumpFreq         = -1.  c     GM_dumpFreq         = -1.
# Line 108  C--   Default values Visbeck Line 111  C--   Default values Visbeck
111        GM_Visbeck_alpha    =    0. _d 0        GM_Visbeck_alpha    =    0. _d 0
112        GM_Visbeck_length   =  200. _d 3        GM_Visbeck_length   =  200. _d 3
113        GM_Visbeck_depth    = 1000. _d 0        GM_Visbeck_depth    = 1000. _d 0
114        GM_Visbeck_maxval_K = 2500. _d 0        GM_Visbeck_minDepth =    0. _d 0
115          GM_Visbeck_maxSlope = UNSET_RL
116          GM_Visbeck_minVal_K =    0. _d 0
117          GM_Visbeck_maxVal_K = 2500. _d 0
118    
119        WRITE(msgBuf,'(A)') ' GM_READPARMS: opening data.gmredi'        WRITE(msgBuf,'(A)') ' GM_READPARMS: opening data.gmredi'
120        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
# Line 129  C     Close the open data file Line 135  C     Close the open data file
135  C     Default value for GM_isopycK is equal to GM_background_K :  C     Default value for GM_isopycK is equal to GM_background_K :
136        IF (GM_isopycK.EQ.-999.) GM_isopycK = GM_background_K        IF (GM_isopycK.EQ.-999.) GM_isopycK = GM_background_K
137    
138    C     Default value for GM_Visbeck_maxSlope is equal to GM_maxSlope :
139          IF ( GM_Visbeck_maxSlope .EQ. UNSET_RL )
140         &     GM_Visbeck_maxSlope = GM_maxSlope
141    
142  C     Some constants  C     Some constants
143        GM_rMaxSlope=0.        GM_rMaxSlope = 0.
144        if (GM_maxSlope.ne.0.) GM_rMaxSlope=1. / GM_maxSlope        if (GM_maxSlope.NE.0.) GM_rMaxSlope = 1. _d 0 / GM_maxSlope
145    
146        IF (GM_AdvForm) THEN        IF (GM_AdvForm) THEN
147          GM_skewflx = 0.          GM_skewflx = 0.
# Line 140  C     Some constants Line 150  C     Some constants
150        ELSE        ELSE
151          GM_skewflx = 1.          GM_skewflx = 1.
152          GM_advect  = 0.          GM_advect  = 0.
153          GM_ExtraDiag = (GM_isopycK.NE.GM_background_K)          GM_ExtraDiag = GM_isopycK.NE.GM_background_K
      &      .or.(GM_SM_Ce.gt.0 _d 0)  
154        ENDIF        ENDIF
155    
156  C     Make sure that we locally honor the global MNC on/off flag  C     Make sure that we locally honor the global MNC on/off flag
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22