/[MITgcm]/MITgcm/pkg/ctrl/ctrl_init.F
ViewVC logotype

Diff of /MITgcm/pkg/ctrl/ctrl_init.F

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

revision 1.1 by heimbach, Sun Mar 25 22:33:55 2001 UTC revision 1.2 by heimbach, Fri Jul 13 13:40:17 2001 UTC
# Line 72  c     == local variables == Line 72  c     == local variables ==
72        integer startrec        integer startrec
73        integer endrec        integer endrec
74    
75          _RL nwetc3d
76    
77        character*(max_len_prec) record        character*(max_len_prec) record
78        character*(max_len_mbuf) msgbuf        character*(max_len_mbuf) msgbuf
79    
# Line 86  c--   Read the namelist input. Line 88  c--   Read the namelist input.
88        namelist /ctrl_nml/        namelist /ctrl_nml/
89       &  xx_theta_file,       &  xx_theta_file,
90       &  xx_salt_file,       &  xx_salt_file,
91         &  xx_tr1_file,
92       &  xx_tauu_file,       &  xx_tauu_file,
93       &  xx_tauv_file,       &  xx_tauv_file,
94       &  xx_sflux_file,       &  xx_sflux_file,
# Line 113  c--   Read the namelist input. Line 116  c--   Read the namelist input.
116  c--     Set default values.  c--     Set default values.
117          xx_theta_file      = ' '          xx_theta_file      = ' '
118          xx_salt_file       = ' '          xx_salt_file       = ' '
119            xx_tr1_file        = ' '
120          xx_tauu_file       = ' '          xx_tauu_file       = ' '
121          xx_tauv_file       = ' '          xx_tauv_file       = ' '
122          xx_sflux_file      = ' '          xx_sflux_file      = ' '
# Line 257  c--   ===================== Line 261  c--   =====================
261        _END_MASTER( mythid )        _END_MASTER( mythid )
262  #endif /* ALLOW_TAUV0_CONTROL */  #endif /* ALLOW_TAUV0_CONTROL */
263    
264    cph(
265    cph    index  7-10 reserved for atmos. state,
266    cph    index 11-14 reserved for open boundaries,
267    cph    index 15-16 reserved for mixing coeff.
268    cph    index 17    reserved for passive tracer TR1
269    cph    index 18,19 reserved for sst, sss
270    cph)
271    
272    #ifdef ALLOW_TR10_CONTROL
273          _BEGIN_MASTER( mythid )
274            ncvarindex(17) = 117
275            ncvarrecs(17)  =   1
276            ncvarxmax(17)  = snx
277            ncvarymax(17)  = sny
278            ncvarnrmax(17) =  nr
279            ncvargrd(17)   = 'c'
280          _END_MASTER( mythid )
281    #endif /* ALLOW_TR10_CONTROL */
282    
283  #ifdef ALLOW_SST0_CONTROL  #ifdef ALLOW_SST0_CONTROL
284        _BEGIN_MASTER( mythid )        _BEGIN_MASTER( mythid )
285          ncvarindex(7) = 107          ncvarindex(18) = 118
286          ncvarrecs(7)  =   1          ncvarrecs(18)  =   1
287          ncvarxmax(7)  = snx          ncvarxmax(18)  = snx
288          ncvarymax(7)  = sny          ncvarymax(18)  = sny
289          ncvarnrmax(7) =   1          ncvarnrmax(18) =   1
290          ncvargrd(7)   = 'c'          ncvargrd(18)   = 'c'
291        _END_MASTER( mythid )        _END_MASTER( mythid )
292  #endif /* ALLOW_SST0_CONTROL */  #endif /* ALLOW_SST0_CONTROL */
293    
294  #ifdef ALLOW_SSS0_CONTROL  #ifdef ALLOW_SSS0_CONTROL
295        _BEGIN_MASTER( mythid )        _BEGIN_MASTER( mythid )
296          ncvarindex(8) = 108          ncvarindex(19) = 119
297          ncvarrecs(8)  =   1          ncvarrecs(19)  =   1
298          ncvarxmax(8)  = snx          ncvarxmax(19)  = snx
299          ncvarymax(8)  = sny          ncvarymax(19)  = sny
300          ncvarnrmax(8) =   1          ncvarnrmax(19) =   1
301          ncvargrd(8)   = 'c'          ncvargrd(19)   = 'c'
302        _END_MASTER( mythid )        _END_MASTER( mythid )
303  #endif /* ALLOW_SSS0_CONTROL */  #endif /* ALLOW_SSS0_CONTROL */
304    
# Line 359  c--   Determine the total number of cont Line 382  c--   Determine the total number of cont
382           endif           endif
383        enddo        enddo
384                
385    cph(
386          print *, 'ph-wet 1: nvarlength = ', nvarlength
387          print *, 'ph-wet 2: surface wet C = ', nwetctile(1,1,1)
388          print *, 'ph-wet 3: surface wet W = ', nwetwtile(1,1,1)
389          print *, 'ph-wet 4: surface wet S = ', nwetstile(1,1,1)
390          nwetc3d = 0
391          do k = 1, Nr
392             nwetc3d = nwetc3d + nwetctile(1,1,k)
393          end do
394          print *, 'ph-wet 5: 3D center wet points = ', nwetc3d
395          do i = 1, 6
396             print *, 'ph-wet 6: no recs for i = ', i, ncvarrecs(i)
397          end do
398          print *, 'ph-wet 7: ',
399         &     2*nwetc3d +
400         &     ncvarrecs(3)*nwetctile(1,1,1) +
401         &     ncvarrecs(4)*nwetctile(1,1,1) +
402         &     ncvarrecs(5)*nwetwtile(1,1,1) +
403         &     ncvarrecs(6)*nwetstile(1,1,1)
404    cph)
405    
406  c  c
407  c     Summation of wet point counters  c     Summation of wet point counters
408  c  c
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22