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

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

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

revision 1.1 by cnh, Tue May 29 19:28:53 2001 UTC revision 1.2 by cnh, Mon Jun 18 17:39:58 2001 UTC
# Line 55  C     mnthNam Line 55  C     mnthNam
55        INTEGER prevMnthIndex        INTEGER prevMnthIndex
56        DATA    prevMnthIndex / 0 /        DATA    prevMnthIndex / 0 /
57        SAVE    prevMnthIndex        SAVE    prevMnthIndex
58        CHARACTER*16 fNam        CHARACTER*17 fNam
59        CHARACTER*3 mnthNam(12)        CHARACTER*3 mnthNam(12)
60        DATA mnthNam /        DATA mnthNam /
61       & 'jan', 'feb', 'mar', 'apr', 'may', 'jun',       & 'jan', 'feb', 'mar', 'apr', 'may', 'jun',
# Line 80  C       Calculate offset into a year Line 80  C       Calculate offset into a year
80  C        New month so load in data  C        New month so load in data
81           prevMnthIndex = mnthIndex           prevMnthIndex = mnthIndex
82  C        o Albedo ( convert % to fraction )  C        o Albedo ( convert % to fraction )
83           WRITE(fNam,'(A,A,A)' ) 'salb.',mnthNam(mnthIndex),'.sun.b'           WRITE(fNam,'(A,A,A)' ) 'salb.',mnthNam(mnthIndex),'.cs.bin'
84           OPEN(1,FILE=fNam(1:14),STATUS='old',FORM='unformatted')             CALL MDSREADFIELD(fNam(1:15),readBinaryPrec,'RS',1,
85           READ(1) aim_albedo       O    aim_albedo,  
86           CLOSE(1)       I    1,myThid)
          DO J=1,aim_nyIo  
           DO I=1,aim_nxIo  
 C          aim_albedo(I,J) = aim_albedo(I,J)/100.  
           ENDDO  
          ENDDO  
87    
88  C        o Surface temperature ( in kelvin )  C        o Surface temperature ( in kelvin )
89           WRITE(fNam,'(A,A,A)' ) 'tsurf.',mnthNam(mnthIndex),'.sun.b'           WRITE(fNam,'(A,A,A)' ) 'tsurf.',mnthNam(mnthIndex),'.cs.bin'
90           OPEN(1,FILE=fNam(1:15),STATUS='old',FORM='unformatted')           CALL MDSREADFIELD(fNam(1:16),readBinaryPrec,'RS',1,
91           READ(1) aim_surfTemp       O    aim_surftemp,  
92           CLOSE(1)       I    1,myThid)
93    
94  C        o Soil moisture ( convert to 20cm bucket fraction )  C        o Soil moisture ( convert to 20cm bucket fraction )
95           WRITE(fNam,'(A,A,A)' ) 'smoist.',mnthNam(mnthIndex),'.sun.b'           WRITE(fNam,'(A,A,A)' ) 'smoist.',mnthNam(mnthIndex),'.cs.bin'
96           OPEN(1,FILE=fNam(1:16),STATUS='old',FORM='unformatted')           CALL MDSREADFIELD(fNam(1:17),readBinaryPrec,'RS',1,
97           READ(1) aim_soilMoisture       O    aim_soilMoisture,  
98           CLOSE(1)       I    1,myThid)
          DO J=1,aim_nyIo  
           DO I=1,aim_nxIo  
 C          aim_soilMoisture(I,J) = aim_soilMoisture(I,J)/20.  
           ENDDO  
          ENDDO  
99    
100          ENDIF          ENDIF
101    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22