Changeset 1151 for LMDZ4/branches/LMDZ4-dev/libf/phylmd
- Timestamp:
- Apr 28, 2009, 6:33:49 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol.F90 ¶
r1150 r1151 143 143 WRITE(lunout,*) 'get_aero iyr=', iyr,' ',cyear 144 144 CALL get_aero_fromfile(cyear, var_1, name_aero(id_aero)) 145 ELSE 146 ! Read in data:145 146 ELSE ! Read data from 2 decades 147 147 ! a) from actual 10-yr-period 148 148 IF (iyr.LT.1900) THEN … … 159 159 WRITE(cyear,'(I4)') iyr1 160 160 WRITE(lunout,*) 'get_aero iyr=', iyr,' ',cyear 161 162 161 CALL get_aero_fromfile(cyear, var_1, name_aero(id_aero)) 163 162 164 ENDIF165 ! If to read two decades:166 IF (.NOT.lonlyone) THEN167 168 ! b) from the next following one169 WRITE(cyear,'(I4)') iyr2170 WRITE(lunout,*) 'get_aero iyr=', iyr,' ',cyear171 172 CALL get_aero_fromfile(cyear, var_2, name_aero(id_aero))173 174 ! Interpolate linarily to the actual year:175 DO it=1,12176 DO k=1,klev177 DO j=1,jjm178 DO i=1,iim179 var_1(i,j,k,it) =&180 var_1(i,j,k,it) - FLOAT(iyr-iyr1)/FLOAT(iyr2-iyr1) * &181 (var_1(i,j,k,it) - var_2(i,j,k,it))163 ! If to read two decades: 164 IF (.NOT.lonlyone) THEN 165 166 ! b) from the next following one 167 WRITE(cyear,'(I4)') iyr2 168 WRITE(lunout,*) 'get_aero iyr=', iyr,' ',cyear 169 170 CALL get_aero_fromfile(cyear, var_2, name_aero(id_aero)) 171 172 ! Interpolate linarily to the actual year: 173 DO it=1,12 174 DO k=1,klev 175 DO j=1,jjm 176 DO i=1,iim 177 var_1(i,j,k,it) = & 178 var_1(i,j,k,it) - FLOAT(iyr-iyr1)/FLOAT(iyr2-iyr1) * & 179 (var_1(i,j,k,it) - var_2(i,j,k,it)) 180 ENDDO 182 181 ENDDO 183 182 ENDDO 184 183 ENDDO 185 ENDDO186 184 187 ENDIF !lonlyone 185 ENDIF ! lonlyone 186 ENDIF ! iyr .LT. 1850 188 187 ENDIF ! aer_type 189 188
Note: See TracChangeset
for help on using the changeset viewer.