Changeset 2868 for trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd

Timestamp:

Jan 16, 2023, 4:47:08 PM (3 years ago)

Author:

jleconte

Message:

Big cleanup of interface_v4. M_TSURF is now P_TSURF.

Location:

trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd

Files:

• callphysiq_mod.F (modified) (1 diff)
• update_inputs_physiq_mod.F (modified) (14 diffs)
• update_outputs_physiq_mod.F (modified) (2 diffs)

trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd/callphysiq_mod.F

@@ -43 +43 @@
 
   call allocate_comm_wrf(klon,llm)
-
 ! Call physics package with required inputs/outputs
   CALL physiq(klon,           & ! ngrid
               llm,            & ! nlayer
               nqtot,          & ! nq
-!              noms,          & ! nametrac
               debut_split,    & ! firstcall
               lafin_split,    & ! lastcall

trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd/update_inputs_physiq_mod.F

@@ -27 +27 @@
   REAL :: sec,nsec
   
-  IF (JULYR .le. 8999) THEN
-    if (tlocked .eqv. .false.) THEN
+  !IF (JULYR .le. 8999) THEN
+  !  if (tlocked .eqv. .false.) THEN
       JH_cur_split = (GMT + elaps/3600.) !! universal time (0<JH_cur_split<1): JH_cur_split=0.5 at 12:00 UT 
       JH_cur_split = MODULO(JH_cur_split,24.)   !! the two arguments of MODULO must be of the same type 
@@ -36 +36 @@
       MY = (JULYR-2000) + (86400*(JULDAY - 1)+3600.0*GMT+elaps)/31968000
       MY = INT(MY)
-    ELSE
-      JH_cur_split = (GMT)! + elaps/420000.) !! universal time (0<JH_cur_split<1): JH_cur_split=0.5 at 12:00 UT 
-      JH_cur_split = MODULO(JH_cur_split,24.)   !! the two arguments of MODULO must be of the same type 
-      JH_cur_split = JH_cur_split / 24.
-      JD_cur = (JULDAY - 1 + INT((3600*GMT)/86400))!+elaps)/1.008e7))
-      JD_cur = MODULO(int(JD_cur),2)
-      MY = (JULYR-2000) + (86400*(JULDAY - 1)+3600*GMT+elaps)/31968000
-      MY = INT(MY)
-    ENDIF
-  ELSE
-    if (tlocked .eqv. .false.) THEN
-      JH_cur_split = lct_input - lon_input / 15. + elaps/1500.0
-      JD_cur =  INT((sec*(lct_input - lon_input / 15.) + elaps)/36000)
-      !ptime = lct_input - lon_input / 15. + elaps/3600.
-      !pday =  INT((3600*(lct_input - lon_input / 15.) + elaps)/86400)
-    ELSE
-      JH_cur_split = lct_input - lon_input / 15. !+ elaps/1500.0
-      !pday =  INT((sec*(lct_input - lon_input / 15.)+ elaps)/36000)
-      JD_cur =  INT((sec*(lct_input - lon_input / 15.))/3600)
-      !print*,'ptime',ptime
-      !print*,'pday',pday 
-    !pday =  INT((3600*(lct_input - lon_input / 15.) + elaps)/86400)
-      JH_cur_split = MODULO(ptime,24.)
-      JH_cur_split = JH_cur_split / 24.
-      JD_cur = MODULO(int(pday),365)
-      MY = 2024
-    ENDIF
-  ENDIF
+  !  ELSE
+  !    JH_cur_split = (GMT)! + elaps/420000.) !! universal time (0<JH_cur_split<1): JH_cur_split=0.5 at 12:00 UT 
+  !    JH_cur_split = MODULO(JH_cur_split,24.)   !! the two arguments of MODULO must be of the same type 
+  !    JH_cur_split = JH_cur_split / 24.
+  !    JD_cur = (JULDAY - 1 + INT((3600*GMT)/86400))!+elaps)/1.008e7))
+  !    JD_cur = MODULO(int(JD_cur),2)
+  !    MY = (JULYR-2000) + (86400*(JULDAY - 1)+3600*GMT+elaps)/31968000
+  !    MY = INT(MY)
+  !  ENDIF
+  !ELSE
+  !  if (tlocked .eqv. .false.) THEN
+  !    JH_cur_split = lct_input - lon_input / 15. + elaps/1500.0
+  !    JD_cur =  INT((sec*(lct_input - lon_input / 15.) + elaps)/36000)
+  !    !ptime = lct_input - lon_input / 15. + elaps/3600.
+  !    !pday =  INT((3600*(lct_input - lon_input / 15.) + elaps)/86400)
+  !  ELSE
+  !    JH_cur_split = lct_input - lon_input / 15. !+ elaps/1500.0
+  !    !pday =  INT((sec*(lct_input - lon_input / 15.)+ elaps)/36000)
+  !    JD_cur =  INT((sec*(lct_input - lon_input / 15.))/3600)
+  !    !print*,'ptime',ptime
+  !    !print*,'pday',pday 
+  !  !pday =  INT((3600*(lct_input - lon_input / 15.) + elaps)/86400)
+  !    JH_cur_split = MODULO(ptime,24.)
+  !    JH_cur_split = JH_cur_split / 24.
+  !    JD_cur = MODULO(int(pday),365)
+  !    MY = 2024
+  !  ENDIF
+  !ENDIF
+
 
 END SUBROUTINE update_inputs_physiq_time
