/[MITgcm]/MITgcm/pkg/autodiff/active_file_control_slice.F
ViewVC logotype

Diff of /MITgcm/pkg/autodiff/active_file_control_slice.F

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

revision 1.8 by heimbach, Tue Jul 26 13:10:46 2005 UTC revision 1.9 by jmc, Mon Oct 8 23:50:53 2007 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "AUTODIFF_OPTIONS.h"  #include "AUTODIFF_OPTIONS.h"
5    
# Line 115  c     >>>>>>>>>>>>>>>>>>>             << Line 117  c     >>>>>>>>>>>>>>>>>>>             <<
117    
118  c       Read the active variable from file.  c       Read the active variable from file.
119    
120          call mdsreadfieldxz(          call mdsreadfieldxz(
121       &                     active_var_file,       &                     active_var_file,
122       &                     prec,       &                     prec,
123       &                     'RL',       &                     'RL',
124       &                     mynr,             &                     mynr,
125       &                     active_var,       &                     active_var,
126       &                     irec,       &                     irec,
127       &                     mythid )       &                     mythid )
128    
129          if (lAdInit) then          if (lAdInit) then
# Line 139  c         adjoint variable's file. These Line 141  c         adjoint variable's file. These
141               enddo               enddo
142            enddo            enddo
143    
144            call mdswritefieldxz(            call mdswritefieldxz(
145       &                           adfname,       &                           adfname,
146       &                           prec,       &                           prec,
147       &                           globalfile,       &                           globalfile,
148       &                           'RL',       &                           'RL',
149       &                           mynr,       &                           mynr,
150       &                           active_data_t,       &                           active_data_t,
151       &                           irec,       &                           irec,
152       &                           myOptimIter,       &                           myOptimIter,
# Line 164  c     >>>>>>>>>>>>>>>>>>>             << Line 166  c     >>>>>>>>>>>>>>>>>>>             <<
166          _BEGIN_MASTER( mythid )          _BEGIN_MASTER( mythid )
167    
168          writeglobalfile = .false.          writeglobalfile = .false.
169          call mdsreadfieldxz(          call mdsreadfieldxz(
170       &                        active_var_file,       &                        active_var_file,
171       &                        prec,       &                        prec,
172       &                        'RL',       &                        'RL',
173       &                        mynr,       &                        mynr,
174       &                        active_data_t,       &                        active_data_t,
175       &                        irec,       &                        irec,
176       &                        mythid )       &                        mythid )
177    
178  c             Add active_var from appropriate location to data.  c             Add active_var from appropriate location to data.
# Line 178  c             Add active_var from approp Line 180  c             Add active_var from approp
180             do bi = 1,nsx             do bi = 1,nsx
181                do k = 1,mynr                do k = 1,mynr
182                   do i = 1,snx                   do i = 1,snx
183                      active_data_t(i,k,bi,bj) =                      active_data_t(i,k,bi,bj) =
184       &                   active_data_t(i,k,bi,bj) +       &                   active_data_t(i,k,bi,bj) +
185       &                   active_var(i,k,bi,bj)       &                   active_var(i,k,bi,bj)
186                   enddo                   enddo
187                enddo                enddo
188             enddo             enddo
189          enddo          enddo
190            
191  c             Store the result on disk.  c             Store the result on disk.
192          call mdswritefieldxz(          call mdswritefieldxz(
193       &                         active_var_file,       &                         active_var_file,
194       &                         prec,       &                         prec,
195       &                         writeglobalfile,       &                         writeglobalfile,
196       &                         'RL',       &                         'RL',
197       &                         mynr,       &                         mynr,
198       &                         active_data_t,       &                         active_data_t,
199       &                         irec,       &                         irec,
200       &                         myOptimIter,       &                         myOptimIter,
201       &                         mythid )       &                         mythid )
202    
# Line 222  c     >>>>>>>>>>>>>>>>>>>             << Line 224  c     >>>>>>>>>>>>>>>>>>>             <<
224    
225  c       Read the active variable from file.  c       Read the active variable from file.
226    
227          call mdsreadfieldxz(          call mdsreadfieldxz(
228       &                     active_var_file,       &                     active_var_file,
229       &                     prec,       &                     prec,
230       &                     'RL',       &                     'RL',
231       &                     mynr,             &                     mynr,
232       &                     active_var,       &                     active_var,
233       &                     irec,       &                     irec,
234       &                     mythid )       &                     mythid )
235    
236          _END_MASTER( mythid )          _END_MASTER( mythid )
# Line 312  c     >>>>>>>>>>>>>>>>>>>             << Line 314  c     >>>>>>>>>>>>>>>>>>>             <<
314    
315          _BEGIN_MASTER( mythid )          _BEGIN_MASTER( mythid )
316    
317            call mdswritefieldxz(            call mdswritefieldxz(
318       &                        active_var_file,       &                        active_var_file,
319       &                        prec,       &                        prec,
320       &                        globalfile,       &                        globalfile,
321       &                        'RL',       &                        'RL',
322       &                        mynr,       &                        mynr,
323       &                        active_var,       &                        active_var,
324       &                        irec,       &                        irec,
325       &                        myOptimIter,       &                        myOptimIter,
326       &                        mythid )       &                        mythid )
327    
# Line 335  c     >>>>>>>>>>>>>>>>>>>             << Line 337  c     >>>>>>>>>>>>>>>>>>>             <<
337    
338          _BEGIN_MASTER( mythid )          _BEGIN_MASTER( mythid )
339    
340          call mdsreadfieldxz(          call mdsreadfieldxz(
341       &                           active_var_file,       &                           active_var_file,
342       &                           prec,       &                           prec,
343       &                           'RL',       &                           'RL',
344       &                            mynr,       &                            mynr,
345       &                            active_data_t,       &                            active_data_t,
346       &                            irec,       &                            irec,
347       &                            mythid )       &                            mythid )
348    
# Line 349  c             Add active_var from approp Line 351  c             Add active_var from approp
351             do bi = 1,nsx             do bi = 1,nsx
352                do k = 1,nr                do k = 1,nr
353                   do i = 1,snx                   do i = 1,snx
354                      active_var(i,k,bi,bj) =                      active_var(i,k,bi,bj) =
355       &                   active_var(i,k,bi,bj) +       &                   active_var(i,k,bi,bj) +
356       &                   active_data_t(i,k,bi,bj)       &                   active_data_t(i,k,bi,bj)
357                      active_data_t(i,k,bi,bj) = 0. _d 0                      active_data_t(i,k,bi,bj) = 0. _d 0
# Line 357  c             Add active_var from approp Line 359  c             Add active_var from approp
359                enddo                enddo
360             enddo             enddo
361          enddo          enddo
362          call mdswritefieldxz(          call mdswritefieldxz(
363       &                            active_var_file,       &                            active_var_file,
364       &                            prec,       &                            prec,
365       &                            globalfile,       &                            globalfile,
366       &                            'RL',       &                            'RL',
367       &                            mynr,       &                            mynr,
368       &                            active_data_t,       &                            active_data_t,
369       &                            irec,       &                            irec,
370       &                            myOptimIter,       &                            myOptimIter,
371       &                            mythid )       &                            mythid )
372    
# Line 380  c     >>>>>>>>>>>>>>>>>>>             << Line 382  c     >>>>>>>>>>>>>>>>>>>             <<
382    
383          _BEGIN_MASTER( mythid )          _BEGIN_MASTER( mythid )
384    
385            call mdswritefieldxz(            call mdswritefieldxz(
386       &                        active_var_file,       &                        active_var_file,
387       &                        prec,       &                        prec,
388       &                        globalfile,       &                        globalfile,
389       &                        'RL',       &                        'RL',
390       &                        mynr,       &                        mynr,
391       &                        active_var,       &                        active_var,
392       &                        irec,       &                        irec,
393       &                        myOptimIter,       &                        myOptimIter,
394       &                        mythid )       &                        mythid )
395    
# Line 499  c     >>>>>>>>>>>>>>>>>>>             << Line 501  c     >>>>>>>>>>>>>>>>>>>             <<
501    
502  c       Read the active variable from file.  c       Read the active variable from file.
503    
504          call mdsreadfieldyz(          call mdsreadfieldyz(
505       &                     active_var_file,       &                     active_var_file,
506       &                     prec,       &                     prec,
507       &                     'RL',       &                     'RL',
508       &                     mynr,             &                     mynr,
509       &                     active_var,       &                     active_var,
510       &                     irec,       &                     irec,
511       &                     mythid )       &                     mythid )
512    
513          if (lAdInit) then          if (lAdInit) then
# Line 523  c         adjoint variable's file. These Line 525  c         adjoint variable's file. These
525               enddo               enddo
526            enddo            enddo
527    
528            call mdswritefieldyz(            call mdswritefieldyz(
529       &                           adfname,       &                           adfname,
530       &                           prec,       &                           prec,
531       &                           globalfile,       &                           globalfile,
532       &                           'RL',       &                           'RL',
533       &                           mynr,       &                           mynr,
534       &                           active_data_t,       &                           active_data_t,
535       &                           irec,       &                           irec,
536       &                           myOptimIter,       &                           myOptimIter,
# Line 548  c     >>>>>>>>>>>>>>>>>>>             << Line 550  c     >>>>>>>>>>>>>>>>>>>             <<
550          _BEGIN_MASTER( mythid )          _BEGIN_MASTER( mythid )
551    
552          writeglobalfile = .false.          writeglobalfile = .false.
553          call mdsreadfieldyz(          call mdsreadfieldyz(
554       &                        active_var_file,       &                        active_var_file,
555       &                        prec,       &                        prec,
556       &                        'RL',       &                        'RL',
557       &                        mynr,       &                        mynr,
558       &                        active_data_t,       &                        active_data_t,
559       &                        irec,       &                        irec,
560       &                        mythid )       &                        mythid )
561    
562  c             Add active_var from appropriate location to data.  c             Add active_var from appropriate location to data.
# Line 562  c             Add active_var from approp Line 564  c             Add active_var from approp
564             do bi = 1,nsx             do bi = 1,nsx
565                do k = 1,mynr                do k = 1,mynr
566                   do j = 1,sny                   do j = 1,sny
567                      active_data_t(j,k,bi,bj) =                      active_data_t(j,k,bi,bj) =
568       &                   active_data_t(j,k,bi,bj) +       &                   active_data_t(j,k,bi,bj) +
569       &                   active_var(j,k,bi,bj)       &                   active_var(j,k,bi,bj)
570                   enddo                   enddo
571                enddo                enddo
572             enddo             enddo
573          enddo          enddo
574            
575  c             Store the result on disk.  c             Store the result on disk.
576          call mdswritefieldyz(          call mdswritefieldyz(
577       &                         active_var_file,       &                         active_var_file,
578       &                         prec,       &                         prec,
579       &                         writeglobalfile,       &                         writeglobalfile,
580       &                         'RL',       &                         'RL',
581       &                         mynr,       &                         mynr,
582       &                         active_data_t,       &                         active_data_t,
583       &                         irec,       &                         irec,
584       &                         myOptimIter,       &                         myOptimIter,
585       &                         mythid )       &                         mythid )
586    
# Line 606  c     >>>>>>>>>>>>>>>>>>>             << Line 608  c     >>>>>>>>>>>>>>>>>>>             <<
608    
609  c       Read the active variable from file.  c       Read the active variable from file.
610    
611          call mdsreadfieldyz(          call mdsreadfieldyz(
612       &                     active_var_file,       &                     active_var_file,
613       &                     prec,       &                     prec,
614       &                     'RL',       &                     'RL',
615       &                     mynr,             &                     mynr,
616       &                     active_var,       &                     active_var,
617       &                     irec,       &                     irec,
618       &                     mythid )       &                     mythid )
619    
620          _END_MASTER( mythid )          _END_MASTER( mythid )
# Line 696  c     >>>>>>>>>>>>>>>>>>>             << Line 698  c     >>>>>>>>>>>>>>>>>>>             <<
698    
699          _BEGIN_MASTER( mythid )          _BEGIN_MASTER( mythid )
700    
701            call mdswritefieldyz(            call mdswritefieldyz(
702       &                        active_var_file,       &                        active_var_file,
703       &                        prec,       &                        prec,
704       &                        globalfile,       &                        globalfile,
705       &                        'RL',       &                        'RL',
706       &                        mynr,       &                        mynr,
707       &                        active_var,       &                        active_var,
708       &                        irec,       &                        irec,
709       &                        myOptimIter,       &                        myOptimIter,
710       &                        mythid )       &                        mythid )
711    
# Line 719  c     >>>>>>>>>>>>>>>>>>>             << Line 721  c     >>>>>>>>>>>>>>>>>>>             <<
721    
722          _BEGIN_MASTER( mythid )          _BEGIN_MASTER( mythid )
723    
724          call mdsreadfieldyz(          call mdsreadfieldyz(
725       &                           active_var_file,       &                           active_var_file,
726       &                           prec,       &                           prec,
727       &                           'RL',       &                           'RL',
728       &                            mynr,       &                            mynr,
729       &                            active_data_t,       &                            active_data_t,
730       &                            irec,       &                            irec,
731       &                            mythid )       &                            mythid )
732    
733  c             Add active_var from appropriate location to data.  c             Add active_var from appropriate location to data.
# Line 733  c             Add active_var from approp Line 735  c             Add active_var from approp
735             do bi = 1,nsx             do bi = 1,nsx
736                do k = 1, mynr                do k = 1, mynr
737                   do j = 1,sny                   do j = 1,sny
738                      active_var(j,k,bi,bj) =                      active_var(j,k,bi,bj) =
739       &                   active_var(j,k,bi,bj) +       &                   active_var(j,k,bi,bj) +
740       &                   active_data_t(j,k,bi,bj)       &                   active_data_t(j,k,bi,bj)
741                      active_data_t(j,k,bi,bj) = 0. _d 0                      active_data_t(j,k,bi,bj) = 0. _d 0
# Line 741  c             Add active_var from approp Line 743  c             Add active_var from approp
743                enddo                enddo
744             enddo             enddo
745          enddo          enddo
746          call mdswritefieldyz(          call mdswritefieldyz(
747       &                            active_var_file,       &                            active_var_file,
748       &                            prec,       &                            prec,
749       &                            globalfile,       &                            globalfile,
750       &                            'RL',       &                            'RL',
751       &                            mynr,       &                            mynr,
752       &                            active_data_t,       &                            active_data_t,
753       &                            irec,       &                            irec,
754       &                            myOptimIter,       &                            myOptimIter,
755       &                            mythid )       &                            mythid )
756    
# Line 764  c     >>>>>>>>>>>>>>>>>>>             << Line 766  c     >>>>>>>>>>>>>>>>>>>             <<
766    
767          _BEGIN_MASTER( mythid )          _BEGIN_MASTER( mythid )
768    
769            call mdswritefieldyz(            call mdswritefieldyz(
770       &                        active_var_file,       &                        active_var_file,
771       &                        prec,       &                        prec,
772       &                        globalfile,       &                        globalfile,
773       &                        'RL',       &                        'RL',
774       &                        mynr,       &                        mynr,
775       &                        active_var,       &                        active_var,
776       &                        irec,       &                        irec,
777       &                        myOptimIter,       &                        myOptimIter,
778       &                        mythid )       &                        mythid )
779    
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22