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

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

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

revision 1.18 by heimbach, Mon Feb 28 17:29:38 2005 UTC revision 1.28 by heimbach, Wed Jun 7 01:55:14 2006 UTC
# Line 112  c--   Assign file names. Line 112  c--   Assign file names.
112        call ctrl_set_fname(xx_tauv_file, fname_tauv, mythid)        call ctrl_set_fname(xx_tauv_file, fname_tauv, mythid)
113        call ctrl_set_fname(xx_atemp_file, fname_atemp, mythid)        call ctrl_set_fname(xx_atemp_file, fname_atemp, mythid)
114        call ctrl_set_fname(xx_aqh_file, fname_aqh, mythid)        call ctrl_set_fname(xx_aqh_file, fname_aqh, mythid)
115          call ctrl_set_fname(xx_precip_file, fname_precip, mythid)
116          call ctrl_set_fname(xx_swflux_file, fname_swflux, mythid)
117          call ctrl_set_fname(xx_swdown_file, fname_swdown, mythid)
118        call ctrl_set_fname(xx_uwind_file, fname_uwind, mythid)        call ctrl_set_fname(xx_uwind_file, fname_uwind, mythid)
119        call ctrl_set_fname(xx_vwind_file, fname_vwind, mythid)        call ctrl_set_fname(xx_vwind_file, fname_vwind, mythid)
120        call ctrl_set_fname(xx_obcsn_file, fname_obcsn, mythid)        call ctrl_set_fname(xx_obcsn_file, fname_obcsn, mythid)
# Line 123  c--   Assign file names. Line 126  c--   Assign file names.
126        call ctrl_set_fname(xx_tr1_file, fname_tr1, mythid)        call ctrl_set_fname(xx_tr1_file, fname_tr1, mythid)
127        call ctrl_set_fname(xx_sst_file, fname_sst, mythid)        call ctrl_set_fname(xx_sst_file, fname_sst, mythid)
128        call ctrl_set_fname(xx_sss_file, fname_sss, mythid)        call ctrl_set_fname(xx_sss_file, fname_sss, mythid)
129        call ctrl_set_fname(xx_hfacc_file, fname_hfacc, mythid)        call ctrl_set_fname(xx_depth_file, fname_depth, mythid)
130        call ctrl_set_fname(xx_efluxy_file, fname_efluxy, mythid)        call ctrl_set_fname(xx_efluxy_file, fname_efluxy, mythid)
131        call ctrl_set_fname(xx_efluxp_file, fname_efluxp, mythid)        call ctrl_set_fname(xx_efluxp_file, fname_efluxp, mythid)
132        call ctrl_set_fname(xx_bottomdrag_file, fname_bottomdrag, mythid)        call ctrl_set_fname(xx_bottomdrag_file, fname_bottomdrag, mythid)
133        call ctrl_set_fname(xx_edtaux_file, fname_edtaux, mythid)        call ctrl_set_fname(xx_edtaux_file, fname_edtaux, mythid)
134        call ctrl_set_fname(xx_edtauy_file, fname_edtauy, mythid)        call ctrl_set_fname(xx_edtauy_file, fname_edtauy, mythid)
135          call ctrl_set_fname(xx_uvel_file, fname_uvel, mythid)
136          call ctrl_set_fname(xx_vvel_file, fname_vvel, mythid)
137          call ctrl_set_fname(xx_etan_file, fname_etan, mythid)
138          call ctrl_set_fname(xx_relaxsst_file, fname_relaxsst, mythid)
139          call ctrl_set_fname(xx_relaxsss_file, fname_relaxsss, mythid)
140    
141  c  c--   Only the master thread will do I/O.
 c--     Only the master thread will do I/O.  
142        _BEGIN_MASTER( mythid )        _BEGIN_MASTER( mythid )
143    
144        if ( first ) then        if ( first ) then
# Line 202  C     place holder of obsolete variable Line 209  C     place holder of obsolete variable
209  #ifdef ALLOW_THETA0_CONTROL  #ifdef ALLOW_THETA0_CONTROL
210            ivartype = 1            ivartype = 1
211            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
212            write(weighttype(1:80),'(a)') "wtheta"            write(weighttype(1:80),'(a)') "wthetaLev"
213            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
214       &         cunit, ivartype, fname_theta(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_theta(ictrlgrad), "maskCtrlC",
215       &         weighttype, wtheta, lxxadxx, mythid)       &         weighttype, wtheta, lxxadxx, mythid)
216  #endif  #endif
217    
218  #ifdef ALLOW_SALT0_CONTROL  #ifdef ALLOW_SALT0_CONTROL
219            ivartype = 2            ivartype = 2
220            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
221            write(weighttype(1:80),'(a)') "wsalt"            write(weighttype(1:80),'(a)') "wsaltLev"
222            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
223       &         cunit, ivartype, fname_salt(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_salt(ictrlgrad), "maskCtrlC",
224       &         weighttype, wsalt, lxxadxx, mythid)       &         weighttype, wsalt, lxxadxx, mythid)
225  #endif  #endif
226    
227  #if (defined (ALLOW_HFLUX_CONTROL) || \  #if (defined (ALLOW_HFLUX_CONTROL) || defined (ALLOW_HFLUX0_CONTROL))
      defined (ALLOW_HFLUX0_CONTROL))  
