/[MITgcm]/MITgcm_contrib/jscott/pkg_atm2d/init_atm2d.F
ViewVC logotype

Diff of /MITgcm_contrib/jscott/pkg_atm2d/init_atm2d.F

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

revision 1.1 by jscott, Fri Aug 11 18:55:50 2006 UTC revision 1.2 by jscott, Tue Aug 22 20:21:39 2006 UTC
# Line 18  c Line 18  c
18    
19  C     === Global Atmosphere Variables ===  C     === Global Atmosphere Variables ===
20  #include "ATMSIZE.h"  #include "ATMSIZE.h"
21  #include "AGRID.COM"  #include "AGRID.h"
22    
23  C     === Global Ocean Variables ===  C     === Global Ocean Variables ===
24  #include "SIZE.h"  #include "SIZE.h"
# Line 44  C     LOCAL VARIABLES: Line 44  C     LOCAL VARIABLES:
44        INTEGER i,j,jj        INTEGER i,j,jj
45        INTEGER ib, ibj1, ibj2        INTEGER ib, ibj1, ibj2
46        INTEGER j_atm, mn        INTEGER j_atm, mn
47          INTEGER dUnit
48        _RL end1, end2        _RL end1, end2
49        _RL totrun_b(1:sNy)        _RL totrun_b(sNy)
50        _RL a1,a2        _RL a1,a2
51        _RS atm_dyG(1:jm0)        _RS atm_dyG(jm0)
52        DATA atm_dyG/2.0,44*4.0,2.0/        DATA atm_dyG/2.0,44*4.0,2.0/
53    
54        dtatmo = dtatm * 3600.        dtatmo = dtatm * 3600.
# Line 156  c Line 157  c
157          ENDDO          ENDDO
158        ENDDO        ENDDO
159    
160          CALL MDSFINDUNIT( dUnit, myThid )
161    
162        IF ( atmosTauuFile .NE. ' '  ) THEN        IF ( atmosTauuFile .NE. ' '  ) THEN
163           OPEN(UNIT=97, FILE=atmosTauuFile,STATUS='old',           OPEN(dUnit, FILE=atmosTauuFile,STATUS='old',
164       &        ACCESS='direct', RECL=8*jm0*nForcingPer,       &        ACCESS='direct', RECL=8*jm0*nForcingPer,
165       &        FORM='unformatted')       &        FORM='unformatted')
166           READ(97,REC=1), atau           READ(dUnit,REC=1), atau
167           CLOSE(97)           CLOSE(dUnit)
168        ENDIF        ENDIF
169    
170        IF ( atmosTauvFile .NE. ' '  ) THEN        IF ( atmosTauvFile .NE. ' '  ) THEN
171           OPEN(UNIT=98, FILE=atmosTauvFile, STATUS='old',           OPEN(dUnit, FILE=atmosTauvFile, STATUS='old',
172       &        ACCESS='direct', RECL=8*jm0*nForcingPer,       &        ACCESS='direct', RECL=8*jm0*nForcingPer,
173       &        FORM='unformatted')       &        FORM='unformatted')
174           READ(98, REC=1), atav           READ(dUnit, REC=1), atav
175           CLOSE(98)           CLOSE(dUnit)
176        ENDIF        ENDIF
177    
178        IF ( atmosWindFile .NE. ' '  ) THEN        IF ( atmosWindFile .NE. ' '  ) THEN
179           OPEN(UNIT=99, FILE=atmosWindFile, STATUS='old',           OPEN(dUnit, FILE=atmosWindFile, STATUS='old',
180       &        ACCESS='direct', RECL=8*jm0*nForcingPer,       &        ACCESS='direct', RECL=8*jm0*nForcingPer,
181       &        FORM='unformatted')       &        FORM='unformatted')
182           READ(99, REC=1), awind           READ(dUnit, REC=1), awind
183           CLOSE(99)           CLOSE(dUnit)
184        ENDIF        ENDIF
185    
186  C The polar data point values for winds are effectively N/A given the  C The polar data point values for winds are effectively N/A given the
# Line 251  C     Initialize 2D diagnostic variables Line 254  C     Initialize 2D diagnostic variables
254              runoff_atm_ta(i,j,mn)= 0. _d 0              runoff_atm_ta(i,j,mn)= 0. _d 0
255              sum_qrel_ta(i,j,mn)= 0. _d 0              sum_qrel_ta(i,j,mn)= 0. _d 0
256              sum_frel_ta(i,j,mn)= 0. _d 0              sum_frel_ta(i,j,mn)= 0. _d 0
257                sum_iceMask_ta(i,j,mn)= 0. _d 0
258                sum_iceHeight_ta(i,j,mn)= 0. _d 0
259                sum_iceTime_ta(i,j,mn)= 0. _d 0
260                sum_oceMxLT_ta(i,j,mn)= 0. _d 0
261                sum_oceMxLS_ta(i,j,mn)= 0. _d 0
262            ENDDO            ENDDO
263            qnet_atm(i,j)= 0. _d 0            qnet_atm(i,j)= 0. _d 0
264            evap_atm(i,j)= 0. _d 0            evap_atm(i,j)= 0. _d 0
# Line 258  C     Initialize 2D diagnostic variables Line 266  C     Initialize 2D diagnostic variables
266            runoff_atm(i,j)= 0. _d 0            runoff_atm(i,j)= 0. _d 0
267            sum_qrel(i,j)= 0. _d 0            sum_qrel(i,j)= 0. _d 0
268            sum_frel(i,j)= 0. _d 0            sum_frel(i,j)= 0. _d 0
269              sum_iceMask(i,j)= 0. _d 0
270              sum_iceHeight(i,j)= 0. _d 0
271              sum_iceTime(i,j)= 0. _d 0
272              sum_oceMxLT(i,j)= 0. _d 0
273              sum_oceMxLS(i,j)= 0. _d 0
274          ENDDO          ENDDO
275        ENDDO        ENDDO
276    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22