@@ -80 +81 @@
       
                  !! tableau dans tracer_mod.F90
-  nqtot=nq
   if ((TRACER_MODE .eq. 1) .or. (TRACER_MODE .ge. 42)) THEN
     nqtot=2
@@ -107 +107 @@
 SUBROUTINE update_inputs_physiq_constants
 
-   !USE module_model_constants
-   !use comcstfi_h, only: omeg,mugaz
-   !use planete_h,  only: year_day,periheli,aphelie, &
-   !                       peri_day,obliquit,emin_turb, &
-   !                       lmixmin
    use planete_mod, only: year_day, periastr, apoastr, peri_day,&
                        obliquit, z0, lmixmin, emin_turb
@@ -117 +112 @@
                          emisice,dtemisice
    !                       z0_default
-   !use comsoil_h, only: volcapa
    use comcstfi_mod, only: rad, omeg, g, mugaz, rcp, cpp, r
-   !! comcstfi_h
    use phys_state_var_mod, only :cloudfrac,totcloudfrac,hice,rnat,pctsrf_sic,tsea_ice
-   !use iniorbit
-   !use time_phylmdz_mod, only: dtphys, daysec,day_ini
-  
-  
-     !open(17,file='controle.txt',form='formatted',status='old')
-     !rewind(17)
-     !read(17,*)
-     !read(17,*) 
-     !read(17,*) day_ini !(tab0+3)
-     !read(17,*) 
-     !read(17,*) !tab0+5)
-     !read(17,*) omeg !(tab0+6)
-     !read(17,*) !(tab0+7)
-     !read(17,*)  !(tab0+8)
-     !read(17,*)  !(tab0+9)
-     !read(17,*) daysec
-     !read(17,*) dtphys !tab0+11)
-     !read(17,*)
-     !read(17,*)
-     !read(17,*) year_day !(tab0+14)
-     !read(17,*) periastr !tab0+15)
-     !read(17,*) apoastr !tab0+16)
-     !read(17,*) peri_day !tab0+17)
-     !read(17,*) obliquit !tab0+18)
-     !read(17,*) z0
-     !read(17,*)
-     !read(17,*)
-     !read(17,*)
-     !read(17,*)
-     !read(17,*) emisice(1)
-     !read(17,*) emisice(2)
-     !read(17,*) emissiv
-     !read(17,*)
-     !read(17,*)
-     !read(17,*)
-     !read(17,*)
-     !read(17,*) iceradius(1)
-     !read(17,*) iceradius(2)
-     !read(17,*) dtemisice(1)
-     !read(17,*) dtemisice(2)
-     !close(17) 
-     !cpp=(8.314511/(mugaz/1000.0))/rcp
-     !print*,'cpp',cpp
-     !print*,'g',g
-
-     !emissiv(:)=EMIS
-     !cloudfrac(:,:)=0.5
-     !totcloudfrac(:)=0.5
-     !hice(:)=0.
-     !rnat(:)=0.
-     !pctsrf_sic(:)=0.
-     !tsea_ice(:)=0.
-     !qsurf(:,:) = 0. 
-     !print*,'iceradius',iceradius,'dtemisice',dtemisice
-     !print*,'apoastr,periastr,year_day,peri_day,obliq',apoastr,periastr,year_day,peri_day,obliquit
-     !print*,'emissiv',emissiv
- 
-     ! SUBROUTINE iniorbit(apoastr,periastr,year_day,peri_day,obliq)
-
+   !JL22 this routine does not do anything for the generic interface
+   ! The various use abave can surely be removed.
 END SUBROUTINE update_inputs_physiq_constants
 
