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

Timestamp:

Apr 22, 2016, 9:02:11 AM (9 years ago)

Author:

emillour

Message:

All models: Further adaptations to keep up with changes in LMDZ5 concerning
physics/dynamics separation:

• dyn3d:
• adapted gcm.F so that all physics initializations are now done in iniphysiq.

• dyn3dpar:
• adapted gcm.F so that all physics initializations are now done in iniphysiq.
• updated calfis_p.F to follow up with changes.
• copied over updated "bands.F90" from LMDZ5.

• dynphy_lonlat:
• calfis_p.F90, mod_interface_dyn_phys.F90, follow up of changes in phy_common/mod_* routines

• phy_common:
• added "geometry_mod.F90" to store information about the grid (replaces phy*/comgeomphy.F90) and give variables friendlier names: rlond => longitude , rlatd => latitude, airephy => cell_area, cuphy => dx , cvphy => dy
• added "physics_distribution_mod.F90"
• updated "mod_grid_phy_lmdz.F90", "mod_phys_lmdz_mpi_data.F90", "mod_phys_lmdz_para.F90", "mod_phys_lmdz_mpi_transfert.F90", "mod_grid_phy_lmdz.F90", "mod_phys_lmdz_omp_data.F90", "mod_phys_lmdz_omp_transfert.F90", "write_field_phy.F90" and "ioipsl_getin_p_mod.F90" to LMDZ5 versions.

• phy[venus/titan/mars/std]:
• removed "init_phys_lmdz.F90", "comgeomphy.F90"; adapted routines to use geometry_mod (longitude, latitude, cell_area, etc.)

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 comgeomphy
+      USE geometry_mod, ONLY: latitude
       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(rlatd(i)*RPI/180.)-2./RPI)
-            zt_eq(i,k) = ztemp + zdt*(cos(rlatd(i)*RPI/180.)-RPI/4.)
+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.)
             
           enddo