Changeset 1918 for trunk/LMDZ.MARS/libf/phymars/albedocaps.F90
- Timestamp:
- Apr 12, 2018, 8:35:48 AM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/albedocaps.F90
r1779 r1918 90 90 use geometry_mod, only: latitude, longitude ! in radians 91 91 use surfdat_h, only: albedice, TESice_Ncoef, TESice_Scoef 92 use datafile_mod, only: datadir 92 93 use netcdf, only: nf90_open, NF90_NOWRITE, NF90_NOERR, & 93 94 nf90_strerror, nf90_inq_varid, nf90_get_var, nf90_close 94 95 95 96 implicit none 96 include"datafile.h"97 97 98 98 ! arguments: … … 151 151 if (firstcall) then 152 152 ! Load TES albedoes for Northern Hemisphere 153 ! Note: datafile() is defined in "datafile.h" 154 ierr=nf90_open(trim(datafile)//"/npsc_albedo.nc",NF90_NOWRITE,nid) 153 ierr=nf90_open(trim(datadir)//"/npsc_albedo.nc",NF90_NOWRITE,nid) 155 154 IF (ierr.NE.NF90_NOERR) THEN 156 155 write(*,*)'Problem opening npsc_albedo.nc (phymars/albedocaps.F90)' 157 write(*,*)'It should be in :',trim(data file),'/'156 write(*,*)'It should be in :',trim(datadir),'/' 158 157 write(*,*)'1) You can change this directory address in callfis.def with' 159 158 write(*,*)' datadir=/path/to/datafiles' … … 163 162 CALL ABORT 164 163 ELSE 165 write(*,*) "albedocaps: using file ",trim(data file)//"/npsc_albedo.nc"164 write(*,*) "albedocaps: using file ",trim(datadir)//"/npsc_albedo.nc" 166 165 ENDIF 167 166 … … 225 224 226 225 ! Load albedoes for Southern Hemisphere 227 ierr=nf90_open(trim(data file)//"/spsc_albedo.nc",NF90_NOWRITE,nid)226 ierr=nf90_open(trim(datadir)//"/spsc_albedo.nc",NF90_NOWRITE,nid) 228 227 IF (ierr.NE.NF90_NOERR) THEN 229 228 write(*,*)'Problem opening spsc_albedo.nc (phymars/albedocaps.F90)' 230 write(*,*)'It should be in :',trim(data file),'/'229 write(*,*)'It should be in :',trim(datadir),'/' 231 230 write(*,*)'1) You can change this directory address in callfis.def with' 232 231 write(*,*)' datadir=/path/to/datafiles' … … 236 235 CALL ABORT 237 236 ELSE 238 write(*,*) "albedocaps: using file ",trim(data file)//"/spsc_albedo.nc"237 write(*,*) "albedocaps: using file ",trim(datadir)//"/spsc_albedo.nc" 239 238 ENDIF 240 239
Note: See TracChangeset
for help on using the changeset viewer.