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

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

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

revision 1.48 by dimitri, Wed Nov 28 00:18:17 2007 UTC revision 1.49 by jmc, Fri May 16 01:53:27 2008 UTC
# Line 23  C     === Global variables === Line 23  C     === Global variables ===
23  #include "EEPARAMS.h"  #include "EEPARAMS.h"
24  #include "PARAMS.h"  #include "PARAMS.h"
25  #include "GRID.h"  #include "GRID.h"
26    #include "SURFACE.h"
27  #include "FFIELDS.h"  #include "FFIELDS.h"
28    
29  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
# Line 53  C-    Initialise all arrays in common bl Line 54  C-    Initialise all arrays in common bl
54            Qsw             (i,j,bi,bj) = 0. _d 0            Qsw             (i,j,bi,bj) = 0. _d 0
55            pLoad           (i,j,bi,bj) = 0. _d 0            pLoad           (i,j,bi,bj) = 0. _d 0
56            sIceLoad        (i,j,bi,bj) = 0. _d 0            sIceLoad        (i,j,bi,bj) = 0. _d 0
57            surfaceForcingU(i,j,bi,bj) = 0. _d 0            surfaceForcingU (i,j,bi,bj) = 0. _d 0
58            surfaceForcingV(i,j,bi,bj) = 0. _d 0            surfaceForcingV (i,j,bi,bj) = 0. _d 0
59            surfaceForcingT(i,j,bi,bj) = 0. _d 0            surfaceForcingT (i,j,bi,bj) = 0. _d 0
60            surfaceForcingS(i,j,bi,bj) = 0. _d 0            surfaceForcingS (i,j,bi,bj) = 0. _d 0
61            surfaceForcingTice(i,j,bi,bj) = 0. _d 0            surfaceForcingTice(i,j,bi,bj) = 0. _d 0
62  #ifndef EXCLUDE_FFIELDS_LOAD  #ifndef EXCLUDE_FFIELDS_LOAD
63            taux0           (i,j,bi,bj) = 0. _d 0            taux0           (i,j,bi,bj) = 0. _d 0
# Line 197  C     CALL PLOT_FIELD_XYRS( pLoad, 'S/R Line 198  C     CALL PLOT_FIELD_XYRS( pLoad, 'S/R
198  C     CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)  C     CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)
199  C     CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)  C     CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)
200    
201    #ifdef ATMOSPHERIC_LOADING
202          IF ( pLoadFile .NE. ' ' .AND. usingPCoords ) THEN
203    C-- This is a hack used to read phi0surf from a file (pLoadFile)
204    C          instead of computing it from bathymetry & density ref. profile.
205    C-  Ocean: The true atmospheric P-loading is not yet implemented for P-coord
206    C          (requires time varying dP(Nr) like dP(k-bottom) with NonLin FS).
207    C-  Atmos: sometime usefull to overwrite phi0surf with fixed-in-time field
208    C          read from file (and anyway, pressure loading is meaningless here)
209            DO bj = myByLo(myThid), myByHi(myThid)
210             DO bi = myBxLo(myThid), myBxHi(myThid)
211              DO j=1-OLy,sNy+OLy
212               DO i=1-OLx,sNx+OLx
213                 phi0surf(i,j,bi,bj) = pLoad(i,j,bi,bj)
214               ENDDO
215              ENDDO
216             ENDDO
217            ENDDO
218          ENDIF
219    #endif /* ATMOSPHERIC_LOADING */
220    
221        RETURN        RETURN
222        END        END
Legend:
Removed from v.1.48  
changed lines
  Added in v.1.49
  ViewVC Help
Powered by ViewVC 1.1.22