Changeset 2459 for trunk

Timestamp:

Feb 16, 2021, 12:51:04 AM (4 years ago)

Author:

aslmd

Message:

MESOSCALE Mars (and actually GCM too): solved inconsistency in types with riceco2 in a cascade of subroutines. harmless when everything is compiled double precision, harmful when everything is compiled simple precision. appeared between r2446 and r2453 included.

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• callradite_mod.F (modified) (1 diff)
• rocketduststorm_mod.F90 (modified) (1 diff)
• topmons_mod.F90 (modified) (1 diff)
• updatereffrad_mod.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/callradite_mod.F

@@ -203 +203 @@
       REAL,INTENT(OUT) :: rice(ngrid,nlayer)   ! Ice geometric mean radius (m)
       REAL,INTENT(OUT) :: nuice(ngrid,nlayer)  ! Estimated effective variance
-      REAL,INTENT(OUT) :: riceco2(ngrid,nlayer) ! CO2 ice mean radius(m)
+      double precision,INTENT(OUT) :: riceco2(ngrid,nlayer) ! CO2 ice mean radius(m)
       REAL,INTENT(OUT) :: nuiceco2(ngrid,nlayer) ! Effective variance
       REAL,INTENT(IN) :: co2ice(ngrid)           ! co2 ice surface layer (kg.m-2)

trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90

@@ -172 +172 @@
       REAL  rice(ngrid,nlayer)
       REAL  nuice(ngrid,nlayer)
-      REAL  riceco2(ngrid,nlayer)
+      DOUBLE PRECISION  riceco2(ngrid,nlayer)
       REAL  nuiceco2(ngrid,nlayer)
 

trunk/LMDZ.MARS/libf/phymars/topmons_mod.F90

@@ -126 +126 @@
       REAL  rice(ngrid,nlayer)
       REAL  nuice(ngrid,nlayer)
-      REAL  riceco2(ngrid,nlayer)
+      DOUBLE PRECISION  riceco2(ngrid,nlayer)
       REAL  nuiceco2(ngrid,nlayer)
 

trunk/LMDZ.MARS/libf/phymars/updatereffrad_mod.F

@@ -71 +71 @@
       REAL, INTENT(in) :: tauscaling(ngrid)         ! Convertion factor for qccn and Nccn
 c     CO2 ice mean radius (m)
-      REAL, INTENT(out) :: riceco2(ngrid,nlayer) ! co2 ice radius
+      double precision, INTENT(out) :: riceco2(ngrid,nlayer) ! co2 ice radius
       REAL, INTENT(out) :: nuiceco2(ngrid,nlayer)
 
@@ -201 +201 @@
           DO l=1,nlayer
             DO ig=1,ngrid
-              call updaterice_microco2(pq(ig,l,igcm_co2_ice),
-     &                              pq(ig,l,igcm_ccnco2_mass),
-     &                              pq(ig,l,igcm_ccnco2_number),
+              call updaterice_microco2(dble(pq(ig,l,igcm_co2_ice)),
+     &                              dble(pq(ig,l,igcm_ccnco2_mass)),
+     &                              dble(pq(ig,l,igcm_ccnco2_number)),
      &                              tauscaling(ig),riceco2(ig,l),
      &                              rhocloudco2(ig,l))