@@ -272 +208 @@
    latitude_deg(:) = plat(:)/DEGRAD
    cell_area(:) = parea(:)
-   !call planetwide_sumval(parea,totarea_planet)
-   !print*,'parea',parea(1)
-   !totarea=SSUM(ngrid,parea,1)
    totarea=ngrid*parea(1)
-   !totarea_planet=SSUM(ngrid,parea,1)
    totarea_planet=ngrid*parea(1)
 
@@ -294 +226 @@
    DO k=1, nlayer
    read(12,*) znw(k)
-   !write(6,*) 'read level ', k,grid%znw(k)
    ENDDO
    close(12)
@@ -319 +250 @@
             JULYR,TRACER_MODE,&
             M_ALBEDO,CST_AL,&
-            M_TSURF,M_EMISS,M_CO2ICE,&
+            P_TSURF,M_EMISS,M_CO2ICE,&
             M_GW,M_Z0,CST_Z0,&
             M_H2OICE,&
@@ -335 +266 @@
    REAL, INTENT(IN  ) :: CST_AL, phisfi_val, CST_Z0
    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN)  :: &
-     M_ALBEDO,M_TSURF,M_EMISS,M_CO2ICE,M_H2OICE,M_Z0
+     M_ALBEDO,P_TSURF,M_EMISS,M_CO2ICE,M_H2OICE,M_Z0
    REAL, DIMENSION( ims:ime, 5, jms:jme ), INTENT(IN   )  :: M_GW  
 
-   !print*,'ALLOCATED(phisfi)',ALLOCATED(phisfi)
-   !print*,'size phisfi',size(phisfi)
    DO j = jps,jpe
    DO i = ips,ipe
@@ -352 +281 @@
      !---------------------!
      phisfi(subs) = phisfi_val
-!     print*,'size phisfi',size(phisfi)
-     !print*,'phisfi',phisfi(subs)
      !---------------!
      ! Ground albedo !
@@ -383 +310 @@
      !----------------------------!
      z0 = CST_Z0
