Changeset 1543 for trunk/LMDZ.TITAN/libf/phytitan/phyetat0.F90

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.TITAN/libf/phytitan/phyetat0.F90 (modified) (5 diffs)

trunk/LMDZ.TITAN/libf/phytitan/phyetat0.F90

@@ -13 +13 @@
       USE iostart
       USE infotrac
-      USE comgeomphy,  only: rlatd,rlond
+      USE geometry_mod,  only: latitude,longitude
       USE time_phylmdz_mod, only: itau_phy, raz_date
 
@@ -76 +76 @@
 
 ! read latitudes
-call get_field("latitude",rlatd,found)
+call get_field("latitude",latitude,found)
       IF (.not.found) THEN
          PRINT*, 'phyetat0: Le champ <latitude> est absent'
@@ -83 +83 @@
 
 ! read longitudes
-call get_field("longitude",rlond,found)
+call get_field("longitude",longitude,found)
       IF (.not.found) THEN
          PRINT*, 'phyetat0: Le champ <longitude> est absent'
@@ -319 +319 @@
       resch4(1) = 0.    ! pole nord = 1 point
       DO i=2,klon
-          if ((rlatd(i).ge.75..and.rlatd(i).le.85.).or.  &
-              (rlatd(i).ge.-85.and.rlatd(i).le.-75.)) then
+          if ((latitude(i).ge.75..and.latitude(i).le.85.).or.  &
+              (latitude(i).ge.-85.and.latitude(i).le.-75.)) then
             resch4(i) = 2.
           else
@@ -352 +352 @@
 
 ! do some more initializations
-call init_iophy_new(rlatd,rlond)
+call init_iophy_new(latitude,longitude)
 
 end subroutine phyetat0