Changeset 1470 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Sep 16, 2015, 1:00:04 PM (9 years ago)

Author:

emillour

Message:

Generic GCM:

• reorganizing the "datadir" structure: aerosol properties should now be in subdirectory 'aerosol_properties' of datadir, and surface.nc files should be in subdirectory 'surface_data'. These subdirectory names are stored in module datafile_mod.
• Made things retro-compatible so that using an 'old' datadir structure (ie. aerosol properties files and surface files in datadir) still works.

EM

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• datafile_mod.F90 (modified) (2 diffs)
• suaer_corrk.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/datafile_mod.F90

@@ -1 +1 @@
 !-----------------------------------------------------------------------
       module datafile_mod
-!  Address of the directory containing tables of data needed by the UCM
+!  Address of the directory containing tables of data needed by the GCM
       implicit none
 
+      ! Main directory: 'datadir':
       ! Default for Berserker @ UChicago:
 !      character(len=300) :: datadir='/home/rwordsworth/datagcm'
@@ -9 +10 @@
 !      character(len=300) :: datadir='/san/home/rdword/gcm/datagcm'
       ! Default for LMD machines:
-      character(len=300) :: datadir='/u/rwlmd/datagcm'
+      character(len=300) :: datadir='/u/lmdz/WWW/planets/LMDZ.GENERIC/datagcm'
+      
+      ! Subdirectories of 'datadir':
+      
+      ! surfdir stores planetary topography, albedo, etc. (surface.nc files)
+      character(len=12),parameter :: surfdir="surface_data"
+      
+      ! aerdir stores aerosol properties files (optprop_*dat files)
+      character(LEN=18),parameter :: aerdir="aerosol_properties" 
 
       end module datafile_mod

trunk/LMDZ.GENERIC/libf/phystd/suaer_corrk.F90

@@ -4 +4 @@
       use radinc_h,    only: L_NSPECTI,L_NSPECTV,nsizemax,iim,jjm,naerkind
       use radcommon_h, only: blamv,blami,lamrefir,lamrefvis
-      use datafile_mod, only: datadir
+      use datafile_mod, only: datadir, aerdir
 
       ! outputs
@@ -205 +205 @@
 
 !$OMP MASTER
-            INQUIRE(FILE=TRIM(datadir)//&
-         '/'//TRIM(file_id(iaer,idomain)),&
-            EXIST=file_ok)
+            INQUIRE(FILE=TRIM(datadir)//'/'//TRIM(aerdir)//&
+                    '/'//TRIM(file_id(iaer,idomain)),&
+                    EXIST=file_ok)
+            IF (file_ok) THEN
+              OPEN(UNIT=file_unit,&
+                   FILE=TRIM(datadir)//'/'//TRIM(aerdir)//&
+                        '/'//TRIM(file_id(iaer,idomain)),&
+                   FORM='formatted')
+            ELSE
+             ! In ye old days these files were stored in datadir;
+             ! let's be retro-compatible
+              INQUIRE(FILE=TRIM(datadir)//&
+                      '/'//TRIM(file_id(iaer,idomain)),&
+                      EXIST=file_ok)
+              IF (file_ok) THEN
+                OPEN(UNIT=file_unit,&
+                     FILE=TRIM(datadir)//&
+                          '/'//TRIM(file_id(iaer,idomain)),&
+                     FORM='formatted')
+              ENDIF              
+            ENDIF
             IF(.NOT.file_ok) THEN
                write(*,*)'suaer_corrk: Problem opening ',&
                TRIM(file_id(iaer,idomain))
-               write(*,*)'It should be in: ',TRIM(datadir)
+               write(*,*)'It should be in: ',TRIM(datadir)//'/'//TRIM(aerdir)
                write(*,*)'1) You can set this directory address ',&
                'in callphys.def with:'
@@ -220 +238 @@
                write(*,*)' can be obtained online at:'
                write(*,*)' http://www.lmd.jussieu.fr/',&
-               '~forget/datagcm/datafile'
+               '~lmdz/planets/LMDZ.GENERIC/datagcm/'
                CALL ABORT 
             ENDIF
-            OPEN(UNIT=file_unit,&
-            FILE=TRIM(datadir)//&
-         '/'//TRIM(file_id(iaer,idomain)),&
-            FORM='formatted')
 
 !     1.2 Allocate the optical property table