Changeset 4153 for trunk/LMDZ.TITAN/libf/phytitan/optcv.F90

Timestamp:

Mar 25, 2026, 2:01:37 PM (2 weeks ago)

Author:

aarfaux

Message:

Titan PCM:

Add a correction to the shortwave heating rates based on htrdr calculations, and
additional related modifications.

File:

• trunk/LMDZ.TITAN/libf/phytitan/optcv.F90 (modified) (3 diffs)

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

@@ -133 +133 @@
   logical,save :: firstcall=.true.
 !$OMP THREADPRIVATE(firstcall)
-
 
   !! AS: to save time in computing continuum (see bilinearbig)
@@ -537 +536 @@
    DO L=1,L_NLAYRAD-1
       K = 2*L+1
-           atemp(L,NW) = ASF_T(K,NW)*DHAZES_T(K,NW) + ASF_T(K+1,NW)*DHAZES_T(K+1,NW)
+      atemp(L,NW) = ASF_T(K,NW)*DHAZES_T(K,NW) + ASF_T(K+1,NW)*DHAZES_T(K+1,NW)
       btemp(L,NW) = DHAZES_T(K,NW) + DHAZES_T(K+1,NW)
-           ctemp(L,NW) = btemp(L,NW) + 0.9999*(TRAY(K,NW) + TRAY(K+1,NW)) ! JVO 2017 : does this 0.999 is really meaningful ?
-           btemp(L,NW) = btemp(L,NW) + TRAY(K,NW) + TRAY(K+1,NW)
+      ctemp(L,NW) = btemp(L,NW) + 0.9999*(TRAY(K,NW) + TRAY(K+1,NW)) ! JVO 2017 : does this 0.999 is really meaningful ?
+      btemp(L,NW) = btemp(L,NW) + TRAY(K,NW) + TRAY(K+1,NW)
       COSBV(L,NW,1:L_NGAUSS) = atemp(L,NW)/btemp(L,NW)
    END DO
@@ -568 +567 @@
       L = L_NLAYRAD
       K = 2*L+1
-           DTAUV(L,nw,ng) = DTAUKV(K,NW,NG)
+      DTAUV(L,nw,ng) = DTAUKV(K,NW,NG)
       WBARV(L,NW,NG) = ctemp(L,NW) / DTAUV(L,NW,NG)
    END DO