Changeset 270 for codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/physiq.F90

Timestamp:

09/05/14 14:42:27 (10 years ago)

Author:

millour

Message:

Added features for the Saturn case:

• Added possibility to run without startfi or restartfi.nc files
• Added reference temperature "temp_profile.txt" profile to start from
• More XIOS outputs, and put them on "presnivs (pressure) vertical coordinate
• Added "-openmp-threadprivate compat" OpenMP option in Ada arch file

EM

File:

• codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/physiq.F90 (modified) (6 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/physiq.F90

@@ -569 +569 @@
 !ym               rnat,pctsrf_sic,tslab, tsea_ice,sea_ice)
 
+         write(*,*) "physiq: firstcall, call phyetat0_academic"
          call phyetat0_academic(ngrid,nlayer,"startfi.nc",0,0,nsoilmx,nq,   &
                day_ini,time_phys,tsurf,tsoil,emis,q2,qsurf,   &
@@ -575 +576 @@
 
 !mi initialising tsurf with pt
+         write(*,*) "Physiq: initializing tsurf(:) to pt(:,1) !!"
          tsurf(:)=pt(:,1)
 
@@ -744 +746 @@
 
          if (ngrid.ne.1) then ! no need to create a restart file in 1d
-           call physdem0("restartfi.nc",long,lati,nsoilmx,ngrid,nlayer,nq, &
-                         ptimestep,pday+nday,time_phys,area, &
-                         albedodat,inertiedat,zmea,zstd,zsig,zgam,zthe)
+! EM: No restart file (for now).
+!           call physdem0("restartfi.nc",long,lati,nsoilmx,ngrid,nlayer,nq, &
+!                         ptimestep,pday+nday,time_phys,area, &
+!                         albedodat,inertiedat,zmea,zstd,zsig,zgam,zthe)
          endif
          
@@ -1972 +1975 @@
 !                    cloudfrac,totcloudfrac,hice,noms)
 !#endif
-              call physdem1("restartfi.nc",nsoilmx,ngrid,nlayer,nq, &
-                      ptimestep,ztime_fin, &
-                      tsurf,tsoil,emis,q2,qsurf_hist, &
-                      cloudfrac,totcloudfrac,hice, &
-                      rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)
+
+! EM: do not write a restart file (for now).
+!              call physdem1("restartfi.nc",nsoilmx,ngrid,nlayer,nq, &
+!                      ptimestep,ztime_fin, &
+!                      tsurf,tsoil,emis,q2,qsurf_hist, &
+!                      cloudfrac,totcloudfrac,hice, &
+!                      rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)
             endif
 
@@ -2239 +2244 @@
       CALL write_xios_field("ASR",fluxabs_sw)
       CALL write_xios_field("OLR",fluxtop_lw)
+      call write_xios_field("input_temp",pt)
+      call write_xios_field("input_u",pu)
+      call write_xios_field("input_v",pv)
+      call write_xios_field("dtrad",dtrad)
+      call write_xios_field("zdtlw",zdtlw)
+      call write_xios_field("zdtsw",zdtsw)
+      call write_xios_field("zdtdyn",zdtdyn/ptimestep)
+      call write_xios_field("zdtdif",zdtdif)
+      call write_xios_field("zdtadj",zdtadj)
+      call write_xios_field("pdt",pdt)
       IF (lastcall) CALL finalize_xios_output
       
@@ -2344 +2359 @@
       endif
 
+      write(*,*) "physiq: done, zday=",zday
       return
     end subroutine physiq