/[MITgcm]/MITgcm/pkg/aim/aim_do_atmos_physics.F
ViewVC logotype

Diff of /MITgcm/pkg/aim/aim_do_atmos_physics.F

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

revision 1.8 by jmc, Tue Sep 25 19:55:34 2001 UTC revision 1.8.2.1 by adcroft, Tue Feb 26 16:04:48 2002 UTC
# Line 173  c_jmc: Because AIM physics LSC is not ap Line 173  c_jmc: Because AIM physics LSC is not ap
173  c      ==> move water wapor from the stratos to the surface level.  c      ==> move water wapor from the stratos to the surface level.
174        DO J = 1-Oly, sNy+Oly        DO J = 1-Oly, sNy+Oly
175         DO I = 1-Olx, sNx+Olx         DO I = 1-Olx, sNx+Olx
176  c       k = k_surf(i,j,bi,bj)          k = ksurfC(i,j,bi,bj)
177  c       salt(I,J,k,bi,bj) = salt(I,J,k,bi,bj)          IF (k.LE.Nr)
178  c    &   + maskC(i,j,Nr,bi,bj)*salt(I,J,Nr,bi,bj)*drF(Nr)*recip_drF(k)       &    salt(I,J,k,bi,bj) = salt(I,J,k,bi,bj)
179         &                      + salt(I,J,Nr,bi,bj)*drF(Nr)*recip_drF(k)
180          salt(I,J,Nr,bi,bj) = 0.          salt(I,J,Nr,bi,bj) = 0.
181         ENDDO         ENDDO
182        ENDDO        ENDDO
# Line 282  C     Physics package needs to know time Line 283  C     Physics package needs to know time
283    
284  C  C
285  C     Load external data needed by physics package  C     Load external data needed by physics package
286  C     1. Albedo  C     1. Albedo                (between 0-1)
287  C     2. Soil moisture  C     2. Soil moisture         (between 0-1)
288  C     3. Surface temperatures  C     3. Surface temperatures  (in situ Temp. [K])
289  C     4. Snow depth            - assume no snow for now  C     4. Snow depth            - assume no snow for now
290  C     5. Sea ice               - assume no sea ice for now  C     5. Sea ice               - assume no sea ice for now
291  C     6. Land sea mask         - infer from exact zeros in soil moisture dataset  C     6. Land sea mask         - infer from exact zeros in soil moisture dataset
# Line 309  C      ENDDO Line 310  C      ENDDO
310          DO I=1,sNx          DO I=1,sNx
311           I2 = (sNx)*(J-1)+I           I2 = (sNx)*(J-1)+I
312           alb0(I2,myThid) = 0.           alb0(I2,myThid) = 0.
313           alb0(I2,myThid) = aim_albedo(I,J,bi,bj)/100.  c        alb0(I2,myThid) = aim_albedo(I,J,bi,bj)/100.
314             alb0(I2,myThid) = aim_albedo(I,J,bi,bj)
315          ENDDO          ENDDO
316         ENDDO         ENDDO
317  C      Read in surface temperature data (input is in absolute temperature)  C      Read in surface temperature data (input is in absolute temperature)
# Line 339  cdj       tmp4 = (tmp4*7.5/20.)*10. Line 341  cdj       tmp4 = (tmp4*7.5/20.)*10.
341          DO I=1,sNx          DO I=1,sNx
342           I2 = (sNx)*(J-1)+I           I2 = (sNx)*(J-1)+I
343           soilq1(I2,myThid) = 0.           soilq1(I2,myThid) = 0.
344           soilq1(I2,myThid) = aim_soilMoisture(I,J,bi,bj)/20.  c        soilq1(I2,myThid) = aim_soilMoisture(I,J,bi,bj)/20.
345             soilq1(I2,myThid) = aim_soilMoisture(I,J,bi,bj)
346          ENDDO          ENDDO
347         ENDDO         ENDDO
348  C_cnh01      ENDIF  C_cnh01      ENDIF
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.8.2.1
  ViewVC Help
Powered by ViewVC 1.1.22