Changeset 1545 for trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.NewtonCool

Timestamp:

May 2, 2016, 8:31:27 AM (9 years ago)

Author:

emillour

Message:

Venus and Titan GCMs:

Adaptation wrt previous changes for Titan and Venus where
longitude and latitude arrays (in phycommon/geometry_mod) were overwritten
with values from startphy.nc files, where values are given in degrees.
For the sake of homegeneity with other physics package, revert to "default"
behaviour: longitude/latitude are in radians and longitude_deg/latitude_deg
are in degrees.
Also added checking in phyetat0 that the longitude/latitude read in the
restartphy.nc files match the ones provided by the dynamics.

EM

File:

• trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.NewtonCool (modified) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.NewtonCool

@@ -24 +24 @@
 c======================================================================
       use dimphy
-      USE geometry_mod, ONLY: latitude
+      USE geometry_mod, ONLY: latitude ! in radians
       USE phys_state_var_mod, only: heat,cool,radsol,
      .      topsw,toplw,solsw,sollw,sollwdown,lwnet,swnet
@@ -98 +98 @@
             ztemp = tempCLee(level)*(1-fact)+tempCLee(level+1)*fact
             zdt   = dt_epCLee(level)*(1-fact)+dt_epCLee(level+1)*fact
-c           zt_eq(i,k) = ztemp + zdt*(cos(latitude(i)*RPI/180.)-2./RPI)
-            zt_eq(i,k) = ztemp + zdt*(cos(latitude(i)*RPI/180.)-RPI/4.)
+c           zt_eq(i,k) = ztemp + zdt*(cos(latitude(i))-2./RPI)
+            zt_eq(i,k) = ztemp + zdt*(cos(latitude(i))-RPI/4.)
             
           enddo