228            ivartype = 3            ivartype = 3
229            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
230            write(weighttype(1:80),'(a)') "whflux"            write(weighttype(1:80),'(a)') "whflux"
231            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
232       &         cunit, ivartype, fname_hflux(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_hflux(ictrlgrad), "maskCtrlC",
233       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
234  #endif  #endif
235    
236  #if (defined (ALLOW_SFLUX_CONTROL) || \  #if (defined (ALLOW_SFLUX_CONTROL) || defined (ALLOW_SFLUX0_CONTROL))
      defined (ALLOW_SFLUX0_CONTROL))  
237            ivartype = 4            ivartype = 4
238            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
239            write(weighttype(1:80),'(a)') "wsflux"            write(weighttype(1:80),'(a)') "wsflux"
240            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
241       &         cunit, ivartype, fname_sflux(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_sflux(ictrlgrad), "maskCtrlC",
242       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
243  #endif  #endif
244    
245  #if (defined (ALLOW_USTRESS_CONTROL) || \  #if (defined (ALLOW_USTRESS_CONTROL) || defined (ALLOW_TAUU0_CONTROL))
      defined (ALLOW_TAUU0_CONTROL))  
246            ivartype = 5            ivartype = 5
247            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
248            write(weighttype(1:80),'(a)') "wtauu"            write(weighttype(1:80),'(a)') "wtauu"
249            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
250       &         cunit, ivartype, fname_tauu(ictrlgrad), "maskW",       &         cunit, ivartype, fname_tauu(ictrlgrad), "maskCtrlW",
251       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
252  #endif  #endif
253    
254  #if (defined (ALLOW_VSTRESS_CONTROL) || \  #if (defined (ALLOW_VSTRESS_CONTROL) || defined (ALLOW_TAUV0_CONTROL))
      defined (ALLOW_TAUV0_CONTROL))  
