Changeset 3756

Timestamp:

May 7, 2025, 5:49:03 PM (4 weeks ago)

Author:

emillour

Message:

Mars PCM:
Code tidying: turn aerave.F and suaer.F90 into modules and modernize
"blackl" routine (enforce "implicit none", make true constants "parameters")
and include it in the aerave module since it is only called there.
EM

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/aerave.F (modified) (2 diffs)
• libf/phymars/blackl.F (deleted)
• libf/phymars/callradite_mod.F (modified) (1 diff)
• libf/phymars/suaer.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4826 +4826 @@
 == 28/04/2025 == EM
 OpenMP bug fix: add some missing threadprivate statements in moldiff_MPF
+
+== 07/05/2025 == EM
+Code tidying: turn aerave.F and suaer.F90 into modules and modernize
+"blackl" routine (enforce "implicit none", make true constants "parameters")
+and include it in the aerave module since it is only called there. 

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

@@ +1 @@
+      MODULE aerave_mod
+      
+      IMPLICIT NONE
+      
+      CONTAINS
+
       SUBROUTINE aerave ( ndata,
      & longdata,epdata,omegdata,gdata,          
@@ -191 +197 @@
 c......................................................................
 c
-      RETURN
-      END
+      END SUBROUTINE aerave
+
+c.......................................................................
+
+      subroutine blackl(blalong,blat,blae)
+
+      implicit none
+c.......................................................................
+      double precision,intent(in) :: blalong,blat
+      double precision,intent(out) :: blae
+c.....parameters
+      double precision,parameter :: sigma=5.6693d-08
+      double precision,parameter :: pi=datan(1.d0)*4.d0
+      double precision,parameter :: c0=2.9979d+08
+      double precision,parameter :: h=6.6262d-34
+      double precision,parameter :: cbol=1.3806d-23
+      double precision,parameter :: rind=1.d0
+      double precision,parameter :: c=c0/rind
+      double precision,parameter :: c1=h*(c**2)
+      double precision,parameter :: c2=h*c/cbol
+c.......................................................................
+      blae=2.d0*pi*c1/blalong**5/(exp(c2/blalong/blat)-1.d0)
+c.......................................................................
+
+      end subroutine blackl
+
+      END MODULE aerave_mod

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

@@ -27 +27 @@
       use lwmain_mod, only: lwmain
       use swmain_mod, only: swmain
+      use suaer_mod, only: suaer
       use dust_param_mod, only: doubleq, active, submicron
       use callkeys_mod, only: water, activice, rdstorm, topflows,

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

@@ +1 @@
+MODULE suaer_mod
+
+IMPLICIT NONE
+
+CONTAINS
+
 SUBROUTINE suaer
 use dimradmars_mod, only: longrefvis, longrefir, nsizemax, long1vis, &
@@ -15 +21 @@
 USE mod_phys_lmdz_transfert_para, ONLY: bcast
 use callkeys_mod, only: dustrefir
+use aerave_mod, only: aerave
 
 IMPLICIT NONE
@@ -568 +575 @@
 
 !==================================================================
-RETURN
-END
+
+END SUBROUTINE suaer
+
+END MODULE suaer_mod