Changeset 1545 for trunk/LMDZ.TITAN/libf/phytitan/physiq.F

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.TITAN/libf/phytitan/physiq.F (modified) (7 diffs)

trunk/LMDZ.TITAN/libf/phytitan/physiq.F

@@ -58 +58 @@
       USE infotrac
       use dimphy
-      USE geometry_mod, ONLY: longitude, latitude, cell_area, dx, dy
+      USE geometry_mod, ONLY: longitude, latitude, ! in radians
+     &                        longitude_deg, latitude_deg, ! in degrees
+     &                        cell_area, dx, dy
       use cpdet_mod, only: cpdet, t2tpot
       USE mod_phys_lmdz_para, only : is_parallel,jj_nb
@@ -861 +863 @@
 
 c dans zenang, Ls en degres ; dans mucorr, Ls en radians
-      call mucorr(klon,zls,latitude,rmu0bar,fractbar)
+      call mucorr(klon,zls,latitude_deg,rmu0bar,fractbar)
       IF (cycle_diurne) THEN
         zdtime=dtime*REAL(radpas) ! pas de temps du rayonnement (s)
-        CALL zenang(zlsdeg,gmtime,zdtime,latitude,longitude,rmu0,fract)
+        CALL zenang(zlsdeg,gmtime,zdtime,latitude_deg,longitude_deg,
+     &              rmu0,fract)
       ELSE
         rmu0  = rmu0bar
@@ -917 +920 @@
      $            paprs,pplay,ppk,radsol,falbe,
      e            solsw, sollw, sollwdown, fder,
-     e            longitude, latitude, dx, dy,   
+     e            longitude_deg, latitude_deg, dx, dy,   
      e            debut, lafin,
      s            d_t_vdf,d_u_vdf,d_v_vdf,d_ts,
@@ -1368 +1371 @@
 c A ADAPTER POUR VENUS ET TITAN!!!
 c            CALL lift_noro(klon,klev,dtime,paprs,pplay,
-c     e                   latitude,zmea,zstd,zpic,zgam,zthe,zpic,zval,
+c     e                   latitude_deg,zmea,zstd,zpic,zgam,zthe,zpic,zval,
 c     e                   igwd,idx,itest,
 c     e                   t_seri, u_seri, v_seri,
@@ -1425 +1428 @@
 c====================================================================
       if (ballons.eq.1) then
-         CALL ballon(30,pdtphys,rjourvrai,gmtime,latitude,longitude,
+         CALL ballon(30,pdtphys,rjourvrai,gmtime,
+     &               latitude_deg,longitude_deg,
 c    C               t,pplay,u,v,pphi)   ! alt above surface (smoothed for GCM)
      C               t,pplay,u,v,zphi)   ! alt above planet average radius
@@ -1459 +1463 @@
       CALL aaam_bud (27,klon,klev,rjourvrai,gmtime,
      C               ra,rg,romega,
-     C               latitude,longitude,pphis,
+     C               latitude_deg,longitude_deg,pphis,
      C               zustrdr,zustrli,zustrcl,
      C               zvstrdr,zvstrli,zvstrcl,
@@ -1545 +1549 @@
 c     DO k = 1, klev
 c     DO i = 1, klon
-c       mang(i,k) = RA*cos(latitude(i)*RPI/180.)
-c    .     *(u_seri(i,k)+RA*cos(latitude(i)*RPI/180.)*ROMEGA)
+c       mang(i,k) = RA*cos(latitude(i))
+c    .     *(u_seri(i,k)+RA*cos(latitude(i))*ROMEGA)
 c    .     *cell_area(i)*(paprs(i,k)-paprs(i,k+1))/RG
 c       mangtot=mangtot+mang(i,k)