/[MITgcm]/MITgcm/model/src/ini_masks_etc.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_masks_etc.F

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

revision 1.21 by adcroft, Tue May 29 14:01:37 2001 UTC revision 1.22 by jmc, Thu Jul 5 21:44:25 2001 UTC
# Line 29  C     myThid -  Number of this instance Line 29  C     myThid -  Number of this instance
29        INTEGER myThid        INTEGER myThid
30  CEndOfInterface  CEndOfInterface
31    
32    C     == Local variables in common ==
33    C     tmpfld  - Temporary array used to compute & write Total Depth
34    C               has to be in common for multi threading
35          COMMON / LOCAL_INI_MASKS_ETC / tmpfld
36          _RS tmpfld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
37    
38  C     == Local variables ==  C     == Local variables ==
39  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
40  C     I,J,K  C     I,J,K
# Line 40  C     I,J,K Line 46  C     I,J,K
46  #endif  #endif
47        _RL hFacCtmp        _RL hFacCtmp
48        _RL hFacMnSz        _RL hFacMnSz
       _RS tmpfld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
49    
50  C- Calculate lopping factor hFacC : over-estimate the part inside of the domain  C- Calculate lopping factor hFacC : over-estimate the part inside of the domain
51  C    taking into account the lower_R Boundary (Bathymetrie / Top of Atmos)  C    taking into account the lower_R Boundary (Bathymetrie / Top of Atmos)
# Line 157  C         Inverse of fluid column thickn Line 162  C         Inverse of fluid column thickn
162         ENDDO         ENDDO
163        ENDDO        ENDDO
164  C     _EXCH_XY_R4(   recip_Rcol, myThid )  C     _EXCH_XY_R4(   recip_Rcol, myThid )
       CALL WRITE_FLD_XY_RS( 'Depth',' ',tmpfld,0,myThid)  
       CALL WRITE_FLD_XYZ_RS( 'hFacC',' ',hFacC,0,myThid)  
 C     CALL MDSWRITEFIELD( 'Depth', writeBinaryPrec, .TRUE.,  
 C    &                    'RS', 1, tmpfld, 1, -1, myThid )  
165    
166  C     hFacW and hFacS (at U and V points)  C     hFacW and hFacS (at U and V points)
167        DO bj=myByLo(myThid), myByHi(myThid)        DO bj=myByLo(myThid), myByHi(myThid)
# Line 194  C     We should really supply a flag for Line 195  C     We should really supply a flag for
195         ENDDO         ENDDO
196        ENDDO        ENDDO
197    
198    C-    Write to disk: Total Column Thickness & hFac(C,W,S):
199          _BARRIER
200          _BEGIN_MASTER( myThid )
201    C     CALL MDSWRITEFIELD( 'Depth', writeBinaryPrec, .TRUE.,
202    C    &                    'RS', 1, tmpfld, 1, -1, myThid )
203          CALL WRITE_FLD_XY_RS( 'Depth',' ',tmpfld,0,myThid)
204          CALL WRITE_FLD_XYZ_RS( 'hFacC',' ',hFacC,0,myThid)
205          CALL WRITE_FLD_XYZ_RS( 'hFacW',' ',hFacW,0,myThid)
206          CALL WRITE_FLD_XYZ_RS( 'hFacS',' ',hFacS,0,myThid)
207          _END_MASTER(myThid)
208    
209        CALL PLOT_FIELD_XYZRS( hFacC, 'hFacC' , Nr, 1, myThid )        CALL PLOT_FIELD_XYZRS( hFacC, 'hFacC' , Nr, 1, myThid )
210        CALL PLOT_FIELD_XYZRS( hFacW, 'hFacW' , Nr, 1, myThid )        CALL PLOT_FIELD_XYZRS( hFacW, 'hFacW' , Nr, 1, myThid )
211        CALL PLOT_FIELD_XYZRS( hFacS, 'hFacS' , Nr, 1, myThid )        CALL PLOT_FIELD_XYZRS( hFacS, 'hFacS' , Nr, 1, myThid )
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22