Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/modis_simulator.F90

Timestamp:

Aug 2, 2024, 2:12:03 PM (7 weeks ago)

Author:

abarral

Message:

Add missing klon on strataer_emiss_mod.F90
Correct various missing explicit declarations
Replace tabs by spaces (tabs are not part of the fortran charset)
Continue cleaning modules
Removed unused arguments and variables

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/modis_simulator.F90 (modified) (6 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/modis_simulator.F90

@@ -163 +163 @@
     cloudMask = retrievedTau(1:nSubCols) >= min_OpticalThickness
     
-    do i = 1, nSubCols
+    DO i = 1, nSubCols
        if(cloudMask(i)) then 
           ! ##################################################################################
@@ -380 +380 @@
     reffIceWRK(1:nPoints,1:nSubCols) = merge(particle_size,R_UNDEF,iceCloudMask)
     reffLiqWRK(1:nPoints,1:nSubCols) = merge(particle_size,R_UNDEF,waterCloudMask)
-    do j=1,nPoints
+    DO j=1,nPoints
 
        ! Fill clear and optically thin subcolumns with fill
@@ -439 +439 @@
     !   layers and use the trapezoidal rule.
     totalTau = 0._wp; totalProduct = 0._wp
-    do i = 2, size(tauIncrement)
+    DO i = 2, size(tauIncrement)
       if(totalTau + tauIncrement(i) > tauLimit) then 
         deltaX = tauLimit - totalTau
@@ -479 +479 @@
     ! Find the extinction-weighted value of f(tau), assuming constant f within each layer
     totalTau = 0._wp; totalProduct = 0._wp
-    do i = 1, size(tauIncrement)
+    DO i = 1, size(tauIncrement)
       if(totalTau + tauIncrement(i) > tauLimit) then 
         deltaX       = tauLimit - totalTau
@@ -712 +712 @@
     cloudMask(1:nLevels) = tau(1:nLevels) > 0. 
     cloudIndicies = pack((/ (i, i = 1, nLevels) /), mask = cloudMask) 
-    do i = 1, size(cloudIndicies)
+    DO i = 1, size(cloudIndicies)
        call two_stream(tau(cloudIndicies(i)), g(cloudIndicies(i)), w0(cloudIndicies(i)), Refl(i), Trans(i))
     end do
@@ -892 +892 @@
     Tran_cumulative(1) = Tran(1)    
     
-    do i=2, npts
+    DO i=2, npts
        ! place (add) previous combined layer(s) reflectance on top of layer i, w/black surface (or ignoring surface):
        Refl_cumulative(i) = Refl_cumulative(i-1) + Refl(i)*(Tran_cumulative(i-1)**2)/(1._wp - Refl_cumulative(i-1) * Refl(i))