Changeset 317 for codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf

Timestamp:

01/14/15 11:10:01 (9 years ago)

Author:

millour

Message:

Add possibility to output Ls (as a surface field). EM.

File:

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

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

@@ -236 +236 @@
 
       real zls                       ! solar longitude (rad)
+      real zls_ngrid(ngrid) ! dummy array to output Ls (as a surface field)
       real zday                      ! date (time since Ls=0, in martian days)
       real zzlay(ngrid,nlayer)   ! altitude at the middle of the layers
@@ -2238 +2239 @@
       icount=icount+1
 
-!!!!!!!!!!!!!!!! section for XIOS output !!!!!!!!!!!!!!!      
+!!!!!!!!!!!!!!!! section for XIOS output !!!!!!!!!!!!!!!    
+      zls_ngrid(1:ngrid)= zls*180./pi
+      call write_xios_field("Ls",zls_ngrid)
       CALL write_xios_field("tsurf",tsurf)
       CALL write_xios_field("ps",ps)