Changeset 2410 for LMDZ5/trunk/libf/phylmd/surf_land_orchidee_mod.F90

Timestamp:

Dec 16, 2015, 9:10:00 AM (8 years ago)

Author:

jghattas

Message:

Added the variable yrmu0(cosine of the solar zenith angle) into ORCHIDEE:

Current version of LMDZ can now be used with ORCHIDEE trunk revisions from rev 2961 and newer. For revision 1078-2960 on ORCHIDEE trunk, a small modification (change coszang into sinang) in surf_land_orchidee_mod.f90 is needed. For older versions than ORCHIDEE trunk revision 1078, the interface in surf_land_orchidee_noopenmp_mod.f90 should be used (add cpp key ORCHIDEE_NOOPENMP).

For details see ticket ​https://forge.ipsl.jussieu.fr/orchidee/ticket/217

File:

• LMDZ5/trunk/libf/phylmd/surf_land_orchidee_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/surf_land_orchidee_mod.F90

@@ -3 +3 @@
 #ifndef ORCHIDEE_NOOPENMP
 !
-! This module controles the interface towards the model ORCHIDEE
+! This module controles the interface towards the model ORCHIDEE.
+!
+! Compatibility with ORCHIDIEE :
+! The current version can be used with ORCHIDEE/trunk from revision 2961.
+! This interface can also be used with ORCHIDEE/trunk revision 1078-2960 if changing 
+! coszang=yrmu0 into sinang=yrmu0 at 2 places later below in this module. 
 !
 ! Subroutines in this module : surf_land_orchidee
@@ -31 +36 @@
 !  
   SUBROUTINE surf_land_orchidee(itime, dtime, date0, knon, &
-       knindex, rlon, rlat, pctsrf, &
+       knindex, rlon, rlat, yrmu0, pctsrf, &
        debut, lafin, &
        plev,  u1_lay, v1_lay, gustiness, temp_air, spechum, epot_air, ccanopy, & 
@@ -115 +120 @@
     REAL, DIMENSION(klon,nbsrf), INTENT(IN)   :: pctsrf
     REAL, DIMENSION(klon), INTENT(IN)         :: rlon, rlat
+    REAL, DIMENSION(klon), INTENT(IN)         :: yrmu0 ! cosine of solar zenith angle
     REAL, DIMENSION(klon), INTENT(IN)         :: plev
     REAL, DIMENSION(klon), INTENT(IN)         :: u1_lay, v1_lay, gustiness
@@ -405 +411 @@
                evap, fluxsens, fluxlat, coastalflow, riverflow, &
                tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0_new, &
-               lon_scat, lat_scat, q2m, t2m)
+               lon_scat, lat_scat, q2m, t2m, coszang=yrmu0)
 #endif         
        ENDIF
@@ -429 +435 @@
             evap(1:knon), fluxsens(1:knon), fluxlat(1:knon), coastalflow(1:knon), riverflow(1:knon), &
             tsol_rad(1:knon), tsurf_new(1:knon), qsurf(1:knon), albedo_out(1:knon,:), emis_new(1:knon), z0_new(1:knon), &
-            lon_scat, lat_scat, q2m, t2m)
+            lon_scat, lat_scat, q2m, t2m, coszang=yrmu0(1:knon))
 #endif       
     ENDIF