Ignore:
Timestamp:
Mar 8, 2017, 10:37:31 AM (8 years ago)
Author:
jvatant
Message:

Re-plug chemistry of old Titan
+ minor correction on reading Titan's startfiles and haze routine
JVO

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90

    r1670 r1672  
    1818                              init_time, daysec, dtphys
    1919  use comcstfi_mod, only: rad, cpp, g, r, rcp, &
    20                           mugaz, pi, avocado
     20                          mugaz, pi, avocado, kbol
    2121  use planete_mod, only: nres
    2222  use planetwide_mod, only: planetwide_sumval
     
    8484  pi=2.*asin(1.)
    8585  avocado = 6.02214179e23   ! added by RW
     86  kbol = 1.38064852e-23  ! added by JVO for Titan chem
    8687
    8788  ! Initialize some "temporal and calendar" related variables
     
    162163     flatten = 0.0
    163164     call getin_p("flatten", flatten)
    164      write(*,*) "flatten = ", flatten
    165          
     165     write(*,*) "flatten = ", flatten         
    166166
    167167     write(*,*) "Needed if oblate=.true.: J2"
     
    300300     call getin_p("graybody",graybody)
    301301     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
    302322
    303323! Soil model
     
    497517  PRINT*,'        or each ',iradia*dtphys,' seconds'
    498518  PRINT*
    499 
     519  PRINT*,'inifis: Chemistry is computed:'
     520  PRINT*,'           each ',ichim,' physical time-step'
     521  PRINT*,'        or each ',ichim*dtphys,' seconds'
     522  PRINT*
    500523
    501524!-----------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.