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

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

@@ -190 +190 @@
      real(wp),dimension(:),pointer :: &
           isccp_totalcldarea => null(), & ! The fraction of model grid box columns with cloud 
-                                          ! somewhere in them. (%)
+                                     ! somewhere in them. (%)
           isccp_meantb => null(),       & ! Mean all-sky 10.5 micron brightness temperature. (K)
           isccp_meantbclr => null(),    & ! Mean clear-sky 10.5 micron brightness temperature. (K)
@@ -203 +203 @@
                                           ! the 49 ISCCP D level cloud types. (%)
      
-     ! MISR outptus                         
+     ! MISR outptus
      real(wp),dimension(:,:,:),pointer ::   & !
           misr_fq => null()          ! Fraction of the model grid box covered by each of the MISR 
@@ -211 +211 @@
      real(wp),dimension(:),pointer ::   & !
           misr_meanztop => null(), & ! Mean MISR cloud top height
-          misr_cldarea => null()     ! Mean MISR cloud cover area                                   
-
-     ! MODIS outptus                
+          misr_cldarea => null()     ! Mean MISR cloud cover area
+
+     ! MODIS outptus
      real(wp),pointer,dimension(:) ::      & !  
           modis_Cloud_Fraction_Total_Mean => null(),       & ! L3 MODIS retrieved cloud fraction (total) 
@@ -233 +233 @@
           modis_Ice_Water_Path_Mean => null()                ! L3 MODIS retrieved ice water path
      real(wp),pointer,dimension(:,:,:) ::  &
-          modis_Optical_Thickness_vs_Cloud_Top_Pressure => null(), & ! Tau/Pressure joint histogram                                 
+          modis_Optical_Thickness_vs_Cloud_Top_Pressure => null(), & ! Tau/Pressure joint histogram
           modis_Optical_Thickness_vs_ReffICE => null(),            & ! Tau/ReffICE joint histogram
           modis_Optical_Thickness_vs_ReffLIQ => null()               ! Tau/ReffLIQ joint histogram
@@ -239 +239 @@
      ! RTTOV outputs
      real(wp),pointer :: &
-          rttov_tbs(:,:) => null() ! Brightness Temperature         
+          rttov_tbs(:,:) => null() ! Brightness Temperature
      
   end type cosp_outputs
@@ -683 +683 @@
        allocate(parasolPix_refl(parasolIN%Npoints,parasolIN%Ncolumns,PARASOL_NREFL))
        ! Call simulator
-       do icol=1,parasolIN%Ncolumns
+       DO icol=1,parasolIN%Ncolumns
           call parasol_subcolumn(parasolIN%npoints, PARASOL_NREFL,                       &
                                  parasolIN%tautot_S_liq(1:parasolIN%Npoints,icol),       &
@@ -700 +700 @@
        ! Allocate space for local variables
        allocate(cloudsatDBZe(cloudsatIN%Npoints,cloudsatIN%Ncolumns,cloudsatIN%Nlevels))
-       do icol=1,cloudsatIN%ncolumns
+       DO icol=1,cloudsatIN%ncolumns
           call quickbeam_subcolumn(cloudsatIN%rcfg,cloudsatIN%Npoints,cloudsatIN%Nlevels,&
                                    cloudsatIN%hgt_matrix/1000._wp,                       &
@@ -721 +721 @@
                    modisRetrievedCloudTopPressure(modisIN%nSunlit,modisIN%nColumns))
           ! Call simulator
-          do i = 1, modisIN%nSunlit
+          DO i = 1, modisIN%nSunlit
              call modis_subcolumn(modisIN%Ncolumns,modisIN%Nlevels,modisIN%pres(i,:),    &
                                   modisIN%tau(int(modisIN%sunlit(i)),:,:),               &
@@ -1367 +1367 @@
        ! Other grid requested. Constant vertical spacing with top at 20 km
        if (.not. luseCSATvgrid) zstep = 20000._wp/Nvgrid
-       do i=1,Nvgrid
+       DO i=1,Nvgrid
           vgrid_zl(Nlvgrid-i+1) = (i-1)*zstep
           vgrid_zu(Nlvgrid-i+1) = i*zstep
@@ -1726 +1726 @@
     endif
 
-        ! RTTOV Inputs
+    ! RTTOV Inputs
     if (cospgridIN%zenang .lt. -90. .OR. cospgridIN%zenang .gt. 90) then
        nError=nError+1