Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ifs_blocking.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/ecrad/driver/ifs_blocking.F90 (modified) (12 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ifs_blocking.F90

@@ -329 +329 @@
   !$OMP PARALLEL DO SCHEDULE(RUNTIME)&
   !$OMP&PRIVATE(IB,IFLD)
-  do ib=1,ngpblks
-    do ifld=1,ifs_config%ifldstot
+  DO ib=1,ngpblks
+    DO ifld=1,ifs_config%ifldstot
       zrgp(:,ifld,ib) = 0._jprb
     enddo
@@ -351 +351 @@
     !$OMP PARALLEL DO SCHEDULE(RUNTIME)&
     !$OMP&PRIVATE(JRL,IBEG,IEND,IL,IB,JAER,JOFF,JLEV,JALB)
-    do jrl=1,ncol,nproma
+    DO jrl=1,ncol,nproma
 
       ibeg=jrl
@@ -363 +363 @@
       zrgp(1:il,ifs_config%iamu0,ib)  =  single_level%cos_sza(ibeg:iend)   ! cosine of solar zenith ang (mu0)
 
-      do jemiss=1,yderad%nlwemiss
+      DO jemiss=1,yderad%nlwemiss
         zrgp(1:il,ifs_config%iemiss+jemiss-1,ib)  =  single_level%lw_emissivity(ibeg:iend,jemiss)
       enddo
@@ -378 +378 @@
       ! zrgp(1:il,islon,ib)    = ???
 
-      do jalb=1,yderad%nsw
+      DO jalb=1,yderad%nsw
         zrgp(1:il,ifs_config%iald+jalb-1,ib)  =  single_level%sw_albedo(ibeg:iend,jalb)
       enddo
 
       if (allocated(single_level%sw_albedo_direct)) then
-        do jalb=1,yderad%nsw
+        DO jalb=1,yderad%nsw
           zrgp(1:il,ifs_config%ialp+jalb-1,ib)  =  single_level%sw_albedo_direct(ibeg:iend,jalb)
         end do
       else
-        do jalb=1,yderad%nsw
+        DO jalb=1,yderad%nsw
           zrgp(1:il,ifs_config%ialp+jalb-1,ib)  =  single_level%sw_albedo(ibeg:iend,jalb)
         end do
       end if
       
-      do jlev=1,nlev
+      DO jlev=1,nlev
         zrgp(1:il,ifs_config%iti+jlev-1,ib)   = temperature_fl(ibeg:iend,jlev) ! full level temperature
         zrgp(1:il,ifs_config%ipr+jlev-1,ib)   = pressure_fl(ibeg:iend,jlev) ! full level pressure
@@ -398 +398 @@
       enddo
 
-      do jlev=1,nlev
+      DO jlev=1,nlev
         zrgp(1:il,ifs_config%iwv+jlev-1,ib)   = gas%mixing_ratio(ibeg:iend,jlev,IH2O) ! this is already in MassMixingRatio units
         if (rad_config%do_clouds) then
@@ -421 +421 @@
       if (yderad%naermacc == 1) then
         joff=ifs_config%iaero
-        do jaer=1,rad_config%n_aerosol_types
-          do jlev=1,nlev
+        DO jaer=1,rad_config%n_aerosol_types
+          DO jlev=1,nlev
             zrgp(1:il,joff,ib) = aerosol%mixing_ratio(ibeg:iend,jlev,jaer)
             joff=joff+1
@@ -429 +429 @@
       endif
 
-      do jlev=1,nlev+1
+      DO jlev=1,nlev+1
         ! zrgp(1:il,ihpr+jlev-1,ib)  = ???
         zrgp(1:il,ifs_config%iaprs+jlev-1,ib) = thermodynamics%pressure_hl(ibeg:iend,jlev)
@@ -453 +453 @@
       ! local workaround variables for standalone input files
       if (rad_config%do_clouds) then
-        do jlev=1,nlev
+        DO jlev=1,nlev
           ! missing full-level temperature and pressure as well as land-sea-mask
           zrgp(1:il,ifs_config%ire_liq+jlev-1,ib) = cloud%re_liq(ibeg:iend,jlev)
           zrgp(1:il,ifs_config%ire_ice+jlev-1,ib) = cloud%re_ice(ibeg:iend,jlev)
         enddo
-        do jlev=1,nlev-1
+        DO jlev=1,nlev-1
           ! for the love of it, I can't figure this one out. Probably to do with
           ! my crude approach of setting PGEMU?
@@ -465 +465 @@
         if(present(iseed)) iseed(1:il,ib) = single_level%iseed(ibeg:iend)
       else
-        do jlev=1,nlev
+        DO jlev=1,nlev
           ! missing full-level temperature and pressure as well as land-sea-mask
           zrgp(1:il,ifs_config%ire_liq+jlev-1,ib) = 0._jprb
           zrgp(1:il,ifs_config%ire_ice+jlev-1,ib) = 0._jprb
         enddo
-        do jlev=1,nlev-1
+        DO jlev=1,nlev-1
           zrgp(1:il,ifs_config%ioverlap+jlev-1,ib) = 0._jprb
         enddo
@@ -531 +531 @@
     !$OMP PARALLEL DO SCHEDULE(RUNTIME)&
     !$OMP&PRIVATE(JRL,IBEG,IEND,IL,IB,JLEV,JG)
-    do jrl=1,ncol,nproma
+    DO jrl=1,ncol,nproma
       ibeg=jrl
       iend=min(ibeg+nproma-1,ncol)
@@ -537 +537 @@
       ib=(jrl-1)/nproma+1
 
-      do jlev=1,nlev+1
+      DO jlev=1,nlev+1
         flux%sw_up(ibeg:iend,jlev) = zrgp(1:il,ifs_config%ifrso+jlev-1,ib)
         flux%lw_up(ibeg:iend,jlev) = zrgp(1:il,ifs_config%ifrth+jlev-1,ib)
@@ -561 +561 @@
       emissivity_out(ibeg:iend) = zrgp(1:il,ifs_config%iemit,ib)
       if (yradiation%yrerad%lapproxswupdate) then
-        do jg=1,yradiation%yrerad%nsw
+        DO jg=1,yradiation%yrerad%nsw
           flux_diffuse_band(ibeg:iend,jg) = zrgp(1:il,ifs_config%iswdiffuseband+jg-1,ib)
           flux_direct_band(ibeg:iend,jg) = zrgp(1:il,ifs_config%iswdirectband+jg-1,ib)