Changeset 2359 for LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Sep 4, 2015, 6:48:31 PM (9 years ago)

Author:

oboucher

Message:

Recoded diurnal cycle in physiq.F90
iflag_diurnal_cycle is the new flag
=0 for no diurnal cycle
=1 for old diurnal cycle
=2 for new diurnal cycle with updates
in between two calls to radiation according
to how rmu0 has changed over the previous
timestep of physiq

File:

• LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -97 +97 @@
 !$OMP THREADPRIVATE(clwcon0th,rnebcon0th)
 ! radiation outputs
+      REAL,ALLOCATABLE,SAVE :: zswdn0(:,:), zswdn(:,:)
+!$OMP THREADPRIVATE(zswdn0,zswdn)
+      REAL,ALLOCATABLE,SAVE :: zswup0(:,:), zswup(:,:)
+!$OMP THREADPRIVATE(zswup0,zswup)
       REAL,ALLOCATABLE,SAVE :: swdn0(:,:), swdn(:,:)
 !$OMP THREADPRIVATE(swdn0,swdn)
@@ -340 +344 @@
       REAL,ALLOCATABLE,SAVE :: lwup0p(:,:), lwupp(:,:)
 !$OMP THREADPRIVATE(lwdn0p, lwdnp, lwup0p, lwupp)
+      REAL,ALLOCATABLE,SAVE :: zswdn0p(:,:), zswdnp(:,:)
+      REAL,ALLOCATABLE,SAVE :: zswup0p(:,:), zswupp(:,:)
+!$OMP THREADPRIVATE(zswdn0p, zswdnp, zswup0p, zswupp)
       REAL,ALLOCATABLE,SAVE :: swdn0p(:,:), swdnp(:,:)
       REAL,ALLOCATABLE,SAVE :: swup0p(:,:), swupp(:,:)
@@ -458 +465 @@
       ALLOCATE(clwcon0th(klon,klev),rnebcon0th(klon,klev))
 ! radiation outputs
+      ALLOCATE(zswdn0(klon,klevp1), zswdn(klon,klevp1))
+      ALLOCATE(zswup0(klon,klevp1), zswup(klon,klevp1))
       ALLOCATE(swdn0(klon,klevp1), swdn(klon,klevp1))
       ALLOCATE(swup0(klon,klevp1), swup(klon,klevp1))
@@ -557 +566 @@
       ALLOCATE(lwdn0p(klon,klevp1), lwdnp(klon,klevp1))
       ALLOCATE(lwup0p(klon,klevp1), lwupp(klon,klevp1))
+      ALLOCATE(zswdn0p(klon,klevp1), zswdnp(klon,klevp1))
+      ALLOCATE(zswup0p(klon,klevp1), zswupp(klon,klevp1))
       ALLOCATE(swdn0p(klon,klevp1), swdnp(klon,klevp1))
       ALLOCATE(swup0p(klon,klevp1), swupp(klon,klevp1))
@@ -611 +622 @@
       deallocate(clwcon0th, rnebcon0th)
 ! radiation outputs
+      deallocate(zswdn0, zswdn)
+      deallocate(zswup0, zswup)
       deallocate(swdn0, swdn)
       deallocate(swup0, swup)
@@ -687 +700 @@
       deallocate(lwdn0p, lwdnp)
       deallocate(lwup0p, lwupp)
+      deallocate(zswdn0p, zswdnp)
+      deallocate(zswup0p, zswupp)
       deallocate(swdn0p, swdnp)
       deallocate(swup0p, swupp)