Ignore:
Timestamp:
Oct 16, 2016, 11:17:32 AM (8 years ago)
Author:
oboucher
Message:

adding new diagnostics for evap/rain/snow weighted by surface types
this parallels what was done for wbilo (net flux)
new diagnostics are wevap_??? wrain_??? wsnow_???
where ??? is ter, oce, lic or sic.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

    r2635 r2670  
    372372      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: fsolsw, wfbils, wfbilo
    373373!$OMP THREADPRIVATE(fsolsw, wfbils, wfbilo)
     374      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: wfevap, wfrain, wfsnow
     375!$OMP THREADPRIVATE(wfevap,wfrain,wfsnow)
    374376      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: t2m, fluxlat, fsollw,evap_pot
    375377!$OMP THREADPRIVATE(t2m, fluxlat, fsollw,evap_pot)
    … …  
    633635      ALLOCATE(fsollw(klon, nbsrf))
    634636      ALLOCATE(fsolsw(klon, nbsrf), wfbils(klon, nbsrf), wfbilo(klon, nbsrf))
     637      ALLOCATE(wfevap(klon, nbsrf), wfrain(klon,nbsrf), wfsnow(klon, nbsrf))
    635638      ALLOCATE(evap_pot(klon, nbsrf))
    636639
    … …  
    854857      DEALLOCATE(fsollw, evap_pot)
    855858      DEALLOCATE(fsolsw, wfbils, wfbilo)
     859      DEALLOCATE(wfevap,wfrain,wfsnow)
    856860
    857861      DEALLOCATE(pmflxr, pmflxs, prfl)
Note: See TracChangeset for help on using the changeset viewer.