/[MITgcm]/MITgcm/pkg/admtlm/admtlm_bypassad.F
ViewVC logotype

Diff of /MITgcm/pkg/admtlm/admtlm_bypassad.F

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

revision 1.2 by heimbach, Thu Apr 27 12:49:02 2006 UTC revision 1.3 by jmc, Mon Oct 8 23:45:51 2007 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
# Line 58  C     ================================== Line 60  C     ==================================
60        jmax = sny        jmax = sny
61        imin = 1        imin = 1
62        imax = snx        imax = snx
63        doglobalread =  .false.        doglobalread =  .false.
64        ladinit =  .false.        ladinit =  .false.
65        equal =  .true.        equal =  .true.
66        if (equal) then        if (equal) then
67          fac = 1.d0          fac = 1.d0
68        else        else
# Line 72  C     ================================== Line 74  C     ==================================
74          DO j = jmin, jmax          DO j = jmin, jmax
75           DO i = imin, imax           DO i = imin, imax
76            DO k=1,Nr            DO k=1,Nr
77              adtheta(i,j,k,bi,bj) =              adtheta(i,j,k,bi,bj) =
78       &       g_objf_state_final(i,j,bi,bj,k)       &       g_objf_state_final(i,j,bi,bj,k)
79              adsalt(i,j,k,bi,bj) =              adsalt(i,j,k,bi,bj) =
80       &       g_objf_state_final(i,j,bi,bj,1*Nr+k)       &       g_objf_state_final(i,j,bi,bj,1*Nr+k)
81              aduvel(i,j,k,bi,bj) =              aduvel(i,j,k,bi,bj) =
82       &       g_objf_state_final(i,j,bi,bj,2*Nr+k)       &       g_objf_state_final(i,j,bi,bj,2*Nr+k)
83              advvel(i,j,k,bi,bj) =              advvel(i,j,k,bi,bj) =
84       &       g_objf_state_final(i,j,bi,bj,3*Nr+k)       &       g_objf_state_final(i,j,bi,bj,3*Nr+k)
85            END DO            END DO
86            adetan(i,j,bi,bj) =            adetan(i,j,bi,bj) =
87       &       g_objf_state_final(i,j,bi,bj,4*Nr+1)       &       g_objf_state_final(i,j,bi,bj,4*Nr+1)
88           END DO           END DO
89          END DO          END DO
# Line 101  c--------------------------------------- Line 103  c---------------------------------------
103          end do          end do
104        end do        end do
105        il = ilnblnk(xx_etan_file)        il = ilnblnk(xx_etan_file)
106        write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)')        write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)')
107       &     xx_etan_file(1:il),'.',optimcycle       &     xx_etan_file(1:il),'.',optimcycle
108        call adactive_read_xy_loc( fnamegeneric,1,doglobalread,ladinit,        call adactive_read_xy_loc( fnamegeneric,1,doglobalread,ladinit,
109       &     optimcycle,mythid,tmpfld2d )       &     optimcycle,mythid,tmpfld2d )
# Line 119  c-- Line 121  c--
121          end do          end do
122        end do        end do
123        il = ilnblnk(xx_vvel_file)        il = ilnblnk(xx_vvel_file)
124        write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)')        write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)')
125       &     xx_vvel_file(1:il),'.',optimcycle       &     xx_vvel_file(1:il),'.',optimcycle
126        call adactive_read_xyz( fnamegeneric,1,doglobalread,ladinit,        call adactive_read_xyz( fnamegeneric,1,doglobalread,ladinit,
127       &     optimcycle,mythid,tmpfld3d )       &     optimcycle,mythid,tmpfld3d )
# Line 137  c-- Line 139  c--
139          end do          end do
140        end do        end do
141        il = ilnblnk(xx_uvel_file)        il = ilnblnk(xx_uvel_file)
142        write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)')        write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)')
143       &     xx_uvel_file(1:il),'.',optimcycle       &     xx_uvel_file(1:il),'.',optimcycle
144        call adactive_read_xyz( fnamegeneric,1,doglobalread,ladinit,        call adactive_read_xyz( fnamegeneric,1,doglobalread,ladinit,
145       &     optimcycle,mythid,tmpfld3d )       &     optimcycle,mythid,tmpfld3d )
# Line 155  c-- Line 157  c--
157          end do          end do
158        end do        end do
159        il = ilnblnk(xx_salt_file)        il = ilnblnk(xx_salt_file)
160        write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)')        write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)')
161       &     xx_salt_file(1:il),'.',optimcycle       &     xx_salt_file(1:il),'.',optimcycle
162        call adactive_read_xyz_loc( fnamegeneric,1,doglobalread,ladinit,        call adactive_read_xyz_loc( fnamegeneric,1,doglobalread,ladinit,
163       &     optimcycle,mythid,tmpfld3d )       &     optimcycle,mythid,tmpfld3d )
# Line 173  c-- Line 175  c--
175          end do          end do
176        end do        end do
177        il = ilnblnk(xx_theta_file)        il = ilnblnk(xx_theta_file)
178        write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)')        write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)')
179       &     xx_theta_file(1:il),'.',optimcycle       &     xx_theta_file(1:il),'.',optimcycle
180        call adactive_read_xyz_loc( fnamegeneric,1,doglobalread,ladinit,        call adactive_read_xyz_loc( fnamegeneric,1,doglobalread,ladinit,
181       &     optimcycle,mythid,tmpfld3d )       &     optimcycle,mythid,tmpfld3d )
182    
183  #endif  #endif
184          
185        end        end
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22