/[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.17 by adcroft, Mon Mar 27 22:25:44 2000 UTC revision 1.18 by heimbach, Mon Sep 11 23:02:31 2000 UTC
# Line 34  C     I,J Line 34  C     I,J
34         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
35          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
36           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
37            fu   (i,j,bi,bj) = 0. _d 0            fu              (i,j,bi,bj) = 0. _d 0
38            fv   (i,j,bi,bj) = 0. _d 0            fv              (i,j,bi,bj) = 0. _d 0
39            Qnet (i,j,bi,bj) = 0. _d 0            Qnet            (i,j,bi,bj) = 0. _d 0
40            EmPmR(i,j,bi,bj) = 0. _d 0            EmPmR           (i,j,bi,bj) = 0. _d 0
41            SST  (i,j,bi,bj) = 0. _d 0            SST             (i,j,bi,bj) = 0. _d 0
42            SSS  (i,j,bi,bj) = 0. _d 0            SSS             (i,j,bi,bj) = 0. _d 0
43            Qsw  (i,j,bi,bj) = 0. _d 0            Qsw             (i,j,bi,bj) = 0. _d 0
44              surfaceTendencyU(i,j,bi,bj) = 0. _d 0
45              surfaceTendencyV(i,j,bi,bj) = 0. _d 0
46              surfaceTendencyT(i,j,bi,bj) = 0. _d 0
47              surfaceTendencyS(i,j,bi,bj) = 0. _d 0
48           ENDDO           ENDDO
49          ENDDO          ENDDO
50         ENDDO         ENDDO
# Line 59  C Line 63  C
63        IF ( EmPmRfile .NE. ' '  ) THEN        IF ( EmPmRfile .NE. ' '  ) THEN
64         CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid )         CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid )
65        ENDIF        ENDIF
66    #ifdef SHORTWAVE_HEATING
67          IF ( surfQswFile .NE. ' '  ) THEN
68           CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )
69          ENDIF
70    #endif
71        IF ( thetaClimFile .NE. ' '  ) THEN        IF ( thetaClimFile .NE. ' '  ) THEN
72         CALL READ_FLD_XY_RS( thetaClimFile, ' ', SST, 0, myThid )         CALL READ_FLD_XY_RS( thetaClimFile, ' ', SST, 0, myThid )
73        ENDIF        ENDIF
74        IF ( saltClimFile .NE. ' '  ) THEN        IF ( saltClimFile .NE. ' '  ) THEN
75         CALL READ_FLD_XY_RS( saltClimFile, ' ', SSS, 0, myThid )         CALL READ_FLD_XY_RS( saltClimFile, ' ', SSS, 0, myThid )
76        ENDIF        ENDIF
77  c     IF ( surfQswFile .NE. ' '  ) THEN  #ifdef SHORTWAVE_HEATING
78  c      CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )        IF ( surfQswFile .NE. ' '  ) THEN
79  c     ENDIF         CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )
80          ENDIF
81    #endif
82        _END_MASTER(myThid)        _END_MASTER(myThid)
83  C  C
84        _EXCH_XY_R4(fu   , myThid )        _EXCH_XY_R4(fu   , myThid )
# Line 81  C Line 92  C
92  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)
93  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)
94    
95    C     Convert units and signs (!)
96          CALL EXTERNAL_FIELDS_SCALE( myThid )
97    
98        RETURN        RETURN
99        END        END
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22