Ignore:
Timestamp:
Apr 16, 2012, 12:22:40 PM (13 years ago)
Author:
jleconte
Message:
  • Added consistency checks for calculations including water and global1d+diurnal.
  • Corrected small bugs in precipitation scheme
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/libf/phystd/inifis.F

    r596 r622  
    438438         endif
    439439
    440 ! Test of incompatibility:
    441 
    442440         write(*,*) "Gravitationnal sedimentation ?"
    443441         sedimentation=.true. ! default value
    444442         call getin("sedimentation",sedimentation)
    445443         write(*,*) " sedimentation = ",sedimentation
    446 
    447 
    448 ! Test of incompatibility:
    449444
    450445         write(*,*) "Compute water cycle ?"
     
    453448         write(*,*) " water = ",water
    454449         
     450! Test of incompatibility:
     451! if water is true, there should be at least a tracer
     452         if (water.and.(.not.tracer)) then
     453           print*,'if water is ON, tracer must be ON too!'
     454           stop
     455         endif
     456
    455457         write(*,*) "Include water condensation ?"
    456458         watercond=.false. ! default value
     
    458460         write(*,*) " watercond = ",watercond
    459461
     462! Test of incompatibility:
     463! if watercond is used, then water should be used too
     464         if (watercond.and.(.not.water)) then
     465           print*,'if watercond is used, water should be used too'
     466           stop
     467         endif
     468
    460469         write(*,*) "Include water precipitation ?"
    461470         waterrain=.false. ! default value
     
    497506         call getin("Tsaldiff",Tsaldiff)
    498507         write(*,*) " Tsaldiff = ",Tsaldiff
    499 
    500 ! Test of incompatibility:
    501 ! if watercond is used, then water should be used too
    502 
    503          if (watercond.and.(.not.watercond)) then
    504            print*,'if watercond is used, water should be used too'
    505            stop
    506          endif
    507508
    508509         write(*,*) "Does user want to force cpp and mugaz?"
Note: See TracChangeset for help on using the changeset viewer.