Changeset 2243 for LMDZ5/trunk/libf/phylmd/etat0_netcdf.F90
- Timestamp:
- Mar 24, 2015, 2:28:51 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/etat0_netcdf.F90
r2239 r2243 54 54 #include "dimsoil.h" 55 55 #include "temps.h" 56 REAL, DIMENSION(klon) :: tsol , qsol56 REAL, DIMENSION(klon) :: tsol 57 57 REAL, DIMENSION(klon) :: sn, rugmer, run_off_lic_0 58 58 REAL, DIMENSION(iip1,jjp1) :: orog, rugo, psol … … 61 61 REAL, DIMENSION(iip1,jjm ,llm) :: vvent 62 62 REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: q3d 63 REAL, DIMENSION(klon,nbsrf) :: qsolsrf, snsrf, evap 64 REAL, DIMENSION(klon,nbsrf) :: frugs, agesno 63 REAL, DIMENSION(klon,nbsrf) :: qsolsrf, snsrf 65 64 REAL, DIMENSION(klon,nsoilmx,nbsrf) :: tsoil 66 65 … … 480 479 DO i=1,nbsrf; ftsol(:,i) = tsol; END DO 481 480 DO i=1,nbsrf; snsrf(:,i) = sn; END DO 482 falb1(:,is_ter) = 0.08; falb1(:,is_lic) = 0.6483 falb1(:,is_oce) = 0.5; falb1(:,is_sic) = 0.6484 falb2 = falb1485 481 !albedo SB >>> 486 482 falb_dir(:,is_ter,:)=0.08; falb_dir(:,is_lic,:)=0.6 487 483 falb_dir(:,is_oce,:)=0.5; falb_dir(:,is_sic,:)=0.6 488 484 !albedo SB <<< 489 evap(:,:) = 0.485 fevap(:,:) = 0. 490 486 DO i=1,nbsrf; qsolsrf(:,i)=150.; END DO 491 487 DO i=1,nbsrf; DO j=1,nsoilmx; tsoil(:,j,i) = tsol; END DO; END DO … … 495 491 q_ancien = 0. 496 492 agesno = 0. 497 frugs(:,is_oce) = rugmer(:) 498 frugs(:,is_ter) = MAX(1.0e-05,zstd(:)*zsig(:)/2.0) 499 frugs(:,is_lic) = MAX(1.0e-05,zstd(:)*zsig(:)/2.0) 500 frugs(:,is_sic) = 0.001 493 494 z0m(:,is_oce) = rugmer(:) 495 z0m(:,is_ter) = MAX(1.0e-05,zstd(:)*zsig(:)/2.0) 496 z0m(:,is_lic) = MAX(1.0e-05,zstd(:)*zsig(:)/2.0) 497 z0m(:,is_sic) = 0.001 498 z0h(:,:)=z0m(:,:) 499 501 500 fder = 0.0 502 501 clwcon = 0.0 … … 526 525 527 526 CALL fonte_neige_init(run_off_lic_0) 528 CALL pbl_surface_init( qsol, fder, snsrf, qsolsrf, evap, frugs, agesno, tsoil )527 CALL pbl_surface_init( fder, snsrf, qsolsrf, tsoil ) 529 528 CALL phyredem( "startphy.nc" ) 530 529
Note: See TracChangeset
for help on using the changeset viewer.