Changeset 2395 for LMDZ5/trunk/libf/phylmd

Timestamp:

Nov 18, 2015, 1:25:35 PM (9 years ago)

Author:

Ehouarn Millour

Message:

Bug fix for aquaplanets: "rlat" and "rlon" were uninitialized when written to startphy.nc; phyredem should write "longitude_deg" and "latitude_deg" to the file. Overall "rlat" and "rlon" should not be used and "latitude_deg" and "longitude_deg" should be used instead in the physics. To be further cleaned up.
Some cleanup on that mater also made in phydev.
Note that this change will make bench test results different for longitudes and latitudes, because of roundoff effects.
EM

File:

• LMDZ5/trunk/libf/phylmd/phyredem.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/phyredem.F90

@@ -6 +6 @@
 ! Purpose: Write restart state for physics.
 !-------------------------------------------------------------------------------
-  USE dimphy
-  USE mod_grid_phy_lmdz
-  USE mod_phys_lmdz_para
+  USE dimphy, ONLY: klon, klev
   USE fonte_neige_mod,  ONLY : fonte_neige_final
   USE pbl_surface_mod,  ONLY : pbl_surface_final
-  USE phys_state_var_mod
-  USE iostart
+  USE phys_state_var_mod, ONLY: radpas, zmasq, pctsrf, ftsol, falb_dir, &
+                                falb_dif, qsol, fevap, radsol, solsw, sollw, &
+                                sollwdown, rain_fall, snow_fall, z0m, z0h, &
+                                agesno, zmea, zstd, zsig, zgam, zthe, zpic, &
+                                zval, rugoro, t_ancien, q_ancien, u_ancien, &
+                                v_ancien, clwcon, rnebcon, ratqs, pbl_tke, &
+                                wake_delta_pbl_tke, zmax0, f0, sig1, w01, &
+                                wake_deltat, wake_deltaq, wake_s, wake_cstar, &
+                                wake_pe, wake_fip, fm_therm, entr_therm, &
+                                detr_therm, Ale_bl, Ale_bl_trig, Alp_bl, &
+                                du_gwd_rando, du_gwd_front
+  USE geometry_mod, ONLY : longitude_deg, latitude_deg
+  USE iostart, ONLY: open_restartphy, close_restartphy, put_field, put_var
   USE traclmdz_mod, ONLY : traclmdz_to_restart
   USE infotrac_phy, ONLY: type_trac, niadv, tname, nbtr, nqo
   USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, co2_send
-  USE indice_sol_mod
-  USE surface_data
+  USE indice_sol_mod, ONLY: nbsrf, is_oce, is_sic, is_ter, is_lic, epsfra
+  USE surface_data, ONLY: type_ocean, version_ocean
   USE ocean_slab_mod, ONLY : tslab, seaice, tice, fsic
   USE time_phylmdz_mod, ONLY: annee_ref, day_end, itau_phy, pdtphys
@@ -90 +99 @@
 
   CALL put_field("longitude", &
-       "Longitudes de la grille physique", rlon)
-
-  CALL put_field("latitude", "Latitudes de la grille physique", rlat)
+       "Longitudes de la grille physique", longitude_deg)
+
+  CALL put_field("latitude", "Latitudes de la grille physique", latitude_deg)
 
   ! PB ajout du masque terre/mer