Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmd/cpl_mod.F90

Timestamp:

Jul 22, 2024, 9:29:09 PM (4 months ago)

Author:

abarral

Message:

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/cpl_mod.F90 (modified) (52 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/cpl_mod.F90

@@ -1 @@
-!
+
 MODULE cpl_mod
-!
+
 ! This module excahanges and transforms all fields that should be recieved or sent to 
 ! coupler. The transformation of the fields are done from the grid 1D-array in phylmd 
 ! to the regular 2D grid accepted by the coupler. Cumulation of the fields for each 
 ! timestep is done in here. 
-!
+
 ! Each type of surface that recevie fields from the coupler have a subroutine named 
 ! cpl_receive_XXX_fields and each surface that have fields to be sent to the coupler 
 ! have a subroutine named cpl_send_XXX_fields.
-!
+
 !*************************************************************************************
 
@@ -129 +129 @@
   
 CONTAINS
-!
+
 !************************************************************************************
-!
+
   SUBROUTINE cpl_init(dtime, rlon, rlat)
     USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, fco2_ocn_day
@@ -178 +178 @@
 !*************************************************************************************
 ! Calculate coupling period
-!
+
 !*************************************************************************************
      
@@ -191 +191 @@
 !*************************************************************************************
 ! Allocate variables
-!
+
 !*************************************************************************************
     error = 0
@@ -324 +324 @@
 !*************************************************************************************
 ! Initialize the allocated varaibles
-!
+
 !*************************************************************************************
     DO ig = 1, klon
@@ -332 +332 @@
 !*************************************************************************************
 ! Initialize coupling
-!
+
 !*************************************************************************************
     idtime = INT(dtime)
@@ -341 +341 @@
 !*************************************************************************************
 ! initialize NetCDF output
-!
+
 !*************************************************************************************
     IF (is_sequential) THEN
@@ -390 +390 @@
 !*************************************************************************************
 ! compatibility test
-!
+
 !*************************************************************************************
     IF (carbon_cycle_cpl .AND. version_ocean=='opa8') THEN
@@ -398 +398 @@
 
   END SUBROUTINE cpl_init
-  
-!
-!*************************************************************************************
-!
- 
+
+!*************************************************************************************
+
   SUBROUTINE cpl_receive_frac(itime, dtime, pctsrf, is_modified)
 ! This subroutine receives from coupler for both ocean and seaice
@@ -442 +440 @@
 ! Start calculation
 ! Get fields from coupler
-!
+
 !*************************************************************************************
 
@@ -510 +508 @@
 !  Transform seaice fraction (read_sic : ocean-seaice mask) into global 
 !  fraction (pctsrf : ocean-seaice-land-landice mask)
-!
+
 !*************************************************************************************
        CALL cpl2gath(read_sic, read_sic1D, klon, unity)
@@ -530 +528 @@
   END SUBROUTINE cpl_receive_frac
 
-!
-!*************************************************************************************
-!
+!*************************************************************************************
 
   SUBROUTINE cpl_receive_ocean_fields(knon, knindex, tsurf_new, u0_new, &
        v0_new, sss)
-!
+
 ! This routine returns the field for the ocean that has been read from the coupler
 ! (done earlier with cpl_receive_frac). The field is the temperature.
 ! The temperature is transformed into 1D array with valid points from index 1 to knon.
-!
+
     USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, fco2_ocn_day
     USE indice_sol_mod
@@ -568 +564 @@
 !*************************************************************************************
 ! Transform read_sst into compressed 1D variable tsurf_new
-!
+
 !*************************************************************************************
     CALL cpl2gath(read_sst, tsurf_new, knon, knindex)
@@ -579 +575 @@
 ! Transform read_co2 into uncompressed 1D variable fco2_ocn_day added directly in 
 ! the module carbon_cycle_mod
-!
+
 !*************************************************************************************
     IF (carbon_cycle_cpl) THEN
@@ -591 +587 @@
 ! The fields received from the coupler have to be weighted with the fraction of ocean 
 ! in relation to the total sea-ice+ocean
-!
+
 !*************************************************************************************
     DO i=1, knon
@@ -599 +595 @@
   END SUBROUTINE cpl_receive_ocean_fields
 
-!
-!*************************************************************************************
-!
+!*************************************************************************************
 
   SUBROUTINE cpl_receive_seaice_fields(knon, knindex, &
        tsurf_new, alb_new, u0_new, v0_new)
-!
+
 ! This routine returns the fields for the seaice that have been read from the coupler
 ! (done earlier with cpl_receive_frac). These fields are the temperature and 
 ! albedo at sea ice surface and fraction of sea ice.
 ! The fields are transformed into 1D arrays with valid points from index 1 to knon. 
-!
 
 ! Input arguments
@@ -631 +624 @@
 !*************************************************************************************
 ! Transform fields read from coupler from 2D into compressed 1D variables
-!
+
 !*************************************************************************************
     CALL cpl2gath(read_sit, tsurf_new, knon, knindex)
@@ -642 +635 @@
 ! The fields received from the coupler have to be weighted with the sea-ice 
 ! concentration (in relation to the total sea-ice + ocean).
-!
+
 !*************************************************************************************
     DO i= 1, knon
@@ -651 +644 @@
   END SUBROUTINE cpl_receive_seaice_fields
 
-!
-!*************************************************************************************
-!
+!*************************************************************************************
 
   SUBROUTINE cpl_send_ocean_fields(itime, knon, knindex, &
@@ -715 +706 @@
 ! Start calculation
 ! The ocean points are saved with second array index=1
-!
+
 !*************************************************************************************
     cpl_index = 1
@@ -721 +712 @@
 !*************************************************************************************
 ! Reset fields to zero in the beginning of a new coupling period 
-!
+
 !*************************************************************************************
     IF (MOD(itime, nexca) == 1) THEN
@@ -751 +742 @@
 !*************************************************************************************
 ! Cumulate at each time-step
-!
+
 !*************************************************************************************    
     DO ig = 1, knon
@@ -802 +793 @@
 ! fields are transformed to the 2D grid. 
 ! No sending to the coupler (it is done from cpl_send_seaice_fields).
-!
+
 !*************************************************************************************
     IF (MOD(itime, nexca) == 0) THEN
@@ -915 +906 @@
   END SUBROUTINE cpl_send_ocean_fields
 
-!
-!*************************************************************************************
-!
+!*************************************************************************************
 
   SUBROUTINE cpl_send_seaice_fields(itime, dtime, knon, knindex, &
@@ -924 +913 @@
        precip_rain, precip_snow, evap, tsurf, fder, albsol, taux, tauy,&
        sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol)
-!
+
 ! This subroutine cumulates some fields for each time-step during a coupling 
 ! period. At last time-step in a coupling period the fields are transformed to the 
 ! grid accepted by the coupler. All fields for all types of surfaces are sent to
 ! the coupler.
-!
+
     USE carbon_cycle_mod, ONLY : carbon_cycle_cpl
     USE indice_sol_mod
@@ -962 +951 @@
 ! Start calulation
 ! The sea-ice points are saved with second array index=2
-!
+
 !*************************************************************************************
     cpl_index = 2
@@ -968 +957 @@
 !*************************************************************************************
 ! Reset fields to zero in the beginning of a new coupling period 
-!
+
 !*************************************************************************************
     IF (MOD(itime, nexca) == 1) THEN
@@ -988 +977 @@
 !*************************************************************************************
 ! Cumulate at each time-step
-!
+
 !*************************************************************************************    
     DO ig = 1, knon
@@ -1022 +1011 @@
 ! If the time-step corresponds to the end of coupling period the 
 ! fields are transformed to the 2D grid and all fields are sent to coupler.
-!
+
 !*************************************************************************************
     IF (MOD(itime, nexca) == 0) THEN
@@ -1119 +1108 @@
   END SUBROUTINE cpl_send_seaice_fields
 
-!
-!*************************************************************************************
-!
+!*************************************************************************************
 
   SUBROUTINE cpl_send_land_fields(itime, knon, knindex, rriv_in, rcoa_in)
-!
+
 ! This subroutine cumulates some fields for each time-step during a coupling 
 ! period. At last time-step in a coupling period the fields are transformed to the 
 ! grid accepted by the coupler. No sending to the coupler will be done from here 
 ! (it is done in cpl_send_seaice_fields).
-!
+
     USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
 
@@ -1148 +1135 @@
 ! Rearrange fields in 2D variables 
 ! First initialize to zero to avoid unvalid points causing problems
-!
+
 !*************************************************************************************
 !$OMP MASTER
@@ -1159 +1146 @@
 !*************************************************************************************
 ! Reset cumulated fields to zero in the beginning of a new coupling period 
-!
+
 !*************************************************************************************
     IF (MOD(itime, nexca) == 1) THEN
@@ -1170 +1157 @@
 !*************************************************************************************
 ! Cumulate : Following fields should be cumulated at each time-step
-!
+
 !*************************************************************************************    
 !$OMP MASTER
@@ -1179 +1166 @@
   END SUBROUTINE cpl_send_land_fields
 
-!
-!*************************************************************************************
-!
+!*************************************************************************************
 
   SUBROUTINE cpl_send_landice_fields(itime, knon, knindex, rlic_in, rlic_in_frac)
@@ -1187 +1172 @@
 ! during a coupling period. This routine will not send to coupler. Sending 
 ! will be done in cpl_send_seaice_fields.
-!
 
     USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
@@ -1208 +1192 @@
 ! Rearrange field in a 2D variable 
 ! First initialize to zero to avoid unvalid points causing problems
-!
+
 !*************************************************************************************
 !$OMP MASTER
@@ -1217 +1201 @@
 !*************************************************************************************
 ! Reset field to zero in the beginning of a new coupling period 
-!
+
 !*************************************************************************************
     IF (MOD(itime, nexca) == 1) THEN
@@ -1227 +1211 @@
 !*************************************************************************************
 ! Cumulate : Melting ice should be cumulated at each time-step
-!
+
 !*************************************************************************************    
 !$OMP MASTER
@@ -1235 +1219 @@
   END SUBROUTINE cpl_send_landice_fields
 
-!
-!*************************************************************************************
-!
+!*************************************************************************************
 
   SUBROUTINE cpl_send_all(itime, dtime, pctsrf, lafin, rlon, rlat)
@@ -1243 +1225 @@
 ! This subroutine should be executed after calculations by the last surface(sea-ice),
 ! all calculations at the different surfaces have to be done before. 
-!    
+
     USE surface_data
     USE carbon_cycle_mod, ONLY : carbon_cycle_cpl
@@ -1253 +1235 @@
 
 ! Some includes
-!    
+
 ! Input arguments
 !*************************************************************************************
@@ -1295 +1277 @@
 ! All fields are stored in a table tab_flds(:,:,:)
 ! First store the fields which are already on the right format
-!
+
 !*************************************************************************************
 !$OMP MASTER
@@ -1329 +1311 @@
 !*************************************************************************************
 ! Transform the fraction of sub-surfaces from 1D to 2D array
-!
+
 !*************************************************************************************
     pctsrf2D(:,:,:) = 0.
@@ -1342 +1324 @@
 ! Calculate the average calving per latitude
 ! Store calving in tab_flds(:,:,19)
-! 
+
 !*************************************************************************************      
     IF (is_omp_root) THEN
@@ -1403 +1385 @@
 ! Calculate total flux for snow, rain and wind with weighted addition using the 
 ! fractions of ocean and seaice.
-!
+
 !*************************************************************************************    
        ! fraction oce+seaice
@@ -1507 +1489 @@
 !*************************************************************************************
 ! NetCDF output of all fields just before sending to coupler.
-!
+
 !*************************************************************************************
     IF (is_sequential) THEN
@@ -1517 +1499 @@
 !*************************************************************************************
 ! Send the table of all fields
-!
+
 !*************************************************************************************
     time_sec=(itime-1)*dtime
@@ -1528 +1510 @@
 !*************************************************************************************
 ! Finish with some dellocate
-!
+
 !*************************************************************************************  
     sum_error=0
@@ -1555 +1537 @@
     
   END SUBROUTINE cpl_send_all
-!
-!*************************************************************************************
-!
+
+!*************************************************************************************
+
   SUBROUTINE cpl2gath(champ_in, champ_out, knon, knindex)
   USE mod_phys_lmdz_para
 ! Cette routine transforme un champs de la grille 2D recu du coupleur sur la grille 
 ! 'gathered' (la grille physiq comprime).
-!
-! 
+
+
 ! input:         
 !   champ_in     champ sur la grille 2D
 !   knon         nombre de points dans le domaine a traiter
 !   knindex      index des points de la surface a traiter
-!
+
 ! output:
 !   champ_out    champ sur la grille 'gatherd'
-!
+
     USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
 
@@ -1588 +1570 @@
 
 !*************************************************************************************
-!
-    
 
 ! Transform from 2 dimensions (nbp_lon,jj_nb) to 1 dimension (klon)
@@ -1605 +1585 @@
 
   END SUBROUTINE cpl2gath
-!
-!*************************************************************************************
-!
+
+!*************************************************************************************
+
   SUBROUTINE gath2cpl(champ_in, champ_out, knon, knindex)
   USE mod_phys_lmdz_para
 ! Cette routine ecrit un champ 'gathered' sur la grille 2D pour le passer
 ! au coupleur.
-!
+
 ! input:         
 !   champ_in     champ sur la grille gathere        
 !   knon         nombre de points dans le domaine a traiter
 !   knindex      index des points de la surface a traiter
-!
+
 ! output:
 !   champ_out    champ sur la grille 2D
-!
+
     USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
     
@@ -1658 +1638 @@
     
   END SUBROUTINE gath2cpl
-!
-!*************************************************************************************
-!
+
+!*************************************************************************************
+
 END MODULE cpl_mod