Ignore:
Timestamp:
Jan 16, 2023, 4:47:10 PM (23 months ago)
Author:
jleconte
Message:

Removed too many prints in CIA interpolation + correction for hot H2H2

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.GENERIC/libf/phystd/interpolateH2H2.F90

    r2837 r2870  
    5858
    5959      character*20 bleh
    60       double precision blah, Ttemp
     60      double precision blah, Ttemp, ztemp
    6161      integer nres
    6262
    6363      integer ind
     64
     65      ztemp=temp
    6466      if ((H2orthopara_mixture .eq. "hot")) then
    65         if (temp .gt. 7000.) then
     67        if (ztemp .gt. 7000.) then
    6668          if (strictboundcia) then
    6769            if (is_master) then
    … …  
    7072            stop
    7173          else
    72             if (is_master) then
    73               print*,'Your temperatures are too high for this H2-H2 CIA dataset (Hot Jupiter case)'
    74               print*,'you have chosen strictboundcia = ', strictboundcia
    75               print*,'*********************************************************'
    76               print*,' we allow model to continue but with temp = 7000          '
    77               print*,'  ...       for H2-H2 CIA dataset          ...           '
    78               print*,'  ... we assume we know what you are doing ...           '
    79               print*,'*********************************************************'
    80             endif !is_master
    81             temp = 7000.
     74            !if (is_master) then
     75            !  print*,'Your temperatures are too high for this H2-H2 CIA dataset (Hot Jupiter case)'
     76            !  print*,'you have chosen strictboundcia = ', strictboundcia
     77            !  print*,'*********************************************************'
     78            !  print*,' we allow model to continue but with temp = 7000          '
     79            !  print*,'  ...       for H2-H2 CIA dataset          ...           '
     80            !  print*,'  ... we assume we know what you are doing ...           '
     81            !  print*,'*********************************************************'
     82            !endif !is_master
     83            ztemp = 7000.
    8284          endif !strictboundcia
    8385        endif !(temp .gt. 7000.)
    8486      else ! if not "hot"
    85         if(temp.gt.400)then
     87        if(ztemp.gt.400)then
    8688          if (strictboundcia) then
    8789            if (is_master) then
    … …  
    9294            stop
    9395          else
    94             if (is_master) then
    95               print*,'Your temperatures are too high for this H2-H2 CIA dataset'
    96               print*,'you have chosen strictboundcia = ', strictboundcia
    97               print*,'*********************************************************'
    98               print*,' we allow model to continue but with temp = 400          '
    99               print*,'  ...       for H2-H2 CIA dataset          ...           '
    100               print*,'  ... we assume we know what you are doing ...           '
    101               print*,'*********************************************************'
    102             endif !is_master
    103             temp = 400
     96            !if (is_master) then
     97            !  print*,'Your temperatures are too high for this H2-H2 CIA dataset'
     98            !  print*,'you have chosen strictboundcia = ', strictboundcia
     99            !  print*,'*********************************************************'
     100            !  print*,' we allow model to continue but with temp = 400          '
     101            !  print*,'  ...       for H2-H2 CIA dataset          ...           '
     102            !  print*,'  ... we assume we know what you are doing ...           '
     103            !  print*,'*********************************************************'
     104            !endif !is_master
     105            ztemp = 400
    104106          endif !of strictboundcia
    105         elseif(temp.lt.40)then     
     107        elseif(ztemp.lt.40)then     
    106108          if (strictboundcia) then
    107109            if (is_master) then
    … …  
    112114            stop
    113115          else
    114             if (is_master) then
    115               print*,'Your temperatures are too low for this H2-H2 CIA dataset'
    116               print*,'you have chosen strictboundcia = ', strictboundcia
    117               print*,'*********************************************************'
    118               print*,' we allow model to continue but with temp = 40           '
    119               print*,'  ...       for H2-H2 CIA dataset          ...           '
    120               print*,'  ... we assume we know what you are doing ...           '
    121               print*,'*********************************************************'
    122             endif !is_master
    123             temp = 40
     116            !if (is_master) then
     117            !  print*,'Your temperatures are too low for this H2-H2 CIA dataset'
     118            !  print*,'you have chosen strictboundcia = ', strictboundcia
     119            !  print*,'*********************************************************'
     120            !  print*,' we allow model to continue but with temp = 40           '
     121            !  print*,'  ...       for H2-H2 CIA dataset          ...           '
     122            !  print*,'  ... we assume we know what you are doing ...           '
     123            !  print*,'*********************************************************'
     124            !endif !is_master
     125            ztemp = 40
    124126          endif !of strictboundcia       
    125127        endif ! of (temp .gt. 400)
    … …  
    216218           
    217219           
    218          call bilinearbig(nS,nT,wn_arr,temp_arr,abs_arr,wn,temp,abcoef,ind)
     220         call bilinearbig(nS,nT,wn_arr,temp_arr,abs_arr,wn,ztemp,abcoef,ind)
    219221
    220222         !print*,'the absorption is ',abcoef,' cm^5 molecule^-2'
Note: See TracChangeset for help on using the changeset viewer.