Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/mod_cosp.F90

Timestamp:

Aug 2, 2024, 2:12:03 PM (3 months 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/cosp/mod_cosp.F90 (modified) (7 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/mod_cosp.F90

@@ -68 +68 @@
   integer,dimension(2) :: ix,iy
   logical :: reff_zero
-  real :: maxp,minp
+  REAL :: maxp,minp
   integer,dimension(:),allocatable :: & ! Dimensions nPoints
                   seed    !  It is recommended that the seed is set to a different value for each model
@@ -95 +95 @@
 
 !++++++++++ Depth of model layers ++++++++++++
-  do i=1,Nlevels-1
+  DO i=1,Nlevels-1
     gbx%dlev(:,i) = gbx%zlev_half(:,i+1) - gbx%zlev_half(:,i)
   enddo
@@ -195 +195 @@
         Niter = gbx%Npoints/gbx%Npoints_it ! Integer division
         if (Niter*gbx%Npoints_it < gbx%Npoints) Niter = Niter + 1
-        do i=1,Niter
+        DO i=1,Niter
             i_first = (i-1)*gbx%Npoints_it + 1
             i_last  = i_first + gbx%Npoints_it - 1
@@ -409 +409 @@
         
         ! Precipitation fraction
-        do j=1,Npoints,1
-        do k=1,Nlevels,1
-            do i=1,Ncolumns,1
+        DO j=1,Npoints,1
+        DO k=1,Nlevels,1
+            DO i=1,Ncolumns,1
                 if (sgx%frac_out (j,i,Nlevels+1-k) == I_LSC) frac_ls(j,k)=frac_ls(j,k)+1.
                 if (sgx%frac_out (j,i,Nlevels+1-k) == I_CVC) frac_cv(j,k)=frac_cv(j,k)+1.
@@ -434 +434 @@
         else
             ! This is done within a loop (unvectorized) over nPoints to save memory
-            do j=1,Npoints
+            DO j=1,Npoints
                 sgx%frac_out(j,:,1:Nlevels)  = sgx%frac_out(j,:,Nlevels:1:-1)
                 sgx%prec_frac(j,:,1:Nlevels) = sgx%prec_frac(j,:,Nlevels:1:-1)
@@ -445 +445 @@
         ! Populate the subgrid arrays
         call construct_cosp_sghydro(Npoints,Ncolumns,Nlevels,Nhydro,sghydro)
-        do k=1,Ncolumns
+        DO k=1,Ncolumns
             !--------- Mixing ratios for clouds and Reff for Clouds and precip -------
             column_frac_out => sgx%frac_out(:,k,:)
@@ -498 +498 @@
         enddo
         ! convert the mixing ratio and precipitation flux from gridbox mean to the fraction-based values
-        do k=1,Nlevels
-            do j=1,Npoints
+        DO k=1,Nlevels
+            DO j=1,Npoints
                 !--------- Clouds -------
                 if (frac_ls(j,k) .ne. 0.) then