/[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.4 by cnh, Tue May 29 19:28:53 2001 UTC revision 1.5 by cnh, Mon Jun 18 17:39:58 2001 UTC
# Line 288  C      WRITE(fNam,'(A,A,A)' ) 'salb.',mn Line 288  C      WRITE(fNam,'(A,A,A)' ) 'salb.',mn
288  C      OPEN(1,FILE=fNam(1:14),STATUS='old',FORM='unformatted')  C      OPEN(1,FILE=fNam(1:14),STATUS='old',FORM='unformatted')
289  C      READ(1) tmp4  C      READ(1) tmp4
290  C      CLOSE(1)  C      CLOSE(1)
291         DO J=1,nYio  C      DO J=1,nYio
292          DO I=1,nXio  C       DO I=1,nXio
293           tmp4(I,J) = aim_albedo(I,J)/100.  C        tmp4(I,J) = aim_albedo(I,J)/100.
294          ENDDO  C       ENDDO
295         ENDDO  C      ENDDO
296         DO J=1,sNy         DO J=1,sNy
297          DO I=1,sNx          DO I=1,sNx
298           I2 = (sNx)*(J-1)+I           I2 = (sNx)*(J-1)+I
299           alb0(I2,myThid) = 0.           alb0(I2,myThid) = 0.
300           IF ( IG0+I-1 .LE. nxIo .AND. JG0+J-1 .LE. nyIo ) THEN           alb0(I2,myThid) = aim_albedo(I,J,bi,bj)/100.
           alb0(I2,myThid) = tmp4(IG0+I-1,JG0+J-1)  
          ENDIF  
301          ENDDO          ENDDO
302         ENDDO         ENDDO
303  C      Read in surface temperature data (input is in absolute temperature)  C      Read in surface temperature data (input is in absolute temperature)
# Line 312  C      CLOSE(1) Line 310  C      CLOSE(1)
310           I2 = (sNx)*(J-1)+I           I2 = (sNx)*(J-1)+I
311           sst1(I2,myThid) = 300.           sst1(I2,myThid) = 300.
312           stl1(I2,myThid) = 300.           stl1(I2,myThid) = 300.
313           IF ( IG0+I-1 .LE. nxIo .AND. JG0+J-1 .LE. nyIo ) THEN           sst1(I2,myThid) = aim_surfTemp(I,J,bi,bj)
314            sst1(I2,myThid) = aim_surfTemp(IG0+I-1,JG0+J-1)           stl1(I2,myThid) = aim_surfTemp(I,J,bi,bj)
           stl1(I2,myThid) = aim_surfTemp(IG0+I-1,JG0+J-1)  
          ENDIF  
315          ENDDO          ENDDO
316         ENDDO         ENDDO
317  C  C
# Line 331  cdj       tmp4 = (tmp4*7.5/20.)*10. Line 327  cdj       tmp4 = (tmp4*7.5/20.)*10.
327          DO I=1,sNx          DO I=1,sNx
328           I2 = (sNx)*(J-1)+I           I2 = (sNx)*(J-1)+I
329           soilq1(I2,myThid) = 0.           soilq1(I2,myThid) = 0.
330           IF ( IG0+I-1 .LE. nxIo .AND. JG0+J-1 .LE. nyIo ) THEN           soilq1(I2,myThid) = aim_soilMoisture(I,J,bi,bj)/20.
           soilq1(I2,myThid) = aim_soilMoisture(IG0+I-1,JG0+J-1)/20.  
          ENDIF  
331          ENDDO          ENDDO
332         ENDDO         ENDDO
333  C_cnh01      ENDIF  C_cnh01      ENDIF
# Line 368  Caja    ENDDO Line 362  Caja    ENDDO
362  Caja   ENDDO  Caja   ENDDO
363  Caja  ENDDO  Caja  ENDDO
364    
365    
366        CALL PDRIVER( tYear, myThid )        CALL PDRIVER( tYear, myThid )
367    
368  #ifdef ALLOW_TIMEAVE  #ifdef ALLOW_TIMEAVE
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22