Changeset 1720 for trunk/LMDZ.MARS/libf/phymars/nucleaCO2.F

Timestamp:

Jul 6, 2017, 5:12:58 PM (7 years ago)

Author:

jaudouard

Message:

Update on CO2 ice clouds scheme

File:

• trunk/LMDZ.MARS/libf/phymars/nucleaCO2.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/nucleaCO2.F

@@ -65 +65 @@
       mtetalocal = dble(mtetaco2)  !! use mtetalocal for better performance
       mtetalocalh=dble(mteta)
-
-cccccccccccccccccccccccccccccccccccccccccccccccccc
-ccccccccccc ESSAIS TN MTETA = F (T) cccccccccccccc
-c      if (temp .gt. 200) then
-c         mtetalocal = mtetalocal
-c      else if (temp .lt. 190) then
-c         mtetalocal = mtetalocal-0.05
-c      else
-c         mtetalocal = mtetalocal - (190-temp)*0.005
-c      endif
-c----------------exp law, see Trainer 2008, J. Phys. Chem. C 2009, 113, 2036\u20132040
-       !mtetalocal = max(mtetalocal - 6005*exp(-0.065*temp),0.1)
-       !mtetalocal = max(mtetalocal - 6005*exp(-0.068*temp),0.1)
-               !print*, mtetalocal, temp
-cccccccccccccccccccccccccccccccccccccccccccccccccc
-cccccccccccccccccccccccccccccccccccccccccccccccccc 
-c      IF (firstcall) THEN
-c          print*, ' '  
-c          print*, 'dear user, please keep in mind that'
-c          print*, 'contact parameter IS constant'
-          !print*, 'contact parameter IS NOT constant:'
-          !print*, 'max(mteta - 6005*exp(-0.065*temp),0.1)'
-          !print*, 'max(mteta - 6005*exp(-0.068*temp),0.1)'
-c          print*, ' '  
-c         firstcall=.false.
-c     END IF
-cccccccccccccccccccccccccccccccccccccccccccccccccc
-cccccccccccccccccccccccccccccccccccccccccccccccccc
-   
-c      write(*,*) "IN nuc, SAT = ",sat
-c      write(*,*) "IN nuc, mtetalocal = ",mtetalocal