Changeset 3966 for trunk/LMDZ.PLUTO/libf

Timestamp:

Nov 20, 2025, 2:17:48 PM (3 weeks ago)

Author:

debatzbr

Message:

Pluto PCM: Add sanity check in the radiatif transfer (usefull for debug mode).
BBT

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• aeroptproperties.F90 (modified) (4 diffs)
• callcorrk.F90 (modified) (1 diff)

trunk/LMDZ.PLUTO/libf/phypluto/aeroptproperties.F90

@@ -656 +656 @@
                   qscatrefVISgrid(j,1,iaer)=qscatrefVISgrid(j,1,iaer) /  &
                                 normd(j,1,iaer,idomain)
-                  omegrefVISgrid(j,1,iaer)=qscatrefVISgrid(j,1,iaer) /   &
+                  if (qrefVISgrid(j,1,iaer).gt.0.d0) then
+                    omegrefVISgrid(j,1,iaer)=qscatrefVISgrid(j,1,iaer) / &
                                qrefVISgrid(j,1,iaer)
+                  else
+                    omegrefVISgrid(j,1,iaer)=0.
+                  endif
                   DO m=1,L_NSPECTV
                     qextVISgrid(j,1,m,iaer)=qextVISgrid(j,1,m,iaer) /    &
@@ -663 +667 @@
                     qscatVISgrid(j,1,m,iaer)=qscatVISgrid(j,1,m,iaer) /  &
                                 normd(j,1,iaer,idomain)
-                    gVISgrid(j,1,m,iaer)=gVISgrid(j,1,m,iaer) /          &
+                    if (qscatVISgrid(j,1,m,iaer).gt.0.d0) then
+                      gVISgrid(j,1,m,iaer)=gVISgrid(j,1,m,iaer) /        &
                                 qscatVISgrid(j,1,m,iaer) /               &
                                 normd(j,1,iaer,idomain)
-
-                    qsqrefVISgrid(j,1,m,iaer)=qextVISgrid(j,1,m,iaer) /  &
+                    else
+                      gVISgrid(j,1,m,iaer)=0.
+                    endif
+
+                    if (qrefVISgrid(j,1,iaer).gt.0.d0) then
+                      qsqrefVISgrid(j,1,m,iaer)=qextVISgrid(j,1,m,iaer) /&
                                 qrefVISgrid(j,1,iaer)
-                    omegVISgrid(j,1,m,iaer)=qscatVISgrid(j,1,m,iaer) /   &
+                    else
+                      qsqrefVISgrid(j,1,m,iaer)=0.
+                    endif
+                    if (qextVISgrid(j,1,m,iaer).gt.0.d0) then
+                      omegVISgrid(j,1,m,iaer)=qscatVISgrid(j,1,m,iaer) / &
                                 qextVISgrid(j,1,m,iaer)
+                    else
+                      omegVISgrid(j,1,m,iaer)=0.
+                    endif
                   ENDDO
                 ELSE                   ! INFRARED DOMAIN ----------
@@ -772 +788 @@
                   qscatrefIRgrid(j,1,iaer)=qscatrefIRgrid(j,1,iaer) /  &
                                 normd(j,1,iaer,idomain)
-                  omegrefIRgrid(j,1,iaer)=qscatrefIRgrid(j,1,iaer) /   &
+                  if (qrefIRgrid(j,1,iaer).gt.0.d0) then
+                    omegrefIRgrid(j,1,iaer)=qscatrefIRgrid(j,1,iaer) / &
                                qrefIRgrid(j,1,iaer)
+                  else
+                    omegrefIRgrid(j,1,iaer)=0.
+                  endif
                   DO m=1,L_NSPECTI
                     qextIRgrid(j,1,m,iaer)=qextIRgrid(j,1,m,iaer) /    &
@@ -779 +799 @@
                     qscatIRgrid(j,1,m,iaer)=qscatIRgrid(j,1,m,iaer) /  &
                                 normd(j,1,iaer,idomain)
-                    gIRgrid(j,1,m,iaer)=gIRgrid(j,1,m,iaer) /          &
+                    
+                    if (qscatIRgrid(j,1,m,iaer).gt.0.d0) then
+                      gIRgrid(j,1,m,iaer)=gIRgrid(j,1,m,iaer) /        &
                                 qscatIRgrid(j,1,m,iaer) /              &
                                 normd(j,1,iaer,idomain)
-
-                    qsqrefIRgrid(j,1,m,iaer)=qextIRgrid(j,1,m,iaer) /  &
+                    else
+                      gIRgrid(j,1,m,iaer)=0.
+                    endif
+
+                    if (qrefIRgrid(j,1,iaer).gt.0.d0) then
+                      qsqrefIRgrid(j,1,m,iaer)=qextIRgrid(j,1,m,iaer) /&
                                 qrefVISgrid(j,1,iaer)
-                    omegIRgrid(j,1,m,iaer)=qscatIRgrid(j,1,m,iaer) /   &
+                    else
+                      qsqrefIRgrid(j,1,m,iaer)=0.
+                    endif
+                    if (qextIRgrid(j,1,m,iaer).gt.0.d0) then
+                      omegIRgrid(j,1,m,iaer)=qscatIRgrid(j,1,m,iaer) / &
                                 qextIRgrid(j,1,m,iaer)
+                    else
+                      omegIRgrid(j,1,m,iaer)=0.
+                    endif
                   ENDDO
                 ENDIF                  ! --------------------------

trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90

@@ -723 +723 @@
                ! As 'dtau_aer' is at reference (visible) wavelenght we scale it as
                ! it will be multplied by qxi/v in optci/v
-               temp=dtau_aer(ig,L_NLAYRAD-k,iaer)/QREFvis3d(ig,L_NLAYRAD-k,iaer)
+               if (QREFvis3d(ig,L_NLAYRAD-k,iaer).gt.0.d0) then
+                  temp=dtau_aer(ig,L_NLAYRAD-k,iaer)/QREFvis3d(ig,L_NLAYRAD-k,iaer)
+               else
+                  temp=0.
+               endif
                tauaero(2*k+2,iaer)=max(temp*pweight,0.d0)
                tauaero(2*k+3,iaer)=max(temp-tauaero(2*k+2,iaer),0.d0)