/[MITgcm]/MITgcm_contrib/gael/pkg/smooth2/smooth_init3D.F
ViewVC logotype

Diff of /MITgcm_contrib/gael/pkg/smooth2/smooth_init3D.F

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

revision 1.2 by jmc, Fri Nov 13 19:52:32 2009 UTC revision 1.3 by gforget, Mon Feb 15 23:14:47 2010 UTC
# Line 58  c vertical smoothing: Line 58  c vertical smoothing:
58        if (smooth3DsizeZ(smoothOpNb).EQ.3) then        if (smooth3DsizeZ(smoothOpNb).EQ.3) then
59        write(fnamegeneric(1:80),'(1a,i3.3)')        write(fnamegeneric(1:80),'(1a,i3.3)')
60       &    'smooth3DscalesZ',smoothOpNb       &    'smooth3DscalesZ',smoothOpNb
61        call mdsreadfield(fnamegeneric,smoothprec,'RL',nR,        CALL READ_REC_3D_RL(fnamegeneric,smoothprec,
62       & smooth3D_Lz,1,mythid)       &           Nr, smooth3D_Lz,1,1,mythid)
63        _EXCH_XYZ_RL( smooth3D_Lz, mythid )        _EXCH_XYZ_RL( smooth3D_Lz, mythid )
64        else        else
65        DO bj=jtlo,jthi        DO bj=jtlo,jthi
# Line 116  c hoizontal smoothing: Line 116  c hoizontal smoothing:
116        if (smooth3DsizeH(smoothOpNb).EQ.3) then        if (smooth3DsizeH(smoothOpNb).EQ.3) then
117        write(fnamegeneric(1:80),'(1a,i3.3)')        write(fnamegeneric(1:80),'(1a,i3.3)')
118       &    'smooth3DscalesH',smoothOpNb       &    'smooth3DscalesH',smoothOpNb
119        call mdsreadfield(fnamegeneric,smoothprec,'RL',nR,        CALL READ_REC_3D_RL(fnamegeneric,smoothprec,
120       & smooth3D_Lx,1,mythid)       &           Nr, smooth3D_Lx,1,1,mythid)
121        call mdsreadfield(fnamegeneric,smoothprec,'RL',nR,        CALL READ_REC_3D_RL(fnamegeneric,smoothprec,
122       & smooth3D_Ly,2,mythid)       &           Nr, smooth3D_Ly,2,1,mythid)
123        _EXCH_XYZ_RL( smooth3D_Lx, mythid )        _EXCH_XYZ_RL( smooth3D_Lx, mythid )
124        _EXCH_XYZ_RL( smooth3D_Ly, mythid )        _EXCH_XYZ_RL( smooth3D_Ly, mythid )
125        else        else
# Line 178  c write diffusion operator to file Line 178  c write diffusion operator to file
178        write(fnamegeneric(1:80),'(1a,i3.3)')        write(fnamegeneric(1:80),'(1a,i3.3)')
179       &    'smooth3Doperator',smoothOpNb       &    'smooth3Doperator',smoothOpNb
180    
181        call mdswritefield(fnamegeneric,smoothprec,.false.,'RL',        CALL WRITE_REC_3D_RL(fnamegeneric,smoothprec,
182       & nR,smooth3D_Kwx,1,1,mythid)       &           Nr,smooth3D_Kwx,1,1,mythid)
183        call mdswritefield(fnamegeneric,smoothprec,.false.,'RL',        CALL WRITE_REC_3D_RL(fnamegeneric,smoothprec,
184       & nR,smooth3D_Kwy,2,1,mythid)       &           Nr,smooth3D_Kwy,2,1,mythid)
185        call mdswritefield(fnamegeneric,smoothprec,.false.,'RL',        CALL WRITE_REC_3D_RL(fnamegeneric,smoothprec,
186       & nR,smooth3D_Kwz,3,1,mythid)       &           Nr,smooth3D_Kwz,3,1,mythid)
187        call mdswritefield(fnamegeneric,smoothprec,.false.,'RL',        CALL WRITE_REC_3D_RL(fnamegeneric,smoothprec,
188       & nR,smooth3D_Kux,4,1,mythid)       &           Nr,smooth3D_Kux,4,1,mythid)
189        call mdswritefield(fnamegeneric,smoothprec,.false.,'RL',        CALL WRITE_REC_3D_RL(fnamegeneric,smoothprec,
190       & nR,smooth3D_Kvy,5,1,mythid)       &           Nr,smooth3D_Kvy,5,1,mythid)
191        call mdswritefield(fnamegeneric,smoothprec,.false.,'RL',        CALL WRITE_REC_3D_RL(fnamegeneric,smoothprec,
192       & nR,smooth3D_Kuz,6,1,mythid)       &           Nr,smooth3D_Kuz,6,1,mythid)
193        call mdswritefield(fnamegeneric,smoothprec,.false.,'RL',        CALL WRITE_REC_3D_RL(fnamegeneric,smoothprec,
194       & nR,smooth3D_Kvz,7,1,mythid)       &           Nr,smooth3D_Kvz,7,1,mythid)
195        call mdswritefield(fnamegeneric,smoothprec,.false.,'RL',        CALL WRITE_REC_3D_RL(fnamegeneric,smoothprec,
196       & nR,smooth3D_Kuy,8,1,mythid)       &           Nr,smooth3D_Kuy,8,1,mythid)
197        call mdswritefield(fnamegeneric,smoothprec,.false.,'RL',        CALL WRITE_REC_3D_RL(fnamegeneric,smoothprec,
198       & nR,smooth3D_Kvx,9,1,mythid)       &           Nr,smooth3D_Kvx,9,1,mythid)
199        call mdswritefield(fnamegeneric,smoothprec,.false.,'RL',        CALL WRITE_REC_3D_RL(fnamegeneric,smoothprec,
200       & nR,smooth3D_kappaR,10,1,mythid)       &           Nr,smooth3D_kappaR,10,1,mythid)
201    
202    
203        END        END
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22