Changeset 2258 for LMDZ5/branches/testing/libf/phylmd/phyetat0.F90
- Timestamp:
- Apr 13, 2015, 10:21:09 AM (9 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk merged: 2218,2221-2237
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/phylmd/phyetat0.F90
r2220 r2258 11 11 USE phys_state_var_mod, ONLY : ancien_ok, clwcon, detr_therm, dtime, & 12 12 du_gwd_rando, dv_gwd_rando, entr_therm, f0, falb1, falb2, fm_therm, & 13 falb_dir, falb_dif, & 13 14 ftsol, pbl_tke, pctsrf, q_ancien, radpas, radsol, rain_fall, ratqs, & 14 15 rlat, rlon, rnebcon, rugoro, sig1, snow_fall, solaire_etat0, sollw, sollwdown, & … … 69 70 INTEGER length 70 71 PARAMETER (length=100) 71 INTEGER it, iiq 72 INTEGER it, iiq, isw 72 73 REAL tab_cntrl(length), tabcntr0(length) 73 74 CHARACTER*7 str7 … … 243 244 ENDIF 244 245 246 !=================================================================== 247 ! Lecture des albedo difus et direct 248 249 DO nsrf = 1, nbsrf 250 DO isw=1, nsw 251 IF (isw.GT.99 .AND. nsrf.GT.99) THEN 252 PRINT*, "Trop de bandes SW ou sous-mailles" 253 call abort_gcm("phyetat0", "", 1) 254 ENDIF 255 WRITE(str7, '(i2.2, "srf", i2.2)') isw, nsrf 256 257 CALL get_field('A_dir_SW'//str7, falb_dir(:, isw, nsrf), found) 258 IF (.NOT. found) THEN 259 PRINT*, "phyetat0: Le champ <A_dir_SW"//str7//"> est absent" 260 PRINT*, " Il prend donc la valeur de surface" 261 DO i=1, klon 262 falb_dir(i, isw, nsrf)=0.2 263 ENDDO 264 ENDIF 265 CALL get_field('A_dif_SW'//str7, falb_dif(:, isw, nsrf), found) 266 IF (.NOT. found) THEN 267 PRINT*, "phyetat0: Le champ <A_dif_SW"//str7//"> est absent" 268 PRINT*, " Il prend donc la valeur de surface" 269 DO i=1, klon 270 falb_dif(i, isw, nsrf)=0.2 271 ENDDO 272 ENDIF 273 ENDDO 274 ENDDO 275 276 !=================================================================== 245 277 ! Lecture des temperatures du sol profond: 246 278 … … 264 296 ENDDO 265 297 298 !=================================================================== 266 299 ! Lecture de l'humidite de l'air juste au dessus du sol: 267 300
Note: See TracChangeset
for help on using the changeset viewer.