/[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.12 by heimbach, Thu Nov 6 22:05:08 2003 UTC revision 1.21 by heimbach, Thu Jul 28 13:47:49 2005 UTC
# Line 1  Line 1 
 C  
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
# Line 43  c     == global variables == Line 42  c     == global variables ==
42  #include "GRID.h"  #include "GRID.h"
43    
44  #include "ctrl.h"  #include "ctrl.h"
45  #include "cost.h"  #include "optim.h"
46    
47    #ifdef ALLOW_COST
48    # include "cost.h"
49    #endif
50  #ifdef ALLOW_ECCO  #ifdef ALLOW_ECCO
51  # include "ecco_cost.h"  # include "ecco_cost.h"
52  #else  #else
53  # include "ctrl_weights.h"  # include "ctrl_weights.h"
54  #endif  #endif
55    
 #ifdef ALLOW_ECCO_OPTIMIZATION  
 # include "optim.h"  
 #endif  
   
56  c     == routine arguments ==  c     == routine arguments ==
57    
58        logical first        logical first
# Line 63  c     == routine arguments == Line 61  c     == routine arguments ==
61  #ifndef EXCLUDE_CTRL_PACK  #ifndef EXCLUDE_CTRL_PACK
62  c     == local variables ==  c     == local variables ==
63    
 #ifndef ALLOW_ECCO_OPTIMIZATION  
       integer optimcycle  
       _RL    fmin  
 #endif  
   
64        _RL    fcloc        _RL    fcloc
65    
66        integer i, j, k        integer i, j, k
# Line 97  c     == external == Line 90  c     == external ==
90  c     == end of interface ==  c     == end of interface ==
91    
92  #ifndef ALLOW_ECCO_OPTIMIZATION  #ifndef ALLOW_ECCO_OPTIMIZATION
       optimcycle = 0  
93        fmin       = 0. _d 0        fmin       = 0. _d 0
94  #endif  #endif
95    
# Line 107  c--   Tiled files are used. Line 99  c--   Tiled files are used.
99  c--   Initialise adjoint variables on active files.  c--   Initialise adjoint variables on active files.
100        ladinit = .false.        ladinit = .false.
101    
102    c--   Initialise global buffer index
103          nbuffglobal = 0
104    
105  c--   Assign file names.  c--   Assign file names.
106    
107        call ctrl_set_fname(xx_theta_file, fname_theta, mythid)        call ctrl_set_fname(xx_theta_file, fname_theta, mythid)
# Line 117  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_uwind_file, fname_uwind, mythid)        call ctrl_set_fname(xx_uwind_file, fname_uwind, mythid)
117        call ctrl_set_fname(xx_vwind_file, fname_vwind, mythid)        call ctrl_set_fname(xx_vwind_file, fname_vwind, mythid)
118        call ctrl_set_fname(xx_obcsn_file, fname_obcsn, mythid)        call ctrl_set_fname(xx_obcsn_file, fname_obcsn, mythid)
# Line 132  c--   Assign file names. Line 128  c--   Assign file names.
128        call ctrl_set_fname(xx_efluxy_file, fname_efluxy, mythid)        call ctrl_set_fname(xx_efluxy_file, fname_efluxy, mythid)
129        call ctrl_set_fname(xx_efluxp_file, fname_efluxp, mythid)        call ctrl_set_fname(xx_efluxp_file, fname_efluxp, mythid)
130        call ctrl_set_fname(xx_bottomdrag_file, fname_bottomdrag, mythid)        call ctrl_set_fname(xx_bottomdrag_file, fname_bottomdrag, mythid)
131          call ctrl_set_fname(xx_edtaux_file, fname_edtaux, mythid)
132          call ctrl_set_fname(xx_edtauy_file, fname_edtauy, mythid)
133          call ctrl_set_fname(xx_uvel_file, fname_uvel, mythid)
134          call ctrl_set_fname(xx_vvel_file, fname_vvel, mythid)
135          call ctrl_set_fname(xx_etan_file, fname_etan, mythid)
136          call ctrl_set_fname(xx_relaxsst_file, fname_relaxsst, mythid)
137          call ctrl_set_fname(xx_relaxsss_file, fname_relaxsss, mythid)
138    
139  c  c--   Only the master thread will do I/O.
 c--     Only the master thread will do I/O.  
140        _BEGIN_MASTER( mythid )        _BEGIN_MASTER( mythid )
141    
142        if ( first .AND. optimcycle .EQ. 0 ) then        if ( first ) then
143  c     >>> Initialise control vector for optimcycle=0 <<<  c     >>> Initialise control vector for optimcycle=0 <<<
144            lxxadxx   = .TRUE.            lxxadxx   = .TRUE.
145            ictrlgrad = 1            ictrlgrad = 1
146            fcloc     = fmin            fcloc     = fmin
147            write(cfile(1:128),'(4a,i4.4)')            write(cfile(1:128),'(4a,i4.4)')
148       &      ctrlname(1:9),'_',yctrlid(1:10),'.opt', optimcycle       &         ctrlname(1:9),'_',yctrlid(1:10),
149         &         yctrlpospack, optimcycle
150              print *, 'ph-pack: packing ', ctrlname(1:9)
151        else        else
152  c     >>> Write gradient vector <<<  c     >>> Write gradient vector <<<
153            lxxadxx   = .FALSE.            lxxadxx   = .FALSE.
154            ictrlgrad = 2            ictrlgrad = 2
155            fcloc     = fc            fcloc     = fc
156            write(cfile(1:128),'(4a,i4.4)')            write(cfile(1:128),'(4a,i4.4)')
157       &    costname(1:9),'_',yctrlid(1:10),'.opt', optimcycle       &         costname(1:9),'_',yctrlid(1:10),
158         &         yctrlpospack, optimcycle
159              print *, 'ph-pack: packing ', costname(1:9)
160         endif         endif
161    
162         call mdsfindunit( cunit, mythid )         call mdsfindunit( cunit, mythid )
# Line 164  c--       Header information. Line 170  c--       Header information.
170            write(cunit) nvarlength            write(cunit) nvarlength
171            write(cunit) yctrlid            write(cunit) yctrlid
172            write(cunit) optimCycle            write(cunit) optimCycle
173            write(cunit) fcloc            write(cunit) fc
174            write(cunit) 1  C     place holder of obsolete variable iG
           write(cunit) 1  
175            write(cunit) 1            write(cunit) 1
176    C     place holder of obsolete variable jG
177            write(cunit) 1            write(cunit) 1
178              write(cunit) nsx
179              write(cunit) nsy
180            write(cunit) (nWetcGlobal(k), k=1,nr)            write(cunit) (nWetcGlobal(k), k=1,nr)
181            write(cunit) (nWetsGlobal(k), k=1,nr)            write(cunit) (nWetsGlobal(k), k=1,nr)
182            write(cunit) (nWetwGlobal(k), k=1,nr)            write(cunit) (nWetwGlobal(k), k=1,nr)
# Line 190  c--       Header information. Line 198  c--       Header information.
198  #endif  #endif
199            write(cunit) (ncvarindex(i), i=1,maxcvars)            write(cunit) (ncvarindex(i), i=1,maxcvars)
200            write(cunit) (ncvarrecs(i),  i=1,maxcvars)            write(cunit) (ncvarrecs(i),  i=1,maxcvars)
201            write(cunit) (nx,  i=1,maxcvars)            write(cunit) (ncvarxmax(i),  i=1,maxcvars)
202            write(cunit) (ny,  i=1,maxcvars)            write(cunit) (ncvarymax(i),  i=1,maxcvars)
203            write(cunit) (ncvarnrmax(i), i=1,maxcvars)            write(cunit) (ncvarnrmax(i), i=1,maxcvars)
204            write(cunit) (ncvargrd(i),   i=1,maxcvars)            write(cunit) (ncvargrd(i),   i=1,maxcvars)
205            write(cunit)            write(cunit)
# Line 201  c--       Header information. Line 209  c--       Header information.
209            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
210            write(weighttype(1:80),'(a)') "wtheta"            write(weighttype(1:80),'(a)') "wtheta"
211            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
212       &         cunit, ivartype, fname_theta(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_theta(ictrlgrad), "maskCtrlC",
213       &         weighttype, wtheta, lxxadxx, mythid)       &         weighttype, wtheta, lxxadxx, mythid)
214  #endif  #endif
215    
# Line 210  c--       Header information. Line 218  c--       Header information.
218            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
219            write(weighttype(1:80),'(a)') "wsalt"            write(weighttype(1:80),'(a)') "wsalt"
220            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
221       &         cunit, ivartype, fname_salt(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_salt(ictrlgrad), "maskCtrlC",
222       &         weighttype, wsalt, lxxadxx, mythid)       &         weighttype, wsalt, lxxadxx, mythid)
223  #endif  #endif
224    
# Line 220  c--       Header information. Line 228  c--       Header information.
228            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
229            write(weighttype(1:80),'(a)') "whflux"            write(weighttype(1:80),'(a)') "whflux"
230            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
231       &         cunit, ivartype, fname_hflux(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_hflux(ictrlgrad), "maskCtrlC",
232       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
233  #endif  #endif
234    
# Line 230  c--       Header information. Line 238  c--       Header information.
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    
# Line 240  c--       Header information. Line 248  c--       Header information.
248            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
249            write(weighttype(1:80),'(a)') "wtauu"            write(weighttype(1:80),'(a)') "wtauu"
250            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
251       &         cunit, ivartype, fname_tauu(ictrlgrad), "maskW",       &         cunit, ivartype, fname_tauu(ictrlgrad), "maskCtrlW",
252       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
253  #endif  #endif
254    
# Line 250  c--       Header information. Line 258  c--       Header information.
258            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
259            write(weighttype(1:80),'(a)') "wtauv"            write(weighttype(1:80),'(a)') "wtauv"
260            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
261       &         cunit, ivartype, fname_tauv(ictrlgrad), "maskS",       &         cunit, ivartype, fname_tauv(ictrlgrad), "maskCtrlS",
262       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
263  #endif  #endif
264    
# Line 259  c--       Header information. Line 267  c--       Header information.
267            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
268            write(weighttype(1:80),'(a)') "watemp"            write(weighttype(1:80),'(a)') "watemp"
269            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
270       &         cunit, ivartype, fname_atemp(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_atemp(ictrlgrad), "maskCtrlC",
271       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
272  #endif  #endif
273    
# Line 268  c--       Header information. Line 276  c--       Header information.
276            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
277            write(weighttype(1:80),'(a)') "waqh"            write(weighttype(1:80),'(a)') "waqh"
278            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
279       &         cunit, ivartype, fname_aqh(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_aqh(ictrlgrad), "maskCtrlC",
280       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
281  #endif  #endif
282    
# Line 277  c--       Header information. Line 285  c--       Header information.
285            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
286            write(weighttype(1:80),'(a)') "wuwind"            write(weighttype(1:80),'(a)') "wuwind"
287            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
288       &         cunit, ivartype, fname_uwind(ictrlgrad), "maskW",       &         cunit, ivartype, fname_uwind(ictrlgrad), "maskCtrlC",
289       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
290  #endif  #endif
291    
# Line 286  c--       Header information. Line 294  c--       Header information.
294            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
295            write(weighttype(1:80),'(a)') "wvwind"            write(weighttype(1:80),'(a)') "wvwind"
296            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
297       &         cunit, ivartype, fname_vwind(ictrlgrad), "maskS",       &         cunit, ivartype, fname_vwind(ictrlgrad), "maskCtrlC",
298       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
299  #endif  #endif
300    
# Line 331  c--       Header information. Line 339  c--       Header information.
339            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
340            write(weighttype(1:80),'(a)') "wdiffkr"            write(weighttype(1:80),'(a)') "wdiffkr"
341            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
342       &         cunit, ivartype, fname_diffkr(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_diffkr(ictrlgrad), "maskCtrlC",
343       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, wunit, lxxadxx, mythid)
344  #endif  #endif
345    
# Line 340  c--       Header information. Line 348  c--       Header information.
348            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
349            write(weighttype(1:80),'(a)') "wkapgm"            write(weighttype(1:80),'(a)') "wkapgm"
350            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
351       &         cunit, ivartype, fname_kapgm(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_kapgm(ictrlgrad), "maskCtrlC",
352       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, wunit, lxxadxx, mythid)
353  #endif  #endif
354    
355  #ifdef ALLOW_TR10_CONTROL  #ifdef ALLOW_PRECIP_CONTROL
356            ivartype = 17            ivartype = 17
357            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
358            write(weighttype(1:80),'(a)') "wtr1"            write(weighttype(1:80),'(a)') "wprecip"
359            call ctrl_set_pack_xyz(            call ctrl_set_pack_xy(
360       &         cunit, ivartype, fname_tr1(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_precip(ictrlgrad), "maskCtrlC",
361       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
362  #endif  #endif
363    
364  #ifdef ALLOW_SST0_CONTROL  #ifdef ALLOW_SST0_CONTROL
# Line 358  c--       Header information. Line 366  c--       Header information.
366            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
367            write(weighttype(1:80),'(a)') "wsst0"            write(weighttype(1:80),'(a)') "wsst0"
368            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
369       &         cunit, ivartype, fname_sst(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_sst(ictrlgrad), "maskCtrlC",
370       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
371  #endif  #endif
372    
# Line 367  c--       Header information. Line 375  c--       Header information.
375            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
376            write(weighttype(1:80),'(a)') "wsss0"            write(weighttype(1:80),'(a)') "wsss0"
377            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
378       &         cunit, ivartype, fname_sss(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_sss(ictrlgrad), "maskCtrlC",
379       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
380  #endif  #endif
381    
# Line 377  c--       Header information. Line 385  c--       Header information.
385            write(weighttype(1:80),'(a)') "whfacc"            write(weighttype(1:80),'(a)') "whfacc"
386  # ifdef ALLOW_HFACC3D_CONTROL  # ifdef ALLOW_HFACC3D_CONTROL
387            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
388       &         cunit, ivartype, fname_hfacc(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_hfacc(ictrlgrad), "maskCtrlC",
389       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, wunit, lxxadxx, mythid)
390  # else  # else
391            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
392       &         cunit, ivartype, fname_hfacc(ictrlgrad), "hFacC",       &         cunit, ivartype, fname_hfacc(ictrlgrad), "maskCtrlC",
393       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
394  # endif  # endif
395  #endif  #endif
# Line 391  c--       Header information. Line 399  c--       Header information.
399            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
400            write(weighttype(1:80),'(a)') "wefluxy0"            write(weighttype(1:80),'(a)') "wefluxy0"
401            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
402       &         cunit, ivartype, fname_efluxy(ictrlgrad), "hFacS",       &         cunit, ivartype, fname_efluxy(ictrlgrad), "maskCtrlS",
403       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, wunit, lxxadxx, mythid)
404  #endif  #endif
405    
# Line 400  c--       Header information. Line 408  c--       Header information.
408            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
409            write(weighttype(1:80),'(a)') "wefluxp0"            write(weighttype(1:80),'(a)') "wefluxp0"
410            call ctrl_set_pack_xyz(            call ctrl_set_pack_xyz(
411       &         cunit, ivartype, fname_efluxp(ictrlgrad), "hFacV",       &         cunit, ivartype, fname_efluxp(ictrlgrad), "maskhFacV",
412       &         weighttype, wunit, lxxadxx, mythid)       &         weighttype, wunit, lxxadxx, mythid)
413  #endif  #endif
414    
# Line 409  c--       Header information. Line 417  c--       Header information.
417            write(weighttype(1:80),'(80a)') ' '            write(weighttype(1:80),'(80a)') ' '
418            write(weighttype(1:80),'(a)') "wbottomdrag"            write(weighttype(1:80),'(a)') "wbottomdrag"
419            call ctrl_set_pack_xy(            call ctrl_set_pack_xy(
420       &         cunit, ivartype, fname_bottomdrag(ictrlgrad), "hFacC",       &      cunit, ivartype, fname_bottomdrag(ictrlgrad), "maskCtrlC",
421         &      weighttype, lxxadxx, mythid)
422    #endif
423    
424    #ifdef ALLOW_EDTAUX_CONTROL
425              ivartype = 25
426              write(weighttype(1:80),'(80a)') ' '
427              write(weighttype(1:80),'(a)') "wedtaux"
428              call ctrl_set_pack_xyz(
429         &         cunit, ivartype, fname_edtaux(ictrlgrad), "maskCtrlW",
430         &         weighttype, wunit, lxxadxx, mythid)
431    #endif
432    
433    #ifdef ALLOW_EDTAUY_CONTROL
434              ivartype = 26
435              write(weighttype(1:80),'(80a)') ' '
436              write(weighttype(1:80),'(a)') "wedtauy"
437              call ctrl_set_pack_xyz(
438         &         cunit, ivartype, fname_edtauy(ictrlgrad), "maskCtrlS",
439         &         weighttype, wunit, lxxadxx, mythid)
440    #endif
441    
442    #ifdef ALLOW_UVEL0_CONTROL
443              ivartype = 27
444              write(weighttype(1:80),'(80a)') ' '
445              write(weighttype(1:80),'(a)') "wuvel"
446              call ctrl_set_pack_xyz(
447         &         cunit, ivartype, fname_uvel(ictrlgrad), "maskCtrlW",
448         &         weighttype, wunit, lxxadxx, mythid)
449    #endif
450    
451    #ifdef ALLOW_VVEL0_CONTROL
452              ivartype = 28
453              write(weighttype(1:80),'(80a)') ' '
454              write(weighttype(1:80),'(a)') "wvvel"
455              call ctrl_set_pack_xyz(
456         &         cunit, ivartype, fname_vvel(ictrlgrad), "maskCtrlS",
457         &         weighttype, wunit, lxxadxx, mythid)
458    #endif
459    
460    #ifdef ALLOW_ETAN0_CONTROL
461              ivartype = 29
462              write(weighttype(1:80),'(80a)') ' '
463              write(weighttype(1:80),'(a)') "wetan"
464              call ctrl_set_pack_xy(
465         &         cunit, ivartype, fname_etan(ictrlgrad), "maskCtrlC",
466         &         weighttype, lxxadxx, mythid)
467    #endif
468    
469    #ifdef ALLOW_RELAXSST_CONTROL
470              ivartype = 30
471              write(weighttype(1:80),'(80a)') ' '
472              write(weighttype(1:80),'(a)') "wrelaxsst"
473              call ctrl_set_pack_xy(
474         &         cunit, ivartype, fname_relaxsst(ictrlgrad), "maskCtrlC",
475         &         weighttype, lxxadxx, mythid)
476    #endif
477    
478    #ifdef ALLOW_RELAXSSS_CONTROL
479              ivartype = 31
480              write(weighttype(1:80),'(80a)') ' '
481              write(weighttype(1:80),'(a)') "wrelaxsss"
482              call ctrl_set_pack_xy(
483         &         cunit, ivartype, fname_relaxsss(ictrlgrad), "maskCtrlC",
484       &         weighttype, lxxadxx, mythid)       &         weighttype, lxxadxx, mythid)
485  #endif  #endif
486    
487    #ifdef ALLOW_TR10_CONTROL
488              ivartype = 32
489              write(weighttype(1:80),'(80a)') ' '
490              write(weighttype(1:80),'(a)') "wtr1"
491              call ctrl_set_pack_xyz(
492         &         cunit, ivartype, fname_tr1(ictrlgrad), "maskCtrlC",
493         &         weighttype, wunit, lxxadxx, mythid)
494    #endif
495    
496            close ( cunit )            close ( cunit )
497    
498          _END_MASTER( mythid )          _END_MASTER( mythid )
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22