Changeset 699 for trunk/LMDZ.GENERIC/libf/phystd/phyetat0.F

Timestamp:

Jun 6, 2012, 12:47:56 PM (12 years ago)

Author:

emillour

Message:

Generic GCM:

• Corrected the polar mesh surface area which was wrong in the physics (changes in phyetat0.F, calfis.F and newstart.F)
• Some cleanup in newstart.F (removed some obsolete "Mars" options: mons_ice,..) and also added option "q=profile" to initialize a tracer with a profile read from file "profile_tracername"

EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/phyetat0.F (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/phystd/phyetat0.F

@@ -11 +11 @@
 #include "dimensions.h"
 #include "dimphys.h"
-#include "comgeomfi.h"
+!#include "comgeomfi.h"
 #include "surfdat.h"
 #include "planete.h"
@@ -111 +111 @@
      .              p_omeg,p_g,p_cpp,p_mugaz,p_daysec,time)
 c
-c Lecture des latitudes (coordonnees):
-c
-      ierr = NF_INQ_VARID (nid, "latitude", nvarid)
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, 'phyetat0: Le champ <latitude> est absent'
-         CALL abort
-      ENDIF
-#ifdef NC_DOUBLE
-      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, lati)
-#else
-      ierr = NF_GET_VAR_REAL(nid, nvarid, lati)
-#endif
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, 'phyetat0: Lecture echouee pour <latitude>'
-         CALL abort
-      ENDIF
-c
-c Lecture des longitudes (coordonnees):
-c
-      ierr = NF_INQ_VARID (nid, "longitude", nvarid)
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, 'phyetat0: Le champ <longitude> est absent'
-         CALL abort
-      ENDIF
-#ifdef NC_DOUBLE
-      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, long)
-#else
-      ierr = NF_GET_VAR_REAL(nid, nvarid, long)
-#endif
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, 'phyetat0: Lecture echouee pour <longitude>'
-         CALL abort
-      ENDIF
-c
-c Lecture des aires des mailles:
-c
-      ierr = NF_INQ_VARID (nid, "area", nvarid)
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, 'phyetat0: Le champ <area> est absent'
-         CALL abort
-      ENDIF
-#ifdef NC_DOUBLE
-      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, area)
-#else
-      ierr = NF_GET_VAR_REAL(nid, nvarid, area)
-#endif
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, 'phyetat0: Lecture echouee pour <area>'
-         CALL abort
-      ENDIF
-      xmin = 1.0E+20
-      xmax = -1.0E+20
-      xmin = MINVAL(area)
-      xmax = MAXVAL(area)
-      PRINT*,'Aires des mailles <area>:', xmin, xmax
+c Read latitudes (coordinates): No need, these are provided by the dynamics
+c
+!      ierr=nf90_inq_varid(nid,"latitude",nvarid)
+!      IF (ierr.NE.nf90_noerr) THEN
+!         PRINT*, 'phyetat0: Le champ <latitude> est absent'
+!         write(*,*)trim(nf90_strerror(ierr))
+!         CALL abort
+!      ENDIF
+!      ierr=nf90_get_var(nid,nvarid,lati)
+!      IF (ierr.NE.nf90_noerr) THEN
+!         PRINT*, 'phyetat0: Lecture echouee pour <latitude>'
+!         write(*,*)trim(nf90_strerror(ierr))
+!         CALL abort
+!      ENDIF
+c
+c read longitudes (coordinates): No need, these are provided by the dynamics
+c
+!      ierr=nf90_inq_varid(nid,"longitude",nvarid)
+!      IF (ierr.NE.nf90_noerr) THEN
+!         PRINT*, 'phyetat0: Le champ <longitude> est absent'
+!         write(*,*)trim(nf90_strerror(ierr))
+!         CALL abort
+!      ENDIF
+!      ierr=nf90_get_var(nid,nvarid,long)
+!      IF (ierr.NE.nf90_noerr) THEN
+!         PRINT*, 'phyetat0: Lecture echouee pour <longitude>'
+!         write(*,*)trim(nf90_strerror(ierr))
+!         CALL abort
+!      ENDIF
+c
+c Read areas of meshes: No need, these are provided by the dynamics
+c
+!      ierr=nf90_inq_varid(nid,"area",nvarid)
+!      IF (ierr.NE.nf90_noerr) THEN
+!         PRINT*, 'phyetat0: Le champ <area> est absent'
+!         write(*,*)trim(nf90_strerror(ierr))
+!         CALL abort
+!      ENDIF
+!      ierr=nf90_get_var(nid,nvarid,area)
+!      IF (ierr.NE.nf90_noerr) THEN
+!         PRINT*, 'phyetat0: Lecture echouee pour <area>'
+!         write(*,*)trim(nf90_strerror(ierr))
+!         CALL abort
+!      ENDIF
+!      xmin = 1.0E+20
+!      xmax = -1.0E+20
+!      xmin = MINVAL(area)
+!      xmax = MAXVAL(area)
+!      PRINT*,'Aires des mailles <area>:', xmin, xmax
 c
 c Lecture du geopotentiel au sol: