/[MITgcm]/MITgcm/verification/natl_box/code/external_fields_load.F
ViewVC logotype

Diff of /MITgcm/verification/natl_box/code/external_fields_load.F

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

revision 1.8 by edhill, Sat Nov 1 04:50:03 2003 UTC revision 1.9 by jmc, Fri Jul 16 01:23:19 2004 UTC
# Line 133  C      data for the period ahead and the Line 133  C      data for the period ahead and the
133       &        'RS', 1, Qnet0, intime0, myThid )       &        'RS', 1, Qnet0, intime0, myThid )
134         CALL MDSREADFIELD ( surfQFile, readBinaryPrec,         CALL MDSREADFIELD ( surfQFile, readBinaryPrec,
135       &        'RS', 1, Qnet1, intime1, myThid )       &        'RS', 1, Qnet1, intime1, myThid )
136          ELSEIF ( surfQnetFile .NE. ' '  ) THEN
137           CALL MDSREADFIELD ( surfQnetFile, readBinaryPrec,
138         &        'RS', 1, Qnet0, intime0, myThid )
139           CALL MDSREADFIELD ( surfQnetFile, readBinaryPrec,
140         &        'RS', 1, Qnet1, intime1, myThid )
141        ENDIF        ENDIF
142        IF ( EmPmRfile .NE. ' '  ) THEN        IF ( EmPmRfile .NE. ' '  ) THEN
143         CALL MDSREADFIELD ( EmPmRfile, readBinaryPrec,         CALL MDSREADFIELD ( EmPmRfile, readBinaryPrec,
# Line 158  C      data for the period ahead and the Line 163  C      data for the period ahead and the
163       &        'RS', 1, Qsw0, intime0, myThid )       &        'RS', 1, Qsw0, intime0, myThid )
164         CALL MDSREADFIELD ( surfQswFile, readBinaryPrec,         CALL MDSREADFIELD ( surfQswFile, readBinaryPrec,
165       &        'RS', 1, Qsw1, intime1, myThid )       &        'RS', 1, Qsw1, intime1, myThid )
166           IF ( surfQFile .NE. ' '  ) THEN
167    C-     Qnet is now (after c54) the net Heat Flux (including SW)
168            DO bj=1,nSy
169             DO bi=1,nSx
170              DO j=1-Oly,sNy+Oly
171               DO i=1-Olx,sNx+Olx
172                Qnet0(i,j,bi,bj) = Qnet0(i,j,bi,bj) + Qsw0(i,j,bi,bj)
173                Qnet1(i,j,bi,bj) = Qnet1(i,j,bi,bj) + Qsw1(i,j,bi,bj)
174               ENDDO
175              ENDDO
176             ENDDO
177            ENDDO
178           ENDIF
179        ENDIF        ENDIF
180  #endif  #endif
181    
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22