Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.TITAN/libf/phytitan/pg2.old

    r1543 r1545  
    4242
    4343      use dimphy
    44       USE geometry_mod, ONLY: latitude
     44      USE geometry_mod, ONLY: latitude_deg
    4545#include "dimensions.h"
    4646#include "microtab.h"
     
    180180      lati(1)     = 0.5*RPI
    181181      DO ig=2,ngrid-1
    182         lati(ig)  = latitude(2+(ig-2)*iim)*RPI/180.
     182        lati(ig)  = latitude_deg(2+(ig-2)*iim)*RPI/180.
    183183      ENDDO
    184184      lati(ngrid) = -0.5*RPI
     
    197197c     print*,"ENTREE PG2 PREMIER APPEL"
    198198c     print*,airetot,' airetot?= ',4.*RPI*RA*RA
    199 c     print*,1,latitude(1),aire(1),aire(1)/airetot,' aires'
     199c     print*,1,latitude_deg(1),aire(1),aire(1)/airetot,' aires'
    200200c     DO ig=2,ngrid-1
    201 c     print*,ig,latitude(2+(ig-2)*iim),aire(ig),aire(ig)/airetot,' aires'
     201c     print*,ig,latitude_deg(2+(ig-2)*iim),aire(ig),aire(ig)/airetot,' aires'
    202202c     ENDDO
    203 c     print*,ngrid,latitude(klon),aire(ngrid),aire(ngrid)/airetot,' aires'
     203c     print*,ngrid,latitude_deg(klon),aire(ngrid),aire(ngrid)/airetot,' aires'
    204204c     stop
    205205     
Note: See TracChangeset for help on using the changeset viewer.