Changeset 1918 for trunk/LMDZ.MARS/libf/aeronomars
- Timestamp:
- Apr 12, 2018, 8:35:48 AM (7 years ago)
- Location:
- trunk/LMDZ.MARS/libf/aeronomars
- Files:
-
- 4 edited
-
inichim_newstart.F90 (modified) (5 diffs)
-
param_read.F (modified) (5 diffs)
-
param_read_e107.F (modified) (10 diffs)
-
read_phototable.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/aeronomars/inichim_newstart.F90
r1660 r1918 5 5 USE vertical_layers_mod, ONLY: aps,bps 6 6 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, nbp_lev 7 USE datafile_mod, ONLY: datadir 7 8 implicit none 8 9 … … 36 37 37 38 include "callkeys.h" 38 include "datafile.h"39 39 40 40 ! inputs : … … 454 454 ! 2.1 open initialization files 455 455 456 open(210, iostat=ierr,file=trim(data file)//'/atmosfera_LMD_may.dat')456 open(210, iostat=ierr,file=trim(datadir)//'/atmosfera_LMD_may.dat') 457 457 if (ierr /= 0) then 458 458 write(*,*)'Error : cannot open file atmosfera_LMD_may.dat ' 459 459 write(*,*)'(in aeronomars/inichim_newstart.F)' 460 write(*,*)'It should be in :', trim(data file),'/'460 write(*,*)'It should be in :', trim(datadir),'/' 461 461 write(*,*)'1) You can change this path in callphys.def with' 462 462 write(*,*)' datadir=/path/to/datafiles/' … … 466 466 stop 467 467 end if 468 open(220, iostat=ierr,file=trim(data file)//'/atmosfera_LMD_min.dat')468 open(220, iostat=ierr,file=trim(datadir)//'/atmosfera_LMD_min.dat') 469 469 if (ierr /= 0) then 470 470 write(*,*)'Error : cannot open file atmosfera_LMD_min.dat ' 471 471 write(*,*)'(in aeronomars/inichim_newstart.F)' 472 write(*,*)'It should be in :', trim(data file),'/'472 write(*,*)'It should be in :', trim(datadir),'/' 473 473 write(*,*)'1) You can change this path in callphys.def with' 474 474 write(*,*)' datadir=/path/to/datafiles/' … … 479 479 end if 480 480 if(flagnitro) then 481 open(230, iostat=ierr,file=trim(data file)//'/atmosfera_LMD_nitr.dat')481 open(230, iostat=ierr,file=trim(datadir)//'/atmosfera_LMD_nitr.dat') 482 482 if (ierr.ne.0) then 483 483 write(*,*)'Error : cannot open file atmosfera_LMD_nitr.dat ' 484 484 write(*,*)'(in aeronomars/inichim_newstart.F)' 485 write(*,*)'It should be in :', datafile486 write(*,*)'1) You can change this directory address in'487 write(*,*)' file phymars/datafile.h'485 write(*,*)'It should be in :', trim(datadir),'/' 486 write(*,*)'1) You can change this path in callphys.def with' 487 write(*,*)' datadir=/path/to/datafiles/' 488 488 write(*,*)'2) If necessary atmosfera_LMD_nitr.dat (and others)' 489 489 write(*,*)' can be obtained online on:' -
trunk/LMDZ.MARS/libf/aeronomars/param_read.F
r1381 r1918 14 14 . fluxtop,ct1,ct2,p1,p2 15 15 16 use datafile_mod, only: datadir 17 16 18 implicit none 17 19 18 19 c common variables and constants20 include 'datafile.h'21 22 20 23 21 c local variables … … 34 32 !Tabulated column amount 35 33 open(210, status = 'old', 36 c $file=trim(data file)//'/EUVDAT/coln.dat',iostat=ierr)37 $file=trim(data file)//'/EUVDAT/param_v5/coln.dat',iostat=ierr)34 c $file=trim(datadir)//'/EUVDAT/coln.dat',iostat=ierr) 35 $file=trim(datadir)//'/EUVDAT/param_v5/coln.dat',iostat=ierr) 38 36 39 37 IF (ierr.NE.0) THEN 40 38 write(*,*)'cant find directory EUVDAT containing param_v5 subdir' 41 39 write(*,*)'(in aeronomars/param_read.F)' 42 write(*,*)'It should be in :', trim(data file),'/'40 write(*,*)'It should be in :', trim(datadir),'/' 43 41 write(*,*)'1) You can change this directory address in ' 44 42 write(*,*)' callphys.def with datadir=/path/to/dir' … … 50 48 51 49 !Tabulated photoabsorption coefficients 52 open(220,file=trim(data file)//'/EUVDAT/param_v5/j2_an.dat')53 open(230,file=trim(data file)//'/EUVDAT/param_v5/j3_an.dat')54 open(240,file=trim(data file)//'/EUVDAT/param_v5/j1_an.dat')55 open(250,file=trim(data file)//'/EUVDAT/param_v5/j2_bn.dat')56 open(260,file=trim(data file)//'/EUVDAT/param_v5/j2_cn.dat')57 open(300,file=trim(data file)//'/EUVDAT//param_v5/j2_dn.dat')58 open(270,file=trim(data file)//'/EUVDAT//param_v5/j1_bn.dat')59 open(280,file=trim(data file)//'/EUVDAT//param_v5/j1_cn.dat')60 open(290,file=trim(data file)//'/EUVDAT//param_v5/j1_dn.dat')61 open(150,file=trim(data file)//'/EUVDAT//param_v5/j4n.dat')62 open(160,file=trim(data file)//'/EUVDAT//param_v5/j5n.dat')63 open(170,file=trim(data file)//'/EUVDAT//param_v5/j6n.dat')64 open(180,file=trim(data file)//'/EUVDAT//param_v5/j7n.dat')65 open(390,file=trim(data file)//'/EUVDAT//param_v5/j8_an.dat')66 open(400,file=trim(data file)//'/EUVDAT//param_v5/j8_bn.dat')67 open(410,file=trim(data file)//'/EUVDAT//param_v5/j9n.dat')68 open(420,file=trim(data file)//'/EUVDAT//param_v5/j10_an.dat')69 open(430,file=trim(data file)//'/EUVDAT//param_v5/j10_bn.dat')70 open(440,file=trim(data file)//'/EUVDAT//param_v5/j10_cn.dat')71 open(450,file=trim(data file)//'/EUVDAT//param_v5/j11_an.dat')72 open(460,file=trim(data file)//'/EUVDAT//param_v5/j11_bn.dat')73 open(470,file=trim(data file)//'/EUVDAT//param_v5/j11_cn.dat')74 open(480,file=trim(data file)//'/EUVDAT//param_v5/j12n.dat')75 open(490,file=trim(data file)//'/EUVDAT//param_v5/j13_an.dat')76 open(500,file=trim(data file)//'/EUVDAT//param_v5/j13_bn.dat')77 open(510,file=trim(data file)//'/EUVDAT//param_v5/j13_cn.dat')50 open(220,file=trim(datadir)//'/EUVDAT/param_v5/j2_an.dat') 51 open(230,file=trim(datadir)//'/EUVDAT/param_v5/j3_an.dat') 52 open(240,file=trim(datadir)//'/EUVDAT/param_v5/j1_an.dat') 53 open(250,file=trim(datadir)//'/EUVDAT/param_v5/j2_bn.dat') 54 open(260,file=trim(datadir)//'/EUVDAT/param_v5/j2_cn.dat') 55 open(300,file=trim(datadir)//'/EUVDAT//param_v5/j2_dn.dat') 56 open(270,file=trim(datadir)//'/EUVDAT//param_v5/j1_bn.dat') 57 open(280,file=trim(datadir)//'/EUVDAT//param_v5/j1_cn.dat') 58 open(290,file=trim(datadir)//'/EUVDAT//param_v5/j1_dn.dat') 59 open(150,file=trim(datadir)//'/EUVDAT//param_v5/j4n.dat') 60 open(160,file=trim(datadir)//'/EUVDAT//param_v5/j5n.dat') 61 open(170,file=trim(datadir)//'/EUVDAT//param_v5/j6n.dat') 62 open(180,file=trim(datadir)//'/EUVDAT//param_v5/j7n.dat') 63 open(390,file=trim(datadir)//'/EUVDAT//param_v5/j8_an.dat') 64 open(400,file=trim(datadir)//'/EUVDAT//param_v5/j8_bn.dat') 65 open(410,file=trim(datadir)//'/EUVDAT//param_v5/j9n.dat') 66 open(420,file=trim(datadir)//'/EUVDAT//param_v5/j10_an.dat') 67 open(430,file=trim(datadir)//'/EUVDAT//param_v5/j10_bn.dat') 68 open(440,file=trim(datadir)//'/EUVDAT//param_v5/j10_cn.dat') 69 open(450,file=trim(datadir)//'/EUVDAT//param_v5/j11_an.dat') 70 open(460,file=trim(datadir)//'/EUVDAT//param_v5/j11_bn.dat') 71 open(470,file=trim(datadir)//'/EUVDAT//param_v5/j11_cn.dat') 72 open(480,file=trim(datadir)//'/EUVDAT//param_v5/j12n.dat') 73 open(490,file=trim(datadir)//'/EUVDAT//param_v5/j13_an.dat') 74 open(500,file=trim(datadir)//'/EUVDAT//param_v5/j13_bn.dat') 75 open(510,file=trim(datadir)//'/EUVDAT//param_v5/j13_cn.dat') 78 76 79 77 … … 228 226 229 227 !Parameters for the variation of the solar flux with 11 years cycle 230 open(100,file=trim(data file)//'/EUVDAT/param_v5/varflujo.dat')228 open(100,file=trim(datadir)//'/EUVDAT/param_v5/varflujo.dat') 231 229 read(100,*) 232 230 do i=1,24 … … 266 264 c CO2, O2, NO 267 265 268 open(120,file=trim(data file)//'/EUVDAT/param_v5/efdis_inter.dat')266 open(120,file=trim(datadir)//'/EUVDAT/param_v5/efdis_inter.dat') 269 267 read(120,*) 270 268 ! do i=1,21 -
trunk/LMDZ.MARS/libf/aeronomars/param_read_e107.F
r1888 r1918 16 16 . fluxtop,ct1,ct2,p1,p2 17 17 18 use datafile_mod, only: datadir 19 18 20 implicit none 19 21 20 22 21 23 c common variables and constants 22 include 'datafile.h'23 24 include "callkeys.h" 24 25 … … 39 40 !Tabulated column amount 40 41 open(210, status = 'old', 41 c $file=trim(data file)//'/EUVDAT/coln.dat',iostat=ierr)42 $file=trim(data file)//'/EUVDAT/param_v6/coln.dat',iostat=ierr)42 c $file=trim(datadir)//'/EUVDAT/coln.dat',iostat=ierr) 43 $file=trim(datadir)//'/EUVDAT/param_v6/coln.dat',iostat=ierr) 43 44 44 45 IF (ierr.NE.0) THEN 45 46 write(*,*)'cant find directory EUVDAT containing param_v6 subdir' 46 47 write(*,*)'(in aeronomars/param_read.F)' 47 write(*,*)'It should be in :', trim(data file),'/'48 write(*,*)'It should be in :', trim(datadir),'/' 48 49 write(*,*)'1) You can change this directory address in ' 49 50 write(*,*)' callphys.def with datadir=/path/to/dir' … … 55 56 56 57 !Tabulated photoabsorption coefficients 57 open(220,file=trim(data file)//'/EUVDAT/param_v6/trans2_an.dat')58 open(230,file=trim(data file)//'/EUVDAT/param_v6/trans3_an.dat')59 open(240,file=trim(data file)//'/EUVDAT/param_v6/trans1_an.dat')60 open(250,file=trim(data file)//'/EUVDAT/param_v6/trans2_bn.dat')61 open(260,file=trim(data file)//'/EUVDAT/param_v6/trans2_cn.dat')62 open(300,file=trim(data file)//'/EUVDAT/param_v6/trans2_dn.dat')63 open(270,file=trim(data file)//'/EUVDAT/param_v6/trans1_bn.dat')64 open(280,file=trim(data file)//'/EUVDAT/param_v6/trans1_cn.dat')65 open(290,file=trim(data file)//'/EUVDAT/param_v6/trans1_dn.dat')66 open(150,file=trim(data file)//'/EUVDAT/param_v6/trans4n.dat')67 open(160,file=trim(data file)//'/EUVDAT/param_v6/trans5n.dat')68 open(170,file=trim(data file)//'/EUVDAT/param_v6/trans6n.dat')69 open(180,file=trim(data file)//'/EUVDAT/param_v6/trans7n.dat')70 open(390,file=trim(data file)//'/EUVDAT/param_v6/trans8_an.dat')71 open(400,file=trim(data file)//'/EUVDAT/param_v6/trans8_bn.dat')72 open(410,file=trim(data file)//'/EUVDAT/param_v6/trans9n.dat')73 open(420,file=trim(data file)//'/EUVDAT/param_v6/trans10_an.dat')74 open(430,file=trim(data file)//'/EUVDAT/param_v6/trans10_bn.dat')75 open(440,file=trim(data file)//'/EUVDAT/param_v6/trans10_cn.dat')76 open(450,file=trim(data file)//'/EUVDAT/param_v6/trans11_an.dat')77 open(460,file=trim(data file)//'/EUVDAT/param_v6/trans11_bn.dat')78 open(470,file=trim(data file)//'/EUVDAT/param_v6/trans11_cn.dat')79 open(480,file=trim(data file)//'/EUVDAT/param_v6/trans12n.dat')80 open(490,file=trim(data file)//'/EUVDAT/param_v6/trans13_an.dat')81 open(500,file=trim(data file)//'/EUVDAT/param_v6/trans13_bn.dat')82 open(510,file=trim(data file)//'/EUVDAT/param_v6/trans13_cn.dat')58 open(220,file=trim(datadir)//'/EUVDAT/param_v6/trans2_an.dat') 59 open(230,file=trim(datadir)//'/EUVDAT/param_v6/trans3_an.dat') 60 open(240,file=trim(datadir)//'/EUVDAT/param_v6/trans1_an.dat') 61 open(250,file=trim(datadir)//'/EUVDAT/param_v6/trans2_bn.dat') 62 open(260,file=trim(datadir)//'/EUVDAT/param_v6/trans2_cn.dat') 63 open(300,file=trim(datadir)//'/EUVDAT/param_v6/trans2_dn.dat') 64 open(270,file=trim(datadir)//'/EUVDAT/param_v6/trans1_bn.dat') 65 open(280,file=trim(datadir)//'/EUVDAT/param_v6/trans1_cn.dat') 66 open(290,file=trim(datadir)//'/EUVDAT/param_v6/trans1_dn.dat') 67 open(150,file=trim(datadir)//'/EUVDAT/param_v6/trans4n.dat') 68 open(160,file=trim(datadir)//'/EUVDAT/param_v6/trans5n.dat') 69 open(170,file=trim(datadir)//'/EUVDAT/param_v6/trans6n.dat') 70 open(180,file=trim(datadir)//'/EUVDAT/param_v6/trans7n.dat') 71 open(390,file=trim(datadir)//'/EUVDAT/param_v6/trans8_an.dat') 72 open(400,file=trim(datadir)//'/EUVDAT/param_v6/trans8_bn.dat') 73 open(410,file=trim(datadir)//'/EUVDAT/param_v6/trans9n.dat') 74 open(420,file=trim(datadir)//'/EUVDAT/param_v6/trans10_an.dat') 75 open(430,file=trim(datadir)//'/EUVDAT/param_v6/trans10_bn.dat') 76 open(440,file=trim(datadir)//'/EUVDAT/param_v6/trans10_cn.dat') 77 open(450,file=trim(datadir)//'/EUVDAT/param_v6/trans11_an.dat') 78 open(460,file=trim(datadir)//'/EUVDAT/param_v6/trans11_bn.dat') 79 open(470,file=trim(datadir)//'/EUVDAT/param_v6/trans11_cn.dat') 80 open(480,file=trim(datadir)//'/EUVDAT/param_v6/trans12n.dat') 81 open(490,file=trim(datadir)//'/EUVDAT/param_v6/trans13_an.dat') 82 open(500,file=trim(datadir)//'/EUVDAT/param_v6/trans13_bn.dat') 83 open(510,file=trim(datadir)//'/EUVDAT/param_v6/trans13_cn.dat') 83 84 84 85 … … 233 234 234 235 !Parameters for the variation of the solar flux with 11 years cycle 235 open(620,file=trim(data file)//'/EUVDAT/param_v6/fit_js_e107.dat')236 open(620,file=trim(datadir)//'/EUVDAT/param_v6/fit_js_e107.dat') 236 237 do i=1,ninter 237 238 read(620,*) … … 274 275 endif 275 276 276 open(640,file=trim(data file)//'/EUVDAT/param_v6/'//filename)277 open(640,file=trim(datadir)//'/EUVDAT/param_v6/'//filename) 277 278 read(640,*) 278 279 do i=1,669 … … 316 317 c CO2, O2, NO 317 318 318 ! open(120,file=trim(data file)//'/EUVDAT/param_v5/efdis_inter.dat')319 ! open(120,file=trim(datadir)//'/EUVDAT/param_v5/efdis_inter.dat') 319 320 ! read(120,*) 320 321 !! do i=1,21 … … 375 376 376 377 !4 possible channels for CO2 ionization 377 open(130,file=trim(data file)//'/EUVDAT'//378 open(130,file=trim(datadir)//'/EUVDAT'// 378 379 $ '/co2ion_branchingratio_schunkandnagy2000_param.dat') 379 380 do inter=1,16 … … 408 409 409 410 !2 possible channels for O2 ionization 410 open(131,file=trim(data file)//'/EUVDAT'//411 open(131,file=trim(datadir)//'/EUVDAT'// 411 412 $ '/o2ion_branchingratio_schunkandnagy2000_param.dat') 412 413 do inter=1,23 … … 430 431 431 432 !2 channels for N2 ionization 432 open(132,file=trim(data file)//'/EUVDAT'//433 open(132,file=trim(datadir)//'/EUVDAT'// 433 434 $ '/n2ion_branchingratio_schunkandnagy2000_param.dat') 434 435 do inter=1,15 … … 453 454 454 455 !3 channels for CO ionization 455 open(133,file=trim(data file)//'/EUVDAT'//456 open(133,file=trim(datadir)//'/EUVDAT'// 456 457 $ '/coion_branchingratio_schunkandnagy2000_param.dat') 457 458 do inter=1,16 -
trunk/LMDZ.MARS/libf/aeronomars/read_phototable.F90
r1430 r1918 21 21 !*********************************************************************** 22 22 23 use ioipsl_getincom 23 use ioipsl_getincom, only: getin 24 use datafile_mod, only: datadir 24 25 25 26 implicit none 26 27 27 #include "chimiedata.h" 28 #include "datafile.h" 28 include "chimiedata.h" 29 29 30 30 !ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc … … 50 50 51 51 open(fic, form = 'formatted', status = 'old', & 52 file =trim(data file)//"/"//trim(phototable),iostat=ierr)52 file =trim(datadir)//"/"//trim(phototable),iostat=ierr) 53 53 54 54 if (ierr /= 0) THEN 55 55 write(*,*)'Error : cannot open photolysis lookup table ', trim(phototable) 56 write(*,*)'It should be in :',trim(data file),'/'56 write(*,*)'It should be in :',trim(datadir),'/' 57 57 write(*,*)'1) You can change this directory in callphys.def' 58 58 write(*,*)' with:'
Note: See TracChangeset
for help on using the changeset viewer.