255            ivartype = 6            ivartype = 6
256            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
257            write(weighttype(1:80),'(a)') "wtauv"            write(weighttype(1:80),'(a)') "wtauv"
258            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
259       &         cunit, ivartype, fname_tauv(ictrlgrad), "maskS",       &         cunit, ivartype, fname_tauv(ictrlgrad), "maskCtrlS",
260       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
261  #endif  #endif
262    
# Line 262  C     place holder of obsolete variable Line 265  C     place holder of obsolete variable
265            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
266            write(weighttype(1:80),'(a)') "watemp"            write(weighttype(1:80),'(a)') "watemp"
267            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
268       &         cunit, ivartype, fname_atemp(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_atemp(ictrlgrad), "maskCtrlC",
269       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
270  #endif  #endif
271    
# Line 271  C     place holder of obsolete variable Line 274  C     place holder of obsolete variable
274            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
275            write(weighttype(1:80),'(a)') "waqh"            write(weighttype(1:80),'(a)') "waqh"
276            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
277       &         cunit, ivartype, fname_aqh(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_aqh(ictrlgrad), "maskCtrlC",
278       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
279  #endif  #endif
280    
# Line 280  C     place holder of obsolete variable Line 283  C     place holder of obsolete variable
283            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
284            write(weighttype(1:80),'(a)') "wuwind"            write(weighttype(1:80),'(a)') "wuwind"
285            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
286       &         cunit, ivartype, fname_uwind(ictrlgrad), "maskW",       &         cunit, ivartype, fname_uwind(ictrlgrad), "maskCtrlC",
287       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
288  #endif  #endif
289    
# Line 289  C     place holder of obsolete variable Line 292  C     place holder of obsolete variable
292            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
293            write(weighttype(1:80),'(a)') "wvwind"            write(weighttype(1:80),'(a)') "wvwind"
294            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
295       &         cunit, ivartype, fname_vwind(ictrlgrad), "maskS",       &         cunit, ivartype, fname_vwind(ictrlgrad), "maskCtrlC",
296       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
297  #endif  #endif
298    
# Line 334  C     place holder of obsolete variable Line 337  C     place holder of obsolete variable
337            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
338            write(weighttype(1:80),'(a)') "wdiffkr"            write(weighttype(1:80),'(a)') "wdiffkr"
339            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
340       &         cunit, ivartype, fname_diffkr(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_diffkr(ictrlgrad), "maskCtrlC",
341       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, wdiffkr, lxxadxx, mythid)
342  #endif  #endif
343    
344  #ifdef ALLOW_KAPGM_CONTROL  #ifdef ALLOW_KAPGM_CONTROL
# Line 343  C     place holder of obsolete variable Line 346  C     place holder of obsolete variable
346            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
347            write(weighttype(1:80),'(a)') "wkapgm"            write(weighttype(1:80),'(a)') "wkapgm"
348            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
349       &         cunit, ivartype, fname_kapgm(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_kapgm(ictrlgrad), "maskCtrlC",
350       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, wkapgm, lxxadxx, mythid)
351  #endif  #endif
352    
353  #ifdef ALLOW_TR10_CONTROL  #ifdef ALLOW_TR10_CONTROL
# Line 352  C     place holder of obsolete variable Line 355  C     place holder of obsolete variable
355            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
356            write(weighttype(1:80),'(a)') "wtr1"            write(weighttype(1:80),'(a)') "wtr1"
357            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
358       &         cunit, ivartype, fname_tr1(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_tr1(ictrlgrad), "maskCtrlC",
359       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, wunit, lxxadxx, mythid)
360  #endif  #endif
361    
362  #ifdef ALLOW_SST0_CONTROL  #if (defined (ALLOW_SST_CONTROL) || defined (ALLOW_SST0_CONTROL))
363            ivartype = 18            ivartype = 18
364            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
365            write(weighttype(1:80),'(a)') "wsst0"            write(weighttype(1:80),'(a)') "wsst"
366            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
367       &         cunit, ivartype, fname_sst(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_sst(ictrlgrad), "maskCtrlC",
368       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
369  #endif  #endif
370    
371  #ifdef ALLOW_SSS0_CONTROL  #if (defined (ALLOW_SSS_CONTROL) || defined (ALLOW_SSS0_CONTROL))
372            ivartype = 19            ivartype = 19
373            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
374            write(weighttype(1:80),'(a)') "wsss0"            write(weighttype(1:80),'(a)') "wsss"
375            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
376       &         cunit, ivartype, fname_sss(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_sss(ictrlgrad), "maskCtrlC",
377       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
378  #endif  #endif
379    
380  #ifdef ALLOW_HFACC_CONTROL  #ifdef ALLOW_DEPTH_CONTROL
381            ivartype = 20            ivartype = 20
382            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
383            write(weighttype(1:80),'(a)') "whfacc"            write(weighttype(1:80),'(a)') "wdepth"
 # ifdef ALLOW_HFACC3D_CONTROL  
           call ctrl_set_pack_xyz(  
      &         cunit, ivartype, fname_hfacc(ictrlgrad), "hFacC",  
      &         weighttype, wunit, lxxadxx, mythid)  
 # else  
384            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
385       &         cunit, ivartype, fname_hfacc(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_depth(ictrlgrad), "maskCtrlC",
386       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
387  # endif  #endif /* ALLOW_DEPTH_CONTROL */
 #endif  
388    
389  #ifdef ALLOW_EFLUXY0_CONTROL  #ifdef ALLOW_EFLUXY0_CONTROL
390            ivartype = 21            ivartype = 21
391            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
392            write(weighttype(1:80),'(a)') "wefluxy0"            write(weighttype(1:80),'(a)') "wefluxy0"
393            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
394       &         cunit, ivartype, fname_efluxy(ictrlgrad), "hFacS",       &         cunit, ivartype, fname_efluxy(ictrlgrad), "maskCtrlS",
395       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, wunit, lxxadxx, mythid)
396  #endif  #endif
397    
# Line 403  C     place holder of obsolete variable Line 400  C     place holder of obsolete variable
400            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
401            write(weighttype(1:80),'(a)') "wefluxp0"            write(weighttype(1:80),'(a)') "wefluxp0"
402            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
403       &         cunit, ivartype, fname_efluxp(ictrlgrad), "hFacV",       &         cunit, ivartype, fname_efluxp(ictrlgrad), "maskhFacV",
404       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, wunit, lxxadxx, mythid)
405  #endif  #endif
406    
# Line 412  C     place holder of obsolete variable Line 409  C     place holder of obsolete variable
409            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
410            write(weighttype(1:80),'(a)') "wbottomdrag"            write(weighttype(1:80),'(a)') "wbottomdrag"
411            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
412       &         cunit, ivartype, fname_bottomdrag(ictrlgrad), "hFacC",       &      cunit, ivartype, fname_bottomdrag(ictrlgrad), "maskCtrlC",
413       &         weighttype, lxxadxx, mythid)       &      weighttype, lxxadxx, mythid)
414  #endif  #endif
415    
416  #ifdef ALLOW_EDTAUX_CONTROL  #ifdef ALLOW_EDTAUX_CONTROL
# Line 421  C     place holder of obsolete variable Line 418  C     place holder of obsolete variable
418            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
419            write(weighttype(1:80),'(a)') "wedtaux"            write(weighttype(1:80),'(a)') "wedtaux"
420            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
421       &         cunit, ivartype, fname_edtaux(ictrlgrad), "maskW",       &         cunit, ivartype, fname_edtaux(ictrlgrad), "maskCtrlW",
422       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, wedtaux, lxxadxx, mythid)
423  #endif  #endif
424    
425  #ifdef ALLOW_EDTAUY_CONTROL  #ifdef ALLOW_EDTAUY_CONTROL
# Line 430  C     place holder of obsolete variable Line 427  C     place holder of obsolete variable
427            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
428            write(weighttype(1:80),'(a)') "wedtauy"            write(weighttype(1:80),'(a)') "wedtauy"
429            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
430       &         cunit, ivartype, fname_edtauy(ictrlgrad), "maskS",       &         cunit, ivartype, fname_edtauy(ictrlgrad), "maskCtrlS",
431         &         weighttype, wedtauy, lxxadxx, mythid)
432    #endif
433    
434    #ifdef ALLOW_UVEL0_CONTROL
435              ivartype = 27
436              write(weighttype(1:80),'(80a)') ' '
437              write(weighttype(1:80),'(a)') "wuvel"
438              call ctrl_set_pack_xyz(
439         &         cunit, ivartype, fname_uvel(ictrlgrad), "maskCtrlW",
440       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, wunit, lxxadxx, mythid)
441  #endif  #endif
442    
443    #ifdef ALLOW_VVEL0_CONTROL
444              ivartype = 28
445              write(weighttype(1:80),'(80a)') ' '
446              write(weighttype(1:80),'(a)') "wvvel"
447              call ctrl_set_pack_xyz(
448         &         cunit, ivartype, fname_vvel(ictrlgrad), "maskCtrlS",
449         &         weighttype, wunit, lxxadxx, mythid)
450    #endif
451    
452    #ifdef ALLOW_ETAN0_CONTROL
453              ivartype = 29
454              write(weighttype(1:80),'(80a)') ' '
455              write(weighttype(1:80),'(a)') "wetan"
456              call ctrl_set_pack_xy(
457         &         cunit, ivartype, fname_etan(ictrlgrad), "maskCtrlC",
458         &         weighttype, lxxadxx, mythid)
459    #endif
460    
461    #ifdef ALLOW_RELAXSST_CONTROL
462              ivartype = 30
463              write(weighttype(1:80),'(80a)') ' '
464              write(weighttype(1:80),'(a)') "wrelaxsst"
465              call ctrl_set_pack_xy(
466         &         cunit, ivartype, fname_relaxsst(ictrlgrad), "maskCtrlC",
467         &         weighttype, lxxadxx, mythid)
468    #endif
469    
470    #ifdef ALLOW_RELAXSSS_CONTROL
471              ivartype = 31
472              write(weighttype(1:80),'(80a)') ' '
473              write(weighttype(1:80),'(a)') "wrelaxsss"
474              call ctrl_set_pack_xy(
475         &         cunit, ivartype, fname_relaxsss(ictrlgrad), "maskCtrlC",
476         &         weighttype, lxxadxx, mythid)
477    #endif
478    
479    #ifdef ALLOW_PRECIP_CONTROL
480              ivartype = 32
481              write(weighttype(1:80),'(80a)') ' '
482              write(weighttype(1:80),'(a)') "wprecip"
483              call ctrl_set_pack_xy(
484         &         cunit, ivartype, fname_precip(ictrlgrad), "maskCtrlC",
485         &         weighttype, lxxadxx, mythid)
486    #endif
487    
488    #ifdef ALLOW_SWFLUX_CONTROL
489              ivartype = 33
490              write(weighttype(1:80),'(80a)') ' '
491              write(weighttype(1:80),'(a)') "wswflux"
492              call ctrl_set_pack_xy(
493         &         cunit, ivartype, fname_swflux(ictrlgrad), "maskCtrlC",
494         &         weighttype, lxxadxx, mythid)
495    #endif
496    
497    #ifdef ALLOW_SWDOWN_CONTROL
498              ivartype = 34
499              write(weighttype(1:80),'(80a)') ' '
500              write(weighttype(1:80),'(a)') "wswdown"
501              call ctrl_set_pack_xy(
502         &         cunit, ivartype, fname_swdown(ictrlgrad), "maskCtrlC",
503         &         weighttype, lxxadxx, mythid)
504    #endif
505    
506            close ( cunit )            close ( cunit )
507    
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.22