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

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

File:

• trunk/LMDZ.GENERIC/libf/phystd/suaer_corrk.F90 (modified) (3 diffs)

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