Changeset 1347 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phys_output_write.h

Timestamp:

Apr 13, 2010, 5:12:56 PM (14 years ago)

Author:

Laurent Fairhead

Message:

Additions to aerosol outputs for CMIP5 exercise
(Needed because of chageset r1346 LF)

---

Additions aux sorties aérosols pour l'exercice CMIP5
(Nécessaires suite au changeset r1346 LF)

Michael, Anne

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phys_output_write.h (modified) (9 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phys_output_write.h

@@ -146 +146 @@
        ENDIF
 
+       IF (o_msnow%flag(iff)<=lev_files(iff)) THEN
+      CALL histwrite_phy(nid_files(iff),o_msnow%name,itau_w,snow_o)
+       ENDIF
+
+       IF (o_fsnow%flag(iff)<=lev_files(iff)) THEN
+      CALL histwrite_phy(nid_files(iff),o_fsnow%name,itau_w,zfra_o)
+       ENDIF
+
        IF (o_evap%flag(iff)<=lev_files(iff)) THEN
       CALL histwrite_phy(nid_files(iff),o_evap%name,itau_w,evap)
@@ -190 +198 @@
        ENDIF
 
+       IF (o_nettop%flag(iff)<=lev_files(iff)) THEN
+      zx_tmp_fi2d(:) = topsw(:)-toplw(:)
+      CALL histwrite_phy(nid_files(iff),
+     $                  o_nettop%name,itau_w,zx_tmp_fi2d)
+       ENDIF
+
        IF (o_SWup200%flag(iff)<=lev_files(iff)) THEN
       CALL histwrite_phy(nid_files(iff),o_SWup200%name,itau_w,SWup200)
@@ -319 +333 @@
        CALL histwrite_phy(nid_files(iff),
      $                   o_fqfonte%name,itau_w,zxfqfonte)
+       ENDIF
+
+       IF (o_taux%flag(iff)<=lev_files(iff)) THEN
+         zx_tmp_fi2d=0.
+         do nsrf=1,nbsrf
+          zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxu(:,1,nsrf)
+         enddo
+         CALL histwrite_phy(nid_files(iff),
+     $                   o_taux%name,itau_w,zx_tmp_fi2d)
+       ENDIF
+
+       IF (o_tauy%flag(iff)<=lev_files(iff)) THEN
+         zx_tmp_fi2d=0.
+         do nsrf=1,nbsrf
+          zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxv(:,1,nsrf)
+         enddo
+         CALL histwrite_phy(nid_files(iff),
+     $                   o_tauy%name,itau_w,zx_tmp_fi2d)
        ENDIF
 
@@ -377 +409 @@
       ENDIF
 
+      IF (o_evap_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
+      zx_tmp_fi2d(1 : klon) = fevap(1 : klon, nsrf)
+      CALL histwrite_phy(nid_files(iff),o_evap_srf(nsrf)%name,
+     $           itau_w,zx_tmp_fi2d)
+      ENDIF
+
        IF (o_sens_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
@@ -532 +570 @@
        ENDIF
 
+       IF (o_ftime_con%flag(iff)<=lev_files(iff)) THEN
+        zx_tmp_fi2d=float(itau_con)/float(itap)
+      CALL histwrite_phy(nid_files(iff),o_ftime_con%name,
+     s                   itau_w,zx_tmp_fi2d)
+       ENDIF
+
+       IF (o_mc%flag(iff)<=lev_files(iff)) THEN
+        if(iflag_thermals.gt.1)then
+         zx_tmp_fi3d=dnwd+dnwd0+upwd+fm_therm
+        else
+         zx_tmp_fi3d=dnwd+dnwd0+upwd
+        endif 
+      CALL histwrite_phy(nid_files(iff),o_mc%name,itau_w,zx_tmp_fi3d)
+       ENDIF
+      
       ENDIF !iflag_con .GE. 3
 
@@ -780 +833 @@
       ENDIF
 
+      IF (o_rh2m_min%flag(iff)<=lev_files(iff)) THEN
+      DO i=1, klon
+       zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
+      ENDDO
+      CALL histwrite_phy(nid_files(iff),o_rh2m_min%name,
+     s               itau_w,zx_tmp_fi2d)
+      ENDIF
+
+      IF (o_rh2m_max%flag(iff)<=lev_files(iff)) THEN
+      DO i=1, klon
+       zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
+      ENDDO
+      CALL histwrite_phy(nid_files(iff),o_rh2m_max%name,
+     s              itau_w,zx_tmp_fi2d)
+      ENDIF
+
+
       IF (o_qsat2m%flag(iff)<=lev_files(iff)) THEN
       CALL histwrite_phy(nid_files(iff),o_qsat2m%name,itau_w,qsat2m)
@@ -1127 +1197 @@
         IF (o_pres%flag(iff)<=lev_files(iff)) THEN
       CALL histwrite_phy(nid_files(iff),o_pres%name,itau_w,pplay)
+        ENDIF
+
+        IF (o_paprs%flag(iff)<=lev_files(iff)) THEN
+      CALL histwrite_phy(nid_files(iff),o_paprs%name,
+     s                    itau_w,paprs(:,1:klev))
         ENDIF
 
@@ -1307 +1382 @@
 
        IF (iflag_thermals.gt.1) THEN
+        IF (o_ftime_th%flag(iff)<=lev_files(iff)) THEN
+! Pour l instant 0 a y reflichir pour les thermiques
+         zx_tmp_fi2d=0. 
+        CALL histwrite_phy(nid_files(iff),o_ftime_th%name,
+     s                     itau_w,zx_tmp_fi2d)
+        ENDIF
+
         IF (o_f_th%flag(iff)<=lev_files(iff)) THEN
         CALL histwrite_phy(nid_files(iff),o_f_th%name,itau_w,fm_therm)
@@ -1336 +1418 @@
        ENDIF
 
-       ENDIF !iflag_thermals
+      ENDIF !iflag_thermals
 
        IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN