Changeset 2939 for trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90

Timestamp:

Apr 14, 2023, 11:36:06 AM (20 months ago)

Author:

llange

Message:

PEM
Fixing bug when reading the mass adsorded written in startfi_PEM (only the mass for islope = 1 was read)
I've cleaned some print * I've introduced in my last commit
LL

File:

• trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90 (modified) (3 diffs)

trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90

@@ -90 +90 @@
 
 write(*,*)'Is start PEM?',startpem_file
-startpem_file = .true.
 
 !1. Run
@@ -275 +274 @@
 
 !4. CO2 & H2O Adsorption
+
  if(adsorption_pem) then
   DO islope=1,nslope
@@ -281 +281 @@
     if((.not.found)) then
        m_co2_regolith_phys(:,:,:) = 0.
+       exit
     endif
-    exit
+    
   ENDDO
 
   DO islope=1,nslope
    write(num,fmt='(i2.2)') islope
-   call get_field("mh2o_reg_ads_slope"//num,m_co2_regolith_phys(:,:,islope),found2)
+   call get_field("mh2o_reg_ads_slope"//num,m_h2o_regolith_phys(:,:,islope),found2)
     if((.not.found2)) then
        m_h2o_regolith_phys(:,:,:) = 0.
+      exit
     endif
-    exit
+    
   ENDDO