Changeset 3301

Timestamp:

Apr 15, 2024, 5:27:36 PM (8 months ago)

Author:

jbclement

Message:

Mars PCM:
Modification of "subslope_mola.F90" to take into account other longitude resolution than 64. It makes "newstart" be able to change the number of sub-slopes for low resolution.
JBC

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/dynphy_lonlat/phymars/subslope_mola.F90 (modified) (6 diffs)
• libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (1 diff)

trunk/LMDZ.MARS/changelog.txt

@@ -4599 +4599 @@
 wind stress and dust devils. Renamed "dustlift" routine "dust_windstress_lift"
 and made it a module; also made dustdevil a module.
+
+== 15/04/2024 == JBC
+Modification of "subslope_mola.F90" to take into account other longitude resolution than 64. It makes "newstart" be able to change the number of sub-slopes for low resolution.

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/subslope_mola.F90

@@ -19 +19 @@
 
       double precision :: resol
-      parameter(resol = 64)
+      parameter(resol = real(iim))
       integer :: jjm_mola, iim_mola 
       parameter(jjm_mola=180*resol, iim_mola=2*jjm_mola)
@@ -60 +60 @@
 !-------------Building of theta_mola and psi_mola 
       ! Assume that the mola file is to ben found in "datadir"
-      CALL mola(trim(datadir)//"/",&
+      CALL mola(trim(datadir)//"/",iim, &
                 ierr,theta_mola,psi_mola,resol,iim_mola,jjm_mola)
 
@@ -111 +111 @@
 !==========================================================================================
 
- SUBROUTINE mola(dset,ierr,theta_mola,psi_mola,resol,iim_mola,jjm_mola)
+ SUBROUTINE mola(dset,iim,ierr,theta_mola,psi_mola,resol,iim_mola,jjm_mola)
 
 
@@ -124 +124 @@
 logical output_messages
       parameter (output_messages=.true.)
-      double precision resol
-!c      parameter(resol=16) ! MOLA pixel/degree resolution
-
-
-      integer jjm_mola, iim_mola   ! # of longitude and latitude MOLA data values 
+      double precision resol ! MOLA pixel/degree resolution
+      integer jjm_mola, iim_mola ! # of longitude and latitude MOLA data values 
 
 ! Arguments
 ! inputs
       character*(*) dset ! Path to MCD datafiles 
+      integer :: iim
 
 ! outputs
@@ -143 +141 @@
 
 ! Local variables
-
       real latitude  ! north latitude (degrees)m
       real longitude ! east longitude (degrees)
-      character*140 molafile ! MOLA datafile 
-!      data molafile/'mola_32.2.nc'/
-!c      data molafile/'mola16.nc'/
-!c      data molafile/'mola32.nc'/
+      character*140 molafile ! MOLA datafile
 !c      real invresol
 !c      parameter(invresol=1./resol)
@@ -168 +162 @@
       double precision lat,lon ! longitude and latitude, local values (in degrees)
       real topogrid(iim_mola,jjm_mola) ! altmola in 'real' version
+      character(2) :: iim_char
 
 
 !C 1.1. Open MOLA file
-         molafile=dset//'mola64.nc'
+         write(iim_char,'(i2)') iim
+         molafile = dset//'mola'//iim_char//'.nc'
          if (output_messages) then
            write(*,*)"Loading MOLA topography from file ", &

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -129 +129 @@
 omeg = 4.*asin(1.)/daysec ! rotation rate (rad.s-1)
 g = 3.72                  ! gravity (m.s-2) ~3.72
-mugaz = 43.49             ! atmosphere mola mass (g.mol-1) ~43.49
+mugaz = 43.49             ! atmosphere molar mass (g.mol-1) ~43.49
 rcp = .256793             ! = r/cp ~0.256793
 r = 8.314511*1000./mugaz