/[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.15 by adcroft, Tue Dec 15 00:20:34 1998 UTC revision 1.16 by adcroft, Wed May 5 18:32:34 1999 UTC
# Line 42  C     I,J,K Line 42  C     I,J,K
42            EmPmR(i,j,bi,bj) = 0. _d 0            EmPmR(i,j,bi,bj) = 0. _d 0
43            SST  (i,j,bi,bj) = 0. _d 0            SST  (i,j,bi,bj) = 0. _d 0
44            SSS  (i,j,bi,bj) = 0. _d 0            SSS  (i,j,bi,bj) = 0. _d 0
45              Qsw  (i,j,bi,bj) = 0. _d 0
46           ENDDO           ENDDO
47          ENDDO          ENDDO
48         ENDDO         ENDDO
49        ENDDO        ENDDO
50  C  C
51          _BEGIN_MASTER(myThid)
52          IF ( zonalWindFile .NE. ' '  ) THEN
53           CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid )
54          ENDIF
55          IF ( meridWindFile .NE. ' '  ) THEN
56           CALL READ_FLD_XY_RS( meridWindFile, ' ', fv, 0, myThid )
57          ENDIF
58          IF ( surfQFile .NE. ' '  ) THEN
59           CALL READ_FLD_XY_RS( surfQFile, ' ', Qnet, 0, myThid )
60          ENDIF
61          IF ( EmPmRfile .NE. ' '  ) THEN
62           CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid )
63          ENDIF
64          IF ( thetaClimFile .NE. ' '  ) THEN
65           CALL READ_FLD_XY_RS( thetaClimFile, ' ', SST, 0, myThid )
66          ENDIF
67          IF ( saltClimFile .NE. ' '  ) THEN
68           CALL READ_FLD_XY_RS( saltClimFile, ' ', SSS, 0, myThid )
69          ENDIF
70    c     IF ( surfQswFile .NE. ' '  ) THEN
71    c      CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )
72    c     ENDIF
73          _END_MASTER(myThid)
74    C
75        _EXCH_XY_R4(fu   , myThid )        _EXCH_XY_R4(fu   , myThid )
76        _EXCH_XY_R4(fv   , myThid )        _EXCH_XY_R4(fv   , myThid )
77        _EXCH_XY_R4(Qnet , myThid )        _EXCH_XY_R4(Qnet , myThid )
78        _EXCH_XY_R4(EmPmR, myThid )        _EXCH_XY_R4(EmPmR, myThid )
79        _EXCH_XY_R4(SST  , myThid )        _EXCH_XY_R4(SST  , myThid )
80        _EXCH_XY_R4(SSS  , myThid )        _EXCH_XY_R4(SSS  , myThid )
81          _EXCH_XY_R4(Qsw  , myThid )
82    
83  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)
84  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)
85    
       DO bj = myByLo(myThid), myByHi(myThid)  
        DO bi = myBxLo(myThid), myBxHi(myThid)  
         DO j=1-OLy,sNy+OLy  
          DO i=1-OLx,sNx+OLx  
           fudot   (i,j,bi,bj) = 0. _d 0  
           fvdot   (i,j,bi,bj) = 0. _d 0  
           Qnetdot (i,j,bi,bj) = 0. _d 0  
           EmPmRdot(i,j,bi,bj) = 0. _d 0  
           SSTdot  (i,j,bi,bj) = 0. _d 0  
           SSSdot  (i,j,bi,bj) = 0. _d 0  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
 C  
       _EXCH_XY_R4(fudot   , myThid )  
       _EXCH_XY_R4(fvdot   , myThid )  
       _EXCH_XY_R4(Qnetdot , myThid )  
       _EXCH_XY_R4(EmPmRdot, myThid )  
       _EXCH_XY_R4(SSTdot  , myThid )  
       _EXCH_XY_R4(SSSdot  , myThid )  
   
86        RETURN        RETURN
87        END        END
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22