Changeset 2912

Timestamp:

Mar 12, 2023, 6:15:33 PM (22 months ago)

Author:

emillour

Message:

Dynamico interface with the PCM physics:

• update disvert_icosa_lmdz.f90 to follow the update of disvert_plugin in dynamico
• enforce using fcm1 when compiling dynamico via the make_icosa_lmdz script
• add an illustrative example script to compile on Irene-Rome

EM

Location:

trunk/ICOSA_LMDZ

Files:

• compile_irene-amd (added)
• make_icosa_lmdz (modified) (1 diff)
• src/disvert_icosa_lmdz.f90 (modified) (3 diffs)

trunk/ICOSA_LMDZ/make_icosa_lmdz

@@ -227 +227 @@
 if [[ -d ../ICOSAGCM ]] ; then
     cd ../ICOSAGCM 
-    ./make_icosa -$compil_mode -parallel $parallel -external_ioipsl -with_xios -arch $arch -arch_path $arch_path -job $job $full_flag || exit 1
+    # for now compile using fcm1
+    ./make_icosa -with_fcm1 -$compil_mode -parallel $parallel -external_ioipsl -with_xios -arch $arch -arch_path $arch_path -job $job $full_flag || exit 1
     cd -
 else

trunk/ICOSA_LMDZ/src/disvert_icosa_lmdz.f90

@@ -5 +5 @@
 contains
 
-  SUBROUTINE disvert_icosa_lmdz(ap,bp,presnivs)
+  SUBROUTINE disvert_icosa_lmdz(ap,bp,presnivs,presinter)
   ! Routine to compute ap(),bp() and presnivs() in the same way as done
   ! in LMDZ for planets (see LMDZ.COMMON/libf/dyn3d_common/disvert_noterre.F)
@@ -23 +23 @@
   REAL(rstd),INTENT(OUT) :: bp(:)
   REAL(rstd),INTENT(OUT) :: presnivs(:)
+  REAL(rstd),INTENT(OUT) :: presinter(:)
 
   INTEGER :: unit
@@ -111 +112 @@
       presnivs(l)=0.5*(ap(l)+bp(l)*preff+ap(l+1)+bp(l+1)*preff)
     ENDDO
+    ! build presinter(), approximative inter-layer pressures
+    DO l=1,llm+1
+      presinter(l)=ap(l)+bp(l)*preff
+    ENDDO
     
     ! tell the world about it