Changeset 2870 for trunk/LMDZ.GENERIC/libf/phystd/interpolateHeCH4.F90

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:

• trunk/LMDZ.GENERIC/libf/phystd/interpolateHeCH4.F90 (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/interpolateHeCH4.F90

@@ -55 +55 @@
 
   character*20 bleh
-  double precision blah, Ttemp
+  double precision blah, Ttemp, ztemp
   integer nres
   integer ind
 
-  if(temp.gt.350)then
+  ztemp=temp
+
+  if(ztemp.gt.350)then
     if (strictboundcia) then
       if (is_master) then
@@ -67 +69 @@
       stop
     else
-      if (is_master) then
-        print*,'Your temperatures are too high for this He-CH4 CIA dataset'
-        print*,'you have chosen strictboundcia = ', strictboundcia
-        print*,'*********************************************************'
-        print*,' we allow model to continue but with temp = 350          '
-        print*,'  ...       for He-CH4 CIA dataset         ...           '
-        print*,'  ... we assume we know what you are doing ...           '
-        print*,'*********************************************************'
-      endif
-      temp = 350
+      !if (is_master) then
+      !  print*,'Your temperatures are too high for this He-CH4 CIA dataset'
+      !  print*,'you have chosen strictboundcia = ', strictboundcia
+      !  print*,'*********************************************************'
+      !  print*,' we allow model to continue but with temp = 350          '
+      !  print*,'  ...       for He-CH4 CIA dataset         ...           '
+      !  print*,'  ... we assume we know what you are doing ...           '
+      !  print*,'*********************************************************'
+      !endif
+      ztemp = 350
     endif
-  elseif(temp.lt.40)then
+  elseif(ztemp.lt.40)then
     if (strictboundcia) then
       if (is_master) then
@@ -86 +88 @@
       stop
     else
-      if (is_master) then
-        print*,'Your temperatures are too low for this He-CH4 CIA dataset'
-        print*,'you have chosen strictboundcia = ', strictboundcia
-        print*,'*********************************************************'
-        print*,' we allow model to continue but with temp = 40           '
-        print*,'  ...       for He-CH4 CIA dataset         ...           '
-        print*,'  ... we assume we know what you are doing ...           '
-        print*,'*********************************************************'
-      endif
-      temp = 40
+      !if (is_master) then
+      !  print*,'Your temperatures are too low for this He-CH4 CIA dataset'
+      !  print*,'you have chosen strictboundcia = ', strictboundcia
+      !  print*,'*********************************************************'
+      !  print*,' we allow model to continue but with temp = 40           '
+      !  print*,'  ...       for He-CH4 CIA dataset         ...           '
+      !  print*,'  ... we assume we know what you are doing ...           '
+      !  print*,'*********************************************************'
+      !endif
+      ztemp = 40
     endif
   endif
@@ -154 +156 @@
   endif
 
-  call bilinearbig(nS,nT,wn_arr,temp_arr,abs_arr,wn,temp,abcoef,ind)
+  call bilinearbig(nS,nT,wn_arr,temp_arr,abs_arr,wn,ztemp,abcoef,ind)
 
 !     print*,'the absorption is ',abcoef,' cm^5 molecule^-2'