Changeset 2951

Timestamp:

Apr 27, 2023, 4:37:46 PM (20 months ago)

Author:

emillour

Message:

Mars PCM:
Fix OpenMP bug on inertiesoil, unnecessary loop and a line
too long for fixed fortran format.
Minor fix for picky gfortran compiler in testphys1d (no .eq. for logicals!)
EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/comsoil_h.F90 (modified) (1 diff)
• libf/phymars/dyn1d/testphys1d.F (modified) (1 diff)
• libf/phymars/soil_settings.F (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -4014 +4014 @@
 written in the startfile).
 Introduce the possibility to prescribe subsurface ice in the 1D model
+
+== 27/04/2023 == EM
+Fix OpenMP bug on inertiesoil, unnecessary loop and a line
+too long for fixed fortran format.
+Minor fix for picky gfortran compiler in testphys1d (no .eq. for logicals!)

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

@@ -15 +15 @@
        !                 soil_settings.F)
 
-!$OMP THREADPRIVATE(layer,mlayer,inertiedat,volcapa)
+!$OMP THREADPRIVATE(layer,mlayer,inertiedat,inertiesoil,volcapa)
 
   ! variables (FC: built in firstcall in soil.F)

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

@@ -992 +992 @@
       write_prof=.false.
       call getin("write_prof",write_prof)
-      IF (write_prof.eq..true.) then
+      IF (write_prof) then
               DO iq = 1, nq
                 call writeprofile(nlayer,q(:,iq),noms(iq),iq)

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

@@ -320 +320 @@
       if (.not.ok) then
         write(*,*)'soil_settings: Field <inertiesoil> not found!'
-        write(*,*)' => Building <inertiesoil> from surface values <inertiedat>'
-        do iloop=1,nsoil
-      do islope=1,nslope
+        write(*,*)' => Building <inertiesoil> from surface values'//
+     & ' <inertiedat>'
+        do islope=1,nslope
                 inertiesoil(:,:,islope)=inertiedat(:,:)
-       enddo
-        enddo
+        enddo
       else ! <inertiesoil> found
        if (interpol) then ! put values in oldtsoil