Changeset 2867 for trunk/LMDZ.GENERIC/libf

Timestamp:

Jan 16, 2023, 4:47:07 PM (22 months ago)

Author:

jleconte

Message:

Small fixes to iniaerosol and physiq_mod

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• iniaerosol.F (modified) (1 diff)
• physiq_mod.F90 (modified) (5 diffs)

trunk/LMDZ.GENERIC/libf/phystd/iniaerosol.F

@@ -35 +35 @@
      & allocate(i_rgcs_ice(aerogeneric)) 
       ! Init of i_rgcs_ice
-      i_rgcs_ice(:) =0.0
+      i_rgcs_ice(:) =0
       ia = 1
       do iq=1,nqtot

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -478 +478 @@
 ! --------------------------------
       if (firstcall) then
-        call getin_p("check_physics_inputs", check_physics_inputs)
-        call getin_p("check_physics_outputs", check_physics_outputs)
-
-        ! Allocate saved arrays (except for 1D model, where this has already
-        ! been done)
+         call getin_p("check_physics_inputs", check_physics_inputs)
+         call getin_p("check_physics_outputs", check_physics_outputs)
+
+         ! Allocate saved arrays (except for 1D model, where this has already
+         ! been done)
 #ifndef MESOSCALE
-        if (ngrid>1) call phys_state_var_init(nq)
+         if (ngrid>1) call phys_state_var_init(nq)
 #endif
 
@@ -511 +511 @@
 
 #ifdef CPP_XIOS
-        ! Initialize XIOS context
-        write(*,*) "physiq: call wxios_context_init"
-        CALL wxios_context_init
+         ! Initialize XIOS context
+         write(*,*) "physiq: call wxios_context_init"
+         CALL wxios_context_init
 #endif
 
@@ -525 +525 @@
                        rnat,pctsrf_sic,tslab, tsea_ice,sea_ice)
 #else
-      emis(:)=0.0
-      q2(:,:)=0.0
-      qsurf(:,:)=0.0
-      day_ini = pday
+         day_ini = pday
 #endif
 
@@ -557 +554 @@
 !        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
          albedo(:,:)=0.0
-          albedo_bareground(:)=0.0
-          albedo_snow_SPECTV(:)=0.0
-          albedo_co2_ice_SPECTV(:)=0.0
+         albedo_bareground(:)=0.0
+         albedo_snow_SPECTV(:)=0.0
+         albedo_co2_ice_SPECTV(:)=0.0
          call surfini(ngrid,nq,qsurf,albedo,albedo_bareground,albedo_snow_SPECTV,albedo_co2_ice_SPECTV) 
          
@@ -2596 +2593 @@
    comm_HR_DYN(1:ngrid,1:nlayer) = zdtdyn(1:ngrid,1:nlayer)
 
-      if (turb_resolved) then
-        open(17,file='lsf.txt',form='formatted',status='old')
-        rewind(17)
-        DO l=1,nlayer
-          read(17,*) lsf_dt(l),lsf_dq(l)
-        ENDDO
-        close(17)
-        do ig=1,ngrid
-          if ((tracer).and.(water)) then
-           pdq(ig,:,igcm_h2o_vap) = pdq(ig,:,igcm_h2o_vap) + lsf_dq(:)
-          endif
-          pdt(ig,:) = pdt(ig,:) + lsf_dt(:)
-          comm_HR_DYN(ig,:) = lsf_dt(:)
-        enddo
-      endif
+!      if (turb_resolved) then
+!        open(17,file='lsf.txt',form='formatted',status='old')
+!        rewind(17)
+!        DO l=1,nlayer
+!          read(17,*) lsf_dt(l),lsf_dq(l)
+!        ENDDO
+!        close(17)
+!        do ig=1,ngrid
+!          if ((tracer).and.(water)) then
+!           pdq(ig,:,igcm_h2o_vap) = pdq(ig,:,igcm_h2o_vap) + lsf_dq(:)
+!          endif
+!          pdt(ig,:) = pdt(ig,:) + lsf_dt(:)
+!          comm_HR_DYN(ig,:) = lsf_dt(:)
+!        enddo
+!      endif
 #endif