/[MITgcm]/MITgcm/pkg/ebm/ebm_load_climatology.F
ViewVC logotype

Diff of /MITgcm/pkg/ebm/ebm_load_climatology.F

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

revision 1.5 by jmc, Sun Aug 28 17:35:38 2011 UTC revision 1.6 by jmc, Sun Aug 28 21:47:33 2011 UTC
# Line 36  C     !LOCAL VARIABLES: Line 36  C     !LOCAL VARIABLES:
36        INTEGER bi,bj,i,j        INTEGER bi,bj,i,j
37        _RL distY,tauX,tauMax,lY        _RL distY,tauX,tauMax,lY
38    
       IF ( saltClimFile .NE. ' ' ) THEN  
        CALL READ_FLD_XY_RS( saltClimFile, ' ', SSS, 0, myThid )  
       ENDIF  
   
       IF ( thetaClimFile .NE. ' ' ) THEN  
        CALL READ_FLD_XY_RS( thetaClimFile, ' ', SST, 0, myThid )  
       ENDIF  
   
39        IF ( RunoffFile .NE. ' ' ) THEN        IF ( RunoffFile .NE. ' ' ) THEN
40         CALL READ_FLD_XY_RS( RunoffFile, ' ', Run, 0, myThid )         CALL READ_FLD_XY_RS( RunoffFile, ' ', Run, 0, myThid )
41        ENDIF        ENDIF
# Line 72  c          tauX  = -tauMax*cos(2. _d 0*P Line 64  c          tauX  = -tauMax*cos(2. _d 0*P
64           ENDDO           ENDDO
65          ENDDO          ENDDO
66         ENDDO         ENDDO
       ELSE  
        CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid )  
       ENDIF  
   
       IF (meridWindFile .EQ. ' ' ) THEN  
        DO bj = myByLo(myThid), myByHi(myThid)  
         DO bi = myBxLo(myThid), myBxHi(myThid)  
          DO j=1,sNy  
           DO i=1,sNx  
            fv(i,j,bi,bj) = 0.0  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
       ELSE  
        CALL READ_FLD_XY_RS( meridWindFile, ' ', fv, 0, myThid )  
67        ENDIF        ENDIF
68    
       _EXCH_XY_RS( SSS, myThid )  
       _EXCH_XY_RS( SST, myThid )  
69        _EXCH_XY_RS( Run, myThid )        _EXCH_XY_RS( Run, myThid )
       CALL EXCH_UV_XY_RS( fu, fv, .TRUE., myThid )  
70    
71        IF ( debugLevel.GE.debLevC ) THEN        IF ( debugLevel.GE.debLevC ) THEN
72          CALL PLOT_FIELD_XYRS( SST, 'Theta Climatology' , 1, myThid )          CALL PLOT_FIELD_XYRS( SST, 'Theta Climatology' , 1, myThid )
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22