-     !IF (JULYR .le. 8999) THEN
-     ! IF (CST_Z0 == 0) THEN
-     !  z0(subs) = M_Z0(i,j)
-     ! ELSE
-     !  z0(subs) = CST_Z0
-     !  IF ( (i == ips) .AND. (j == jps) ) PRINT *,'** Mars ** SET CONSTANT SURF ROUGHNESS (m) ',CST_Z0
-     ! ENDIF
-     !ELSE
-     ! IF ( (i == ips) .AND. (j == jps) ) PRINT *,'** Mars ** IDEALIZED SIMULATION z0 (m) ', CST_Z0
-     ! z0(subs)=CST_Z0
-     !ENDIF
-     !!!!! ADDITIONAL SECURITY. THIS MIGHT HAPPEN WITH OLD INIT FILES.
-     !IF (z0(subs) == 0.) THEN
-     !  IF ( (i == ips) .AND. (j == jps) ) PRINT *, 'WELL, z0 is 0, this is no good. Setting to old defaults value 0.01 m'
-     !  z0(subs) = 0.01
-     !ENDIF
-     !!!!! ADDITIONAL SECURITY. INTERP+SMOOTH IN GEOGRID MIGHT YIELD NEGATIVE Z0 !!!
-     !IF (z0(subs) < 0.) THEN
-     !  PRINT *, 'WELL, z0 is NEGATIVE, this is impossible. better stop here.'
-     !  PRINT *, 'advice --> correct interpolation / smoothing of z0 in WPS'
-     !  PRINT *, '           -- or check the constant value set in namelist.input'
-     !  STOP
-     !ENDIF
-
+     
      !-----------------------------------------------!
      ! Ground temperature, emissivity, CO2 ice cover !
      !-----------------------------------------------!
-     tsurf(subs) = M_TSURF(i,j)
+     tsurf(subs) = P_TSURF(i,j)
      emis(subs) = M_EMISS(i,j)     
      !do i=1,noceanmx
@@ -466 +370 @@
             M_TI,CST_TI,&
             M_ISOIL,M_DSOIL,&
-            M_TSOIL,M_TSURF)
+            M_TSOIL,P_TSURF)
 
    use comsoil_h, only: inertiedat,mlayer,layer,volcapa
@@ -476 +380 @@
    REAL, INTENT(IN  ) :: CST_TI 
    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN)  :: &
-     M_TI, M_TSURF
+     M_TI, P_TSURF
    REAL, DIMENSION( ims:ime, nsoil, jms:jme ), INTENT(IN)  :: &
      M_TSOIL, M_ISOIL, M_DSOIL
@@ -541 +445 @@
        IF ( (i == ips) .AND. (j == jps) ) PRINT *,'** Mars ** no tsoil. set it to tsurf.'
        do k=1,nsoil
-        !print*,'M_TSURF(i,j)',M_TSURF(1,:)
-        !print*,'size M_TSURF',size(M_TSURF)
-        !print*,'size tsoil',size(tsoil)
-        tsoil(subs,k) = M_TSURF(i,j)
-        !print*,'tsoil(subs,k)',tsoil(subs,k)
+        tsoil(subs,k) = P_TSURF(i,j)
        enddo
      ENDIF

trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd/update_outputs_physiq_mod.F

@@ -9 +9 @@
             ips,ipe,jps,jpe,&
             TRACER_MODE,&
-            M_TSURF,M_CO2ICE,&
+            P_TSURF,M_CO2ICE,&
             M_H2OICE)
 
@@ -20 +20 @@
    INTEGER :: i,j,subs
    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)  :: &
-     M_TSURF,M_CO2ICE,M_H2OICE
-
-   DO j = jps,jpe
-   DO i = ips,ipe
-
-     !-----------------------------------!
-     ! 1D subscript for physics "cursor" !
-     !-----------------------------------!
-     subs = (j-jps)*(ipe-ips+1)+(i-ips+1)
-
-     !-------------------------------------------------------!
-     ! Save key variables for restart and output and nesting !  
-     !-------------------------------------------------------!
-     M_TSURF(i,j) = tsurf(subs)
+     P_TSURF,M_CO2ICE,M_H2OICE
+
+   DO j = jps,jpe
+   DO i = ips,ipe
+
+     !-----------------------------------!
+     ! 1D subscript for physics "cursor" !
+     !-----------------------------------!
+     subs = (j-jps)*(ipe-ips+1)+(i-ips+1)
+
+     !-------------------------------------------------------!
+     ! Save key variables for restart and output and nesting !  
+     !-------------------------------------------------------!
+     P_TSURF(i,j) = tsurf(subs)
 
    ENDDO