Changeset 3814 for dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/write_histday_seri.h

Timestamp:

Apr 15, 2015, 6:49:07 PM (10 years ago)

Author:

ymipsl

Message:

remove all dynamic dependency in LMDZ physics except for the include "dimensions.h"

YM

File:

• dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/write_histday_seri.h (modified) (2 diffs)

dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/write_histday_seri.h

@@ -51 +51 @@
       CALL histwrite(nid_day_seri,"ecin",itau_w, &
                      zx_tmp_2d,iim*jjmp1,ndex2d) 
-!
-!IM 151004 BEG
-      IF(1.EQ.0) THEN
-!
-      DO k=1, klev
-      DO i=1, klon
-       zx_tmp_fi3d(i,k)=u_seri(i,k)*RA*cos(pir* rlat(i))
-      ENDDO
-      ENDDO
-!
-      CALL moyglo_pondaima(klon, klev, zx_tmp_fi3d,  &
-           airephy, paprs, moyglo)
-      zx_tmp_fi2d(1:klon)=moyglo
-!
-      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
-      CALL histwrite(nid_day_seri,"momang",itau_w,zx_tmp_2d, &
-                     iim*jjmp1,ndex2d)
-!
-! friction torque
-!
-      DO i=1, klon
-       zx_tmp_fi2d(i)=zxfluxu(i,1)*RA* cos(pir* rlat(i))
-      ENDDO
-!
-      ok_msk=.FALSE.
-      CALL moyglo_pondaire(klon, zx_tmp_fi2d, airephy,  &
-           ok_msk, msk, moyglo)
-      zx_tmp_fi2d(1:klon)=moyglo
-!
-      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
-      CALL histwrite(nid_day_seri,"frictor",itau_w,zx_tmp_2d, &
-                     iim*jjmp1,ndex2d)
-!
-! mountain torque
-!
-!IM 190504 BEG
-      CALL gr_fi_dyn(1,klon,iim+1,jjm+1,airephy,airedyn)
-      CALL gr_fi_dyn(klev+1,klon,iim+1,jjm+1,paprs,padyn)
-      CALL gr_fi_dyn(1,klon,iim+1,jjm+1,rlat,rlatdyn)
-      mountor=0.
-      airetot=0.
-      DO j = 1, jjmp1
-       DO i = 1, iim+1
-        ij=i+(iim+1)*(j-1)
-        zx_tmp(ij)=0.
-        DO k = 1, klev
-         zx_tmp(ij)=zx_tmp(ij)+dudyn(i,j,k)*airedyn(i,j)* &
-                    (padyn(i,j,k+1)-padyn(i,j,k))/RG
-         airetot=airetot+airedyn(i,j)
-        ENDDO
-!IM 190504 mountor=mountor+zx_tmp(ij)*airedyn(i,j)*RA*
-        mountor=mountor+zx_tmp(ij)*RA* &
-                 cos(pir* rlatdyn(i,j))
-       ENDDO
-      ENDDO
-!IM 151004 BEG
-      IF(itap.EQ.1) PRINT*,'airetot=',airetot,airetot/klev
-!IM 151004 END
-!IM 190504      mountor=mountor/(airetot*airetot)
-      mountor=mountor/airetot
-!
-!IM 190504 END
-      zx_tmp_2d(1:iim,1:jjmp1)=mountor
-      CALL histwrite(nid_day_seri,"mountor",itau_w,zx_tmp_2d, &
-                     iim*jjmp1,ndex2d)
-!
-      ENDIF !(1.EQ.0) THEN
-!
-!
-      CALL gr_fi_dyn(1,klon,iim+1,jjm+1,airephy,airedyn)
+
+
+!#ifdef _YM_UNUSED_TO_SUPRESS_
+!ym !
+!ym!IM 151004 BEG
+!ym      IF(1.EQ.0) THEN
+!ym!
+!ym      DO k=1, klev
+!ym      DO i=1, klon
+!ym       zx_tmp_fi3d(i,k)=u_seri(i,k)*RA*cos(pir* rlat(i))
+!ym      ENDDO
+!ym      ENDDO
+!ym!
+!ym      CALL moyglo_pondaima(klon, klev, zx_tmp_fi3d,  &
+!ym           airephy, paprs, moyglo)
+!ym      zx_tmp_fi2d(1:klon)=moyglo
+!ym!
+!ym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
+!ym      CALL histwrite(nid_day_seri,"momang",itau_w,zx_tmp_2d, &
+!ym                     iim*jjmp1,ndex2d)
+!ym!
+!ym! friction torque
+!ym!
+!ym      DO i=1, klon
+!ym       zx_tmp_fi2d(i)=zxfluxu(i,1)*RA* cos(pir* rlat(i))
+!ym      ENDDO
+!ym!
+!ym      ok_msk=.FALSE.
+!ym      CALL moyglo_pondaire(klon, zx_tmp_fi2d, airephy,  &
+!ym           ok_msk, msk, moyglo)
+!ym      zx_tmp_fi2d(1:klon)=moyglo
+!ym!
+!ym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
+!ym      CALL histwrite(nid_day_seri,"frictor",itau_w,zx_tmp_2d, &
+!ym                     iim*jjmp1,ndex2d)
+!ym!
+!ym! mountain torque
+!ym!
+!ym!IM 190504 BEG
+!ym      CALL gr_fi_dyn(1,klon,iim+1,jjm+1,airephy,airedyn)
+!ym      CALL gr_fi_dyn(klev+1,klon,iim+1,jjm+1,paprs,padyn)
+!ym      CALL gr_fi_dyn(1,klon,iim+1,jjm+1,rlat,rlatdyn)
+!ym      
+!ym
+!ym      mountor=0.
+!ym      airetot=0.
+!ym      DO j = 1, jjmp1
+!ym       DO i = 1, iim+1
+!ym        ij=i+(iim+1)*(j-1)
+!ym        zx_tmp(ij)=0.
+!ym        DO k = 1, klev
+!ym         zx_tmp(ij)=zx_tmp(ij)+dudyn(i,j,k)*airedyn(i,j)* &
+!ym                    (padyn(i,j,k+1)-padyn(i,j,k))/RG
+!ym         airetot=airetot+airedyn(i,j)
+!ym        ENDDO
+!ym!IM 190504 mountor=mountor+zx_tmp(ij)*airedyn(i,j)*RA*
+!ym        mountor=mountor+zx_tmp(ij)*RA* &
+!ym                 cos(pir* rlatdyn(i,j))
+!ym       ENDDO
+!ym      ENDDO
+!ym!IM 151004 BEG
+!ym      IF(itap.EQ.1) PRINT*,'airetot=',airetot,airetot/klev
+!ym!IM 151004 END
+!ym!IM 190504      mountor=mountor/(airetot*airetot)
+!ym      mountor=mountor/airetot
+!ym!
+!ym!IM 190504 END
+!ym      zx_tmp_2d(1:iim,1:jjmp1)=mountor
+!ym      CALL histwrite(nid_day_seri,"mountor",itau_w,zx_tmp_2d, &
+!ym                     iim*jjmp1,ndex2d)
+!ym!
+!ym      ENDIF !(1.EQ.0) THEN
+!ym!
+!ym!
+!ym
+!ym      CALL gr_fi_dyn(1,klon,iim+1,jjm+1,airephy,airedyn)
+!#endif      
+      
       CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
       airetot=0.
@@ -133 +141 @@
 !     ENDDO !j
 !
+
+
+      airetot=0.
       DO i=1, klon
        airetot=airetot+airephy(i)