/[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.1 by dimitri, Fri May 30 21:51:24 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,
48         &          GM_InMomAsStress,
49       &          GM_isopycK,       &          GM_isopycK,
50       &          GM_background_K,       &          GM_background_K,
51       &          GM_taper_scheme,       &          GM_taper_scheme,
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_Scrit, GM_Sd,
61       &          GM_MNC       &          GM_MNC
62  c     &          GM_dumpFreq,  C SM (3)
63  c     &          GM_taveFreq       &          ,GM_SM_Ce,
64         &          GM_SM_Lf,
65         &          GM_SM_tau
66                  
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.
77        GM_AdvSeparate      = .FALSE.        GM_AdvSeparate      = .FALSE.
78          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
92    CBFK Submesoscale Paramters
93    C     GM_SM_Ce is the efficiency factor.  Nondimensional, typically 0.06
94    C     Will be converted to a dimensional factor later
95    C     Set to zero to kill Submeso at runtime.
96          GM_SM_Ce            = 0.06  _d 0
97    C     GM_SM_Lf is the minimum Front Width scale.  Typically 1-10km
98          GM_SM_Lf            = 5.00  _d 3
99    C     GM_SM_Lmax is the maximum upscale.  Typically 100km
100          GM_SM_Lmax            = 1.00  _d 5
101    C     GM_SM_tau is the momentum mixing timescale: 1/d to 1/week
102          GM_SM_tau           = 8.64  _d 4
103    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 89  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 110  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.
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22