Changeset 1672 for trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90
- Timestamp:
- Mar 8, 2017, 10:37:31 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90
r1670 r1672 18 18 init_time, daysec, dtphys 19 19 use comcstfi_mod, only: rad, cpp, g, r, rcp, & 20 mugaz, pi, avocado 20 mugaz, pi, avocado, kbol 21 21 use planete_mod, only: nres 22 22 use planetwide_mod, only: planetwide_sumval … … 84 84 pi=2.*asin(1.) 85 85 avocado = 6.02214179e23 ! added by RW 86 kbol = 1.38064852e-23 ! added by JVO for Titan chem 86 87 87 88 ! Initialize some "temporal and calendar" related variables … … 162 163 flatten = 0.0 163 164 call getin_p("flatten", flatten) 164 write(*,*) "flatten = ", flatten 165 165 write(*,*) "flatten = ", flatten 166 166 167 167 write(*,*) "Needed if oblate=.true.: J2" … … 300 300 call getin_p("graybody",graybody) 301 301 write(*,*)" graybody = ",graybody 302 303 ! Chemistry 304 305 write(*,*) "Run with or without chemistry?" 306 callchim=.false. ! default value 307 call getin_p("callchim",callchim) 308 write(*,*) " callchim = ",callchim 309 310 ! sanity check 311 if (callchim.and.(.not.tracer)) then 312 print*,"You are running chemistry without tracer" 313 print*,"Please start again with tracer =.true." 314 stop 315 endif 316 317 write(*,*)"Chemistry is computed every ichim", & 318 " physical timestep" 319 ichim=1 ! default value 320 call getin_p("ichim",ichim) 321 write(*,*)" ichim = ",ichim 302 322 303 323 ! Soil model … … 497 517 PRINT*,' or each ',iradia*dtphys,' seconds' 498 518 PRINT* 499 519 PRINT*,'inifis: Chemistry is computed:' 520 PRINT*,' each ',ichim,' physical time-step' 521 PRINT*,' or each ',ichim*dtphys,' seconds' 522 PRINT* 500 523 501 524 !-----------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.