Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/cosp_optics.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/cosp2/cosp_optics.F90 (modified) (10 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/cosp_optics.F90

@@ -69 +69 @@
     
     varOUT(1:dim1,1:dim2,1:dim3) = 0._wp
-    do j=1,dim2
+    DO j=1,dim2
        where(flag(:,j,:) .eq. 1)
           varOUT(:,j,:) = varIN2
@@ -133 +133 @@
     
     
-    do i=1,npoints
+    DO i=1,npoints
        where(cloudIce(i,:, :) <= 0.) 
           fracL(:, :) = 1._wp
@@ -168 +168 @@
     w0(1:nPoints,1:nSubCols,1:nLevels) = 0._wp
     
-    do j =1,nPoints
-       do i=1,nSubCols
+    DO j =1,nPoints
+       DO i=1,nSubCols
           water_g(1:nLevels)  = get_g_nir(  phaseIsLiquid, sizeLIQ(j,i,1:nLevels)) 
           water_w0(1:nLevels) = get_ssa_nir(phaseIsLiquid, sizeLIQ(j,i,1:nLevels))
@@ -187 +187 @@
     
     ! Compute the total optical thickness and the proportion due to liquid in each cell
-    do i=1,npoints
+    DO i=1,npoints
        where(tauLIQ(i,1:nSubCols,1:nLevels) + tauICE(i,1:nSubCols,1:nLevels) > 0.) 
           fracLIQ(i,1:nSubCols,1:nLevels) = tauLIQ(i,1:nSubCols,1:nLevels)/ &
@@ -317 +317 @@
     ! Altitude at half pressure levels:
     zheight(1:npoints,nlev+1) = 0._wp
-    do k=nlev,1,-1
+    DO k=nlev,1,-1
        zheight(1:npoints,k) = zheight(1:npoints,k+1) &
             -(presf(1:npoints,k)-presf(1:npoints,k+1))/(rhoair(1:npoints,k)*grav)
@@ -349 +349 @@
     ! ##############################################################################
     ! Polynomials kp_lidar derived from Mie theory
-    do i = 1, npart
+    DO i = 1, npart
        where (rad_part(1:npoints,1:nlev,i) .gt. 0.0)
           kp_part(1:npoints,1:nlev,i) = &
@@ -363 +363 @@
     
     ! Loop over all subcolumns
-    do icol=1,ncolumns
+    DO icol=1,ncolumns
        ! ##############################################################################
        ! Mixing ratio particles in each subcolum
@@ -376 +376 @@
        ! ##############################################################################
        ! Alpha of particles in each subcolumn:
-       do i = 1, npart
+       DO i = 1, npart
           where (rad_part(1:npoints,1:nlev,i) .gt. 0.0)
              alpha_part(1:npoints,icol,1:nlev,i) = 3._wp/4._wp * Qscat &
@@ -388 +388 @@
        ! Optical thicknes
        tau_part(1:npoints,icol,1:nlev,1:npart) = rdiffm * alpha_part(1:npoints,icol,1:nlev,1:npart)
-       do i = 1, npart
+       DO i = 1, npart
           ! Optical thickness of each layer (particles)
           tau_part(1:npoints,icol,1:nlev,i) = tau_part(1:npoints,icol,1:nlev,i) &
                & * (zheight(1:npoints,1:nlev)-zheight(1:npoints,2:nlev+1) )
           ! Optical thickness from TOA to layer k (particles)
-          do k=2,nlev
+          DO k=2,nlev
              tau_part(1:npoints,icol,k,i) = tau_part(1:npoints,icol,k,i) + tau_part(1:npoints,icol,k-1,i)
           enddo
@@ -434 +434 @@
     tautot_S_liq(1:npoints,1:ncolumns) = 0._wp
     tautot_S_ice(1:npoints,1:ncolumns) = 0._wp
-    do icol=1,ncolumns    
+    DO icol=1,ncolumns
        tautot_S_liq(1:npoints,icol) = tautot_S_liq(1:npoints,icol)+tau_part(1:npoints,icol,nlev,1)+tau_part(1:npoints,icol,nlev,3)
        tautot_S_ice(1:npoints,icol) = tautot_S_ice(1:npoints,icol)+tau_part(1:npoints,icol,nlev,2)+tau_part(1:npoints,icol,nlev,4)