Changeset 3491 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Timestamp:

Nov 4, 2024, 4:15:02 PM (7 weeks ago)

Author:

emillour

Message:

Mars PCM:
Get rid of "start2archive_SSO.F" and adapt start2archive.F to allow adding
sub-grid-scale fields in start_archive.nc. This optional behavior is
triggered at run time by specifying "start2archive.e --add-sso".
EM

Location:

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Files:

• start2archive.F (modified) (5 diffs)
• start2archive_SSO.F (deleted)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F

@@ -38 +38 @@
       USE comcstfi_h, only: pi
       use surfini_mod, only: surfini
+! SSO parameters:
+      USE surfdat_h, ONLY: phisfi, albedodat, z0, z0_default,
+     &    zmea, zstd, zsig, zgam, zthe, hmons, summit, base
       implicit none
 
@@ -100 +103 @@
       REAL, ALLOCATABLE :: subslope_distS(:,:)
 
+! For SSO parameters:
+      REAL zmeaS(ip1jmp1)
+      REAL zsigS(ip1jmp1)
+      REAL zstdS(ip1jmp1)
+      REAL zgamS(ip1jmp1)
+      REAL ztheS(ip1jmp1)
+      REAL albedodatS(ip1jmp1)
+      REAL z0S(ip1jmp1)
+      REAL hmonsS(ip1jmp1)
+      REAL summitS(ip1jmp1)
+      REAL baseS(ip1jmp1)
+
+      logical :: add_sso_fields=.false. ! default, don't include SSO fields
+
+
 c Variables intermediaires : vent naturel, mais pas coord scalaire
 c----------------------------------------------------------------
@@ -128 +146 @@
       integer varid,dimid,timelen 
       INTEGER nid,nid1
+
+C get command line arguments
+C here we assume and check that if there is an argument #1 then
+C it should be --add-sso to signal adding SSO fileds to start_archive.nc
+
+      CALL get_command_argument(1,txt,j,ierr)
+      ! will return ierr==0 if there is an argument #1 to command line
+      IF (ierr==0) THEN
+        ! Check that argument is indeed "--add-sso" or signal the error
+        IF (trim(txt)=="--add-sso") THEN
+          add_sso_fields=.true.
+          write(*,*) "SSO fields will be included in start_archive"
+        ELSE
+          write(*,*) "start2archive error: unexpected command line "//
+     &    "argument: ",trim(txt)
+          write(*,*) " (only --add-sso currently accepted)"
+          write(*,*) "Might as well stop here."
+          stop
+        ENDIF
+      ENDIF ! of IF (ierr==0)
 
 c-----------------------------------------------------------------------
@@ -326 +364 @@
      &    watercaptagS(:))
 
+! SSO parameters, if needed:
+      if (add_sso_fields) then
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,zmea,zmeaS)
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,zstd,zstdS)
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,zsig,zsigS)
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,zthe,ztheS)
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,zgam,zgamS)
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,albedodat,albedodatS)
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,z0,z0S)
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,hmons,hmonsS)
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,summit,summitS)
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,base,baseS)
+      endif
+
 c=======================================================================
 c Info pour controler
@@ -445 +497 @@
       call write_archive(nid,ntime,'q2atm','wind variance','m2.s-2',3,
      .              q2S(1,2))
+! SSO parameters, if needed
+      if (add_sso_fields) then
+        call write_archive(nid,ntime,'ZMEA','zmea',' ',2,zmeaS)
+        call write_archive(nid,ntime,'ZSTD','zstd',' ',2,zstdS)
+        call write_archive(nid,ntime,'ZSIG','zsig',' ',2,zsigS)
+        call write_archive(nid,ntime,'ZTHE','zthe',' ',2,ztheS)
+        call write_archive(nid,ntime,'ZGAM','zgam',' ',2,zgamS)
+        call write_archive(nid,ntime,'albedodat','albedodat',
+     &                             ' ',2,albedodatS)
+        call write_archive(nid,ntime,'z0','z0',' ',2,z0S)
+        call write_archive(nid,ntime,'summit','summit',
+     &                             ' ',2,summitS)
+        call write_archive(nid,ntime,'hmons','hmons',' ',2,hmonsS)
+        call write_archive(nid,ntime,'base','base',' ',2,baseS)
+      endif
 
 c-----------------------------------------------------------------------