Changeset 2351 for LMDZ5/trunk/libf/phylmd/phyaqua_mod.F90

Timestamp:

Aug 25, 2015, 5:14:59 PM (9 years ago)

Author:

Ehouarn Millour

Message:

More on physics/dynamics separation and cleanup:

• Set things up so that all physics-related initializations are done via iniphysiq.
• Created a "geometry_mod.F90" module in phy_common to store information on the loacl grid (i.e. replaces comgeomphy) and moreover give these variables more obvious names (e.g.: rlond => longitude, rlatd => latitude, airephy => cell_area).
• removed obsolete comgeomphy.h and comgeomphy.F90

EM

File:

• LMDZ5/trunk/libf/phylmd/phyaqua_mod.F90 (modified) (7 diffs)

LMDZ5/trunk/libf/phylmd/phyaqua_mod.F90

@@ -5 +5 @@
 CONTAINS
 
-  SUBROUTINE iniaqua(nlon, latfi, lonfi, iflag_phys)
+  SUBROUTINE iniaqua(nlon, iflag_phys)
 
     ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -19 +19 @@
     ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-    USE comgeomphy, ONLY: rlatd, rlond
     USE dimphy, ONLY: klon
+    USE geometry_mod, ONLY : latitude
     USE surface_data, ONLY: type_ocean, ok_veget
     USE pbl_surface_mod, ONLY: pbl_surface_init
@@ -38 +38 @@
     INTEGER, INTENT (IN) :: nlon, iflag_phys
     ! IM ajout latfi, lonfi
-    REAL, INTENT (IN) :: lonfi(nlon), latfi(nlon)
+!    REAL, INTENT (IN) :: lonfi(nlon), latfi(nlon)
 
     INTEGER type_profil, type_aqua
@@ -117 +117 @@
     INTEGER l, ierr, aslun
 
-    REAL longitude, latitude
+!    REAL longitude, latitude
     REAL paire
 
-    DATA latitude, longitude/48., 0./
+!    DATA latitude, longitude/48., 0./
 
     ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -149 +149 @@
 
     ! IM ajout latfi, lonfi
-    rlatd = latfi
-    rlond = lonfi
-    rlat = rlatd*180./pi
-    rlon = rlond*180./pi
+!    rlatd = latfi
+!    rlond = lonfi
+!    rlat = rlatd*180./pi
+!    rlon = rlond*180./pi
 
     ! -----------------------------------------------------------------------
@@ -160 +160 @@
     day_ini = day_ref
     day_end = day_ini + ndays
-    airefi = 1.
-    zcufi = 1.
-    zcvfi = 1.
+!    airefi = 1.
+!    zcufi = 1.
+!    zcvfi = 1.
     !$OMP MASTER
     nbapp_rad_omp = 24
@@ -260 +260 @@
     END DO
     ! IM calcul profil sst
-    CALL profil_sst(nlon, rlatd, type_profil, phy_sst)
+    CALL profil_sst(nlon, latitude, type_profil, phy_sst)
 
     CALL writelim(klon, phy_nat, phy_alb, phy_sst, phy_bil, phy_rug, phy_ice, &