Ignore:
Timestamp:
Apr 21, 2022, 5:34:49 PM (3 years ago)
Author:
aslmd
Message:

CIA : modified the loading of file to use high temperature dataset

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r2662 r2667  
    6262
    6363      integer ind
    64      
    65       if(temp.gt.400)then
    66         if (strictboundcia) then
    67           if (is_master) then
    68             print*,'Your temperatures are too high for this H2-H2 CIA dataset. If you '
    69             print*,'really want to run simulations with hydrogen at T > 400 K, contact'
    70             print*,'Robin Wordsworth [rwordsworth@uchicago.edu].'
    71           endif
    72           stop
    73         else
    74           if (is_master) then
    75             print*,'Your temperatures are too high for this H2-H2 CIA dataset'
    76             print*,'you have chosen strictboundcia = ', strictboundcia
    77             print*,'*********************************************************'
    78             print*,' we allow model to continue but with temp = 400          '
    79             print*,'  ...       for H2-H2 CIA dataset          ...           '
    80             print*,'  ... we assume we know what you are doing ...           '
    81             print*,'*********************************************************'
    82           endif
    83           temp = 400
    84         endif
    85       elseif(temp.lt.40)then     
    86         if (strictboundcia) then
    87           if (is_master) then
    88             print*,'Your temperatures are too low for this H2-H2 CIA dataset. If you '
    89             print*,'really want to run simulations with hydrogen at T < 40 K, contact'
    90             print*,'Robin Wordsworth [rwordsworth@uchicago.edu].'
    91           endif
    92           stop
    93         else
    94           if (is_master) then
    95             print*,'Your temperatures are too low 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 = 40           '
    99             print*,'  ...       for H2-H2 CIA dataset          ...           '
    100             print*,'  ... we assume we know what you are doing ...           '
    101             print*,'*********************************************************'
    102           endif
    103           temp = 40
    104         endif           
    105       endif
    106 
     64      if ((H2orthopara_mixture .eq. "hot")) then
     65        ! print*,"We're in the Hot Jupiter case"
     66        if (temp .gt. 3000.) then
     67          if (strictboundcia) then
     68            if (is_master) then
     69              print*,'Your temperatures are too high for this H2-H2 CIA dataset (Hot Jupiter case). If you '
     70              print*,'really want to run simulations with hydrogen at T > 400 K, contact'
     71              print*,'Robin Wordsworth [rwordsworth@uchicago.edu].'
     72            endif !is_master
     73            stop
     74          else
     75            if (is_master) then
     76              print*,'Your temperatures are too high for this H2-H2 CIA dataset (Hot Jupiter case)'
     77              print*,'you have chosen strictboundcia = ', strictboundcia
     78              print*,'*********************************************************'
     79              print*,' we allow model to continue but with temp = 3000          '
     80              print*,'  ...       for H2-H2 CIA dataset          ...           '
     81              print*,'  ... we assume we know what you are doing ...           '
     82              print*,'*********************************************************'
     83            endif !is_master
     84            temp = 3000.
     85          endif !strictboundcia
     86        endif !(temp .gt. 3000.)
     87      else ! if not Hot Jupiter
     88        if(temp.gt.400)then
     89          if (strictboundcia) then
     90            if (is_master) then
     91              print*,'Your temperatures are too high for this H2-H2 CIA dataset. If you '
     92              print*,'really want to run simulations with hydrogen at T > 400 K, contact'
     93              print*,'Robin Wordsworth [rwordsworth@uchicago.edu].'
     94            endif !is_master
     95            stop
     96          else
     97            if (is_master) then
     98              print*,'Your temperatures are too high for this H2-H2 CIA dataset'
     99              print*,'you have chosen strictboundcia = ', strictboundcia
     100              print*,'*********************************************************'
     101              print*,' we allow model to continue but with temp = 400          '
     102              print*,'  ...       for H2-H2 CIA dataset          ...           '
     103              print*,'  ... we assume we know what you are doing ...           '
     104              print*,'*********************************************************'
     105            endif !is_master
     106            temp = 400
     107          endif !of strictboundcia
     108        elseif(temp.lt.40)then     
     109          if (strictboundcia) then
     110            if (is_master) then
     111              print*,'Your temperatures are too low for this H2-H2 CIA dataset. If you '
     112              print*,'really want to run simulations with hydrogen at T < 40 K, contact'
     113              print*,'Robin Wordsworth [rwordsworth@uchicago.edu].'
     114            endif ! is_master
     115            stop
     116          else
     117            if (is_master) then
     118              print*,'Your temperatures are too low for this H2-H2 CIA dataset'
     119              print*,'you have chosen strictboundcia = ', strictboundcia
     120              print*,'*********************************************************'
     121              print*,' we allow model to continue but with temp = 40           '
     122              print*,'  ...       for H2-H2 CIA dataset          ...           '
     123              print*,'  ... we assume we know what you are doing ...           '
     124              print*,'*********************************************************'
     125            endif !is_master
     126            temp = 40
     127          endif !of strictboundcia       
     128        endif ! of (temp .gt. 400)
     129      endif ! of ((H2orthopara_mixture .eq. "hot").and. (temp .gt. 3000.))
    107130      amagat = (273.15/temp)*(pres/101325.0)
    108131
    … …  
    119142           else if (H2orthopara_mixture.eq."equilibrium") then
    120143             dt_file=TRIM(datadir)//'/continuum_data/H2-H2_eq_2011.cia'
     144           else if (H2orthopara_mixture.eq."hot") then
     145            dt_file=TRIM(datadir)//'/continuum_data/H2-H2_2011.cia'
     146            ns = 9981
    121147           endif
    122148         else if (versH2H2cia.eq.2018) then
Note: See TracChangeset for help on using the changeset viewer.