Changeset 4737 for LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90

Timestamp:

Oct 24, 2023, 12:13:53 PM (9 months ago)

Author:

oboucher

Message:

moving the calculation of wfbilo wfrain and wfsnow inside phys_output_write_mod so that it uses the timestep values of rain and snow instead of the values from the previous timestep as it was the case in pbl_surface. Thus there is exact mass balance in the output.

File:

• LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90

@@ -455 +455 @@
       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: proba_notrig, random_notrig
 !$OMP THREADPRIVATE(proba_notrig, random_notrig)
-      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: fsolsw, wfbils, wfbilo
-!$OMP THREADPRIVATE(fsolsw, wfbils, wfbilo)
-      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: wfevap, wfrain, wfsnow
-!$OMP THREADPRIVATE(wfevap,wfrain,wfsnow)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: fsolsw, wfbils
+!$OMP THREADPRIVATE(fsolsw, wfbils)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: wfevap
+!$OMP THREADPRIVATE(wfevap)
       REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: t2m, fluxlat, fsollw,evap_pot
 !$OMP THREADPRIVATE(t2m, fluxlat, fsollw,evap_pot)
@@ -990 +990 @@
       ALLOCATE(t2m(klon, nbsrf), fluxlat(klon, nbsrf))
       ALLOCATE(fsollw(klon, nbsrf))
-      ALLOCATE(fsolsw(klon, nbsrf), wfbils(klon, nbsrf), wfbilo(klon, nbsrf))
-      ALLOCATE(wfevap(klon, nbsrf), wfrain(klon,nbsrf), wfsnow(klon, nbsrf))
+      ALLOCATE(fsolsw(klon, nbsrf), wfbils(klon, nbsrf))
+      ALLOCATE(wfevap(klon, nbsrf))
       ALLOCATE(evap_pot(klon, nbsrf))
 ! FC
@@ -1361 +1361 @@
       DEALLOCATE(t2m, fluxlat)
       DEALLOCATE(fsollw, evap_pot)
-      DEALLOCATE(fsolsw, wfbils, wfbilo)
-      DEALLOCATE(wfevap,wfrain,wfsnow)
+      DEALLOCATE(fsolsw, wfbils)
+      DEALLOCATE(wfevap)
 
       DEALLOCATE(pmflxr, pmflxs)