Changeset 374 for trunk/LMDZ.GENERIC/libf/phystd/setspv.F90

Timestamp:

Nov 10, 2011, 3:08:51 PM (13 years ago)

Author:

emillour

Message:

Generic GCM: Upgrade: The location of the 'datagcm' directory can now be given in the callphys.def file ( datadir = /absolute/path/to/datagcm ). Changed "datafile.h" into a F90 module "datafile_mod.F90" and spread this change to all routines that used to use "datafile.h".
EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/setspv.F90 (modified) (8 diffs)

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

@@ -26 +26 @@
       use radcommon_h, only: BWNV,BLAMV,WNOV,DWNV,WAVEV, &
                              STELLARF,TAURAY,gastype
+      use datafile_mod, only: datadir
 
       implicit none
@@ -31 +32 @@
 #include "comcstfi.h"
 #include "callkeys.h"
-#include "datafile.h"
 
       logical file_ok
@@ -52 +52 @@
       write(temp1,'(i2.2)') L_NSPECTV
       file_id='/corrk_data/'//trim(adjustl(banddir))//'/narrowbands_VI.in' 
-      file_path=datafile(1:LEN_TRIM(datafile))//file_id(1:LEN_TRIM(file_id))
+      file_path=TRIM(datadir)//TRIM(file_id)
 
       ! check that the file exists
       inquire(FILE=file_path,EXIST=file_ok)
       if(.not.file_ok) then
-         write(*,*)'The file ',file_path(1:LEN_TRIM(file_path))
+         write(*,*)'The file ',TRIM(file_path)
          write(*,*)'was not found by setspv.F90, exiting.'
+         write(*,*)'Check that your path to datagcm:',trim(datadir)
+         write(*,*)' is correct. You can change it in callphys.def with:'
+         write(*,*)' datadir = /absolute/path/to/datagcm'
+         write(*,*)'Also check that the corrkdir you chose in callphys.def exists.'
          call abort
       endif
@@ -70 +74 @@
 
       if(file_entries-1.ne.L_NSPECTV) then
-         write(*,*) 'L_NSPECTV = ',L_NSPECTV, 'in the model, but there are '
+         write(*,*) 'setspv: L_NSPECTV = ',L_NSPECTV, 'in the model, but there are '
          write(*,*) file_entries-1,'entries in ', &
-              file_path(1:LEN_TRIM(file_path)),', exiting.'
+              TRIM(file_path),', exiting.'
          call abort
       endif
@@ -88 +92 @@
 
 
-      print*,'VI band limits:'
+      print*,'setspv: VI band limits:'
       do M=1,L_NSPECTV+1
          print*,m,'-->',BWNV(M),' cm^-1'
@@ -109 +113 @@
 !     Set up stellar spectrum
 
-      write(*,*)'Interpolating stellar spectrum from the hires data...'
+      write(*,*)'setspv: Interpolating stellar spectrum from the hires data...'
       call ave_stelspec(STELLAR)
 
@@ -118 +122 @@
          sum         = sum+STELLARF(N)
       end do
-      write(6,'(" Stellar flux at 1 AU = ",f7.2," W m-2")') sum
+      write(6,'("setspv: Stellar flux at 1 AU = ",f7.2," W m-2")') sum
       print*,' '
 
@@ -130 +134 @@
          call calc_rayleigh
       else
-         print*,'No Rayleigh scattering, check for NaN in output!'
+         print*,'setspv: No Rayleigh scattering, check for NaN in output!'
          do N=1,L_NSPECTV
             TAURAY(N) = 1E-16