Changeset 3659

Timestamp:

Feb 27, 2025, 2:53:09 PM (4 months ago)

Author:

emoisan

Message:

Titan CRM:
Deal with cloud fraction Fcloudy=1.
EMo

File:

• trunk/LMDZ.TITAN/libf/phytitan/callcorrk.F90 (modified) (4 diffs)

trunk/LMDZ.TITAN/libf/phytitan/callcorrk.F90

@@ -438 +438 @@
          if (tlevrad(k).lt.tplanckmin) then
             print*,'Minimum temperature is outside the boundaries for'
-            print*,'Planck function integration set in callphys.def, aborting.'
+            print*,'Planck function integration set in callphys.def,'
+            print*,' aborting.'
             print*,"k=",k," tlevrad(k)=",tlevrad(k)
             print*,"tplanckmin=",tplanckmin
@@ -537 +538 @@
          wbarv = (1-Fcloudy) * wbarv_cc + Fcloudy * wbarv_dc
          cosbv = (1-Fcloudy) * cosbv_cc + Fcloudy * cosbv_dc
-         
+        
          ! Diagnostics for clouds :
          if (callclouds) then
@@ -546 +547 @@
                diag_opttv_dc(:,:,1) = 1.d-30
             endwhere
-            diag_opttv(:,:,1) = -log( (1-Fcloudy)*exp(-diag_opttv_cc(:,:,1)) + Fcloudy*exp(-diag_opttv_dc(:,:,1)) )
-            diag_opttv(:,:,2) = (1-Fcloudy) * diag_opttv_cc(:,:,2) + Fcloudy * diag_opttv_dc(:,:,2)
-            diag_opttv(:,:,3) = (1-Fcloudy) * diag_opttv_cc(:,:,3) + Fcloudy * diag_opttv_dc(:,:,3)
-            diag_opttv(:,:,4) = -log( (1-Fcloudy)*exp(-diag_opttv_cc(:,:,4)) + Fcloudy*exp(-diag_opttv_dc(:,:,4)) )
-            diag_opttv(:,:,5) = -log( (1-Fcloudy)*exp(-diag_opttv_cc(:,:,5)) + Fcloudy*exp(-diag_opttv_dc(:,:,5)) )
-            diag_opttv(:,:,6) = -log( (1-Fcloudy)*exp(-diag_opttv_cc(:,:,6)) + Fcloudy*exp(-diag_opttv_dc(:,:,6)) )
+            if (Fcloudy .eq. 1) then !avoid log(0) errors
+                diag_opttv(:,:,1) = diag_opttv_dc(:,:,1)
+                diag_opttv(:,:,2) = diag_opttv_dc(:,:,2)
+                diag_opttv(:,:,3) = diag_opttv_dc(:,:,3)
+                diag_opttv(:,:,4) = diag_opttv_dc(:,:,4)
+                diag_opttv(:,:,5) = diag_opttv_dc(:,:,5)
+                diag_opttv(:,:,6) = diag_opttv_dc(:,:,6)
+            else !correct for cloud fraction
+                diag_opttv(:,:,1) = -log( (1-Fcloudy)*exp(-diag_opttv_cc(:,:,1)) + Fcloudy*exp(-diag_opttv_dc(:,:,1)) )
+                diag_opttv(:,:,2) = (1-Fcloudy) * diag_opttv_cc(:,:,2) + Fcloudy * diag_opttv_dc(:,:,2)
+                diag_opttv(:,:,3) = (1-Fcloudy) * diag_opttv_cc(:,:,3) + Fcloudy * diag_opttv_dc(:,:,3)
+                diag_opttv(:,:,4) = -log( (1-Fcloudy)*exp(-diag_opttv_cc(:,:,4)) + Fcloudy*exp(-diag_opttv_dc(:,:,4)) )
+                diag_opttv(:,:,5) = -log( (1-Fcloudy)*exp(-diag_opttv_cc(:,:,5)) + Fcloudy*exp(-diag_opttv_dc(:,:,5)) )
+                diag_opttv(:,:,6) = -log( (1-Fcloudy)*exp(-diag_opttv_cc(:,:,6)) + Fcloudy*exp(-diag_opttv_dc(:,:,6)) )
+            endif
          endif
          
@@ -644 +654 @@
                diag_optti_dc(:,:,1) = 1.d-30
             endwhere
-            diag_optti(:,:,1) = -log( (1-Fcloudy)*exp(-diag_optti_cc(:,:,1)) + Fcloudy*exp(-diag_optti_dc(:,:,1)) )
-            diag_optti(:,:,2) = (1-Fcloudy) * diag_optti_cc(:,:,2) + Fcloudy * diag_optti_dc(:,:,2)
-            diag_optti(:,:,3) = (1-Fcloudy) * diag_optti_cc(:,:,3) + Fcloudy * diag_optti_dc(:,:,3)
-            diag_optti(:,:,4) = -log( (1-Fcloudy)*exp(-diag_optti_cc(:,:,4)) + Fcloudy*exp(-diag_optti_dc(:,:,4)) )
-            diag_optti(:,:,5) = -log( (1-Fcloudy)*exp(-diag_optti_cc(:,:,5)) + Fcloudy*exp(-diag_optti_dc(:,:,5)) )
-            diag_optti(:,:,6) = -log( (1-Fcloudy)*exp(-diag_optti_cc(:,:,6)) + Fcloudy*exp(-diag_optti_dc(:,:,6)) )
+            if (Fcloudy .eq. 1) then !avoid log(0) errors
+                diag_optti(:,:,1) = diag_optti_dc(:,:,1)
+                diag_optti(:,:,2) = diag_optti_dc(:,:,2)
+                diag_optti(:,:,3) = diag_optti_dc(:,:,3)
+                diag_optti(:,:,4) = diag_optti_dc(:,:,4)
+                diag_optti(:,:,5) = diag_optti_dc(:,:,5)
+                diag_optti(:,:,6) = diag_optti_dc(:,:,6)
+            else !correct for cloud fraction
+                diag_optti(:,:,1) = -log( (1-Fcloudy)*exp(-diag_optti_cc(:,:,1)) + Fcloudy*exp(-diag_optti_dc(:,:,1)) )
+                diag_optti(:,:,2) = (1-Fcloudy) * diag_optti_cc(:,:,2) + Fcloudy * diag_optti_dc(:,:,2)
+                diag_optti(:,:,3) = (1-Fcloudy) * diag_optti_cc(:,:,3) + Fcloudy * diag_optti_dc(:,:,3)
+                diag_optti(:,:,4) = -log( (1-Fcloudy)*exp(-diag_optti_cc(:,:,4)) + Fcloudy*exp(-diag_optti_dc(:,:,4)) )
+                diag_optti(:,:,5) = -log( (1-Fcloudy)*exp(-diag_optti_cc(:,:,5)) + Fcloudy*exp(-diag_optti_dc(:,:,5)) )
+                diag_optti(:,:,6) = -log( (1-Fcloudy)*exp(-diag_optti_cc(:,:,6)) + Fcloudy*exp(-diag_optti_dc(:,:,6)) )
+            endif
          endif