Changeset 146 for trunk/LMDZ.MARS/libf/phymars/datareadnc.F

Timestamp:

Jun 8, 2011, 12:04:26 PM (13 years ago)

Author:

emillour

Message:

Mars GCM:

minor bug fix in lect_start_archive.F (using wrong surface temp. array).

+ swiched output messages to english and added that tracers not found
in file must be initialized by user.

minor bug fix in datareadnc.F : 'datafile' path must be initialized.

EM

File:

• trunk/LMDZ.MARS/libf/phymars/datareadnc.F (modified) (6 diffs)

trunk/LMDZ.MARS/libf/phymars/datareadnc.F

@@ -41 +41 @@
 c                   (rlatvd de 89 a -89)
 c=======================================================================
+
+! to use  'getin'
+       use ioipsl_getincom 
 
       implicit none
@@ -77 +80 @@
       REAL        zthe(imdp1*jmdp1)
 
-      INTEGER   lnblnk, ierr
-      EXTERNAL    lnblnk
+      INTEGER     ierr
 
       INTEGER   unit,nvarid
@@ -95 +97 @@
       DIMENSION string(4)
 
-      CHARACTER*80    file
 
 #include "lmdstd.h"
@@ -109 +110 @@
 c    Lecture NetCDF des donnees latitude et longitude
 c-----------------------------------------------------------------------
-      write(*,*) 'ouverture du fichier surface.nc'
-
-      ierr = NF_OPEN (datafile(1:lnblnk(datafile))//'/surface.nc',
+      write(*,*) 'datareadnc: opening file surface.nc'
+
+      datafile="/u/forget/WWW/datagcm/datafile" ! default path to surface.nc
+      call getin("datadir",datafile) ! but users may specify another path
+      
+      ierr = NF_OPEN (trim(datafile)//'/surface.nc',
      &  NF_NOWRITE,unit)
       IF (ierr.NE.NF_NOERR) THEN
         write(*,*)'Error : cannot open file surface.nc '
         write(*,*)'(in phymars/datareadnc.F)'
-        write(*,*)'It should be in :',datafile(1:lnblnk(datafile)),'/'
-        write(*,*)'1) You can change this directory address in '
-        write(*,*)'   file phymars/datafile.h'
+        write(*,*)'It should be in :',trim(datafile),'/'
+        write(*,*)'1) You can set this path in the callphys.def file:'
+        write(*,*)'   datadir=/path/to/the/datafiles'
         write(*,*)'2) If necessary surface.nc (and other datafiles)'
         write(*,*)'   can be obtained online on:'
@@ -176 +180 @@
       string(2) = 'thermal'
       if (relief.ne.'pla') then
-        write(*,*) ' La topographie est celle de MOLA'
+        write(*,*) ' MOLA topography'
         relief = 'MOL'
           string(3) = 'z'//relief
@@ -282 +286 @@
             phisinit(i)=1000.*phisinit(i)
       ENDDO
-      CALL dump2d(iimp1,jjp1,phisinit,'Altitude en m')
+      CALL dump2d(iimp1,jjp1,phisinit,'Altitude in m')
       CALL multscal(iimp1*jjp1,phisinit,g,phisinit)