Changeset 694 for LMDZ4/trunk/libf/phylmd
- Timestamp:
- Apr 4, 2006, 5:42:55 PM (19 years ago)
- Location:
- LMDZ4/trunk/libf/phylmd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/physiq.F
r687 r694 42 42 #define histmth 43 43 #define histins 44 c #define histISCCP 45 c #define histmthNMC 44 #define histmthNMC 45 #define histISCCP 46 46 c====================================================================== 47 47 c modif ( P. Le Van , 12/10/98 ) … … 1152 1152 REAL zx_tmp_fi2d(klon) ! variable temporaire grille physique 1153 1153 REAL zx_tmp_fi3d(klon,klev) ! variable temporaire pour champs 3D 1154 #ifdef histmthNMC 1155 REAL zx_tmp_NC(iim,jjmp1,nlevSTD) 1156 REAL zx_tmp_fiNC(klon,nlevSTD) 1157 #endif 1154 1158 REAL*8 zx_tmp2_fi3d(klon,klev) ! variable temporaire pour champs 3D 1155 1159 REAL zx_tmp_2d(iim,jjmp1), zx_tmp_3d(iim,jjmp1,klev) -
LMDZ4/trunk/libf/phylmd/write_histmthNMC.h
r684 r694 14 14 ccc 15 15 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,tsumSTD(:,:,2), 16 $ zx_tmp_ 3d)17 CALL histwrite(nid_nmc,"temp",itau_w,zx_tmp_ 3d,16 $ zx_tmp_NC) 17 CALL histwrite(nid_nmc,"temp",itau_w,zx_tmp_NC, 18 18 $ iim*jjmp1*nlevSTD,ndex3d) 19 19 c 20 20 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,phisumSTD(:,:,2), 21 $ zx_tmp_ 3d)22 CALL histwrite(nid_nmc,"phi",itau_w,zx_tmp_ 3d,21 $ zx_tmp_NC) 22 CALL histwrite(nid_nmc,"phi",itau_w,zx_tmp_NC, 23 23 $ iim*jjmp1*nlevSTD,ndex3d) 24 24 c 25 25 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,qsumSTD(:,:,2), 26 $ zx_tmp_ 3d)27 CALL histwrite(nid_nmc,"q",itau_w,zx_tmp_ 3d,26 $ zx_tmp_NC) 27 CALL histwrite(nid_nmc,"q",itau_w,zx_tmp_NC, 28 28 $ iim*jjmp1*nlevSTD,ndex3d) 29 29 c 30 30 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,rhsumSTD(:,:,2), 31 $ zx_tmp_ 3d)32 CALL histwrite(nid_nmc,"rh",itau_w,zx_tmp_ 3d,31 $ zx_tmp_NC) 32 CALL histwrite(nid_nmc,"rh",itau_w,zx_tmp_NC, 33 33 $ iim*jjmp1*nlevSTD,ndex3d) 34 34 c 35 35 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,usumSTD(:,:,2), 36 $ zx_tmp_ 3d)37 CALL histwrite(nid_nmc,"u",itau_w,zx_tmp_ 3d,36 $ zx_tmp_NC) 37 CALL histwrite(nid_nmc,"u",itau_w,zx_tmp_NC, 38 38 $ iim*jjmp1*nlevSTD,ndex3d) 39 39 c 40 40 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vsumSTD(:,:,2), 41 $ zx_tmp_ 3d)42 CALL histwrite(nid_nmc,"v",itau_w,zx_tmp_ 3d,41 $ zx_tmp_NC) 42 CALL histwrite(nid_nmc,"v",itau_w,zx_tmp_NC, 43 43 $ iim*jjmp1*nlevSTD,ndex3d) 44 44 ccc … … 46 46 ccc 47 47 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,tsumSTD(:,:,3), 48 $ zx_tmp_ 3d)49 CALL histwrite(nid_nmc,"temp",itau_w,zx_tmp_ 3d,48 $ zx_tmp_NC) 49 CALL histwrite(nid_nmc,"temp",itau_w,zx_tmp_NC, 50 50 $ iim*jjmp1*nlevSTD,ndex3d) 51 51 c 52 52 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,phisumSTD(:,:,3), 53 $ zx_tmp_ 3d)54 CALL histwrite(nid_nmc,"phi",itau_w,zx_tmp_ 3d,53 $ zx_tmp_NC) 54 CALL histwrite(nid_nmc,"phi",itau_w,zx_tmp_NC, 55 55 $ iim*jjmp1*nlevSTD,ndex3d) 56 56 c 57 57 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,qsumSTD(:,:,3), 58 $ zx_tmp_ 3d)59 CALL histwrite(nid_nmc,"q",itau_w,zx_tmp_ 3d,58 $ zx_tmp_NC) 59 CALL histwrite(nid_nmc,"q",itau_w,zx_tmp_NC, 60 60 $ iim*jjmp1*nlevSTD,ndex3d) 61 61 c 62 62 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,rhsumSTD(:,:,3), 63 $ zx_tmp_ 3d)64 CALL histwrite(nid_nmc,"rh",itau_w,zx_tmp_ 3d,63 $ zx_tmp_NC) 64 CALL histwrite(nid_nmc,"rh",itau_w,zx_tmp_NC, 65 65 $ iim*jjmp1*nlevSTD,ndex3d) 66 66 c 67 67 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,usumSTD(:,:,3), 68 $ zx_tmp_ 3d)69 CALL histwrite(nid_nmc,"u",itau_w,zx_tmp_ 3d,68 $ zx_tmp_NC) 69 CALL histwrite(nid_nmc,"u",itau_w,zx_tmp_NC, 70 70 $ iim*jjmp1*nlevSTD,ndex3d) 71 71 c 72 72 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vsumSTD(:,:,3), 73 $ zx_tmp_ 3d)74 CALL histwrite(nid_nmc,"v",itau_w,zx_tmp_ 3d,73 $ zx_tmp_NC) 74 CALL histwrite(nid_nmc,"v",itau_w,zx_tmp_NC, 75 75 $ iim*jjmp1*nlevSTD,ndex3d) 76 76 c 77 77 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,wsumSTD(:,:,3), 78 $ zx_tmp_ 3d)79 CALL histwrite(nid_nmc,"w",itau_w,zx_tmp_ 3d,78 $ zx_tmp_NC) 79 CALL histwrite(nid_nmc,"w",itau_w,zx_tmp_NC, 80 80 $ iim*jjmp1*nlevSTD,ndex3d) 81 81 c … … 83 83 DO i=1, klon 84 84 IF(tnondef(i,k,3).NE.1.E+20) THEN 85 zx_tmp_fi 3d(i,k) = (100.*tnondef(i,k,3))/ecrit_hf2mth85 zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,3))/ecrit_hf2mth 86 86 ELSE 87 zx_tmp_fi 3d(i,k) = 1.E+2087 zx_tmp_fiNC(i,k) = 1.E+20 88 88 ENDIF 89 89 ENDDO 90 90 ENDDO !k=1, nlevSTD 91 91 c 92 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,zx_tmp_fi 3d,zx_tmp_3d)93 CALL histwrite(nid_nmc,"psbg",itau_w,zx_tmp_ 3d,92 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,zx_tmp_fiNC,zx_tmp_NC) 93 CALL histwrite(nid_nmc,"psbg",itau_w,zx_tmp_NC, 94 94 $ iim*jjmp1*nlevSTD,ndex3d) 95 95 c 96 96 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,uvsumSTD(:,:,3), 97 $ zx_tmp_ 3d)98 CALL histwrite(nid_nmc,"uv",itau_w,zx_tmp_ 3d,97 $ zx_tmp_NC) 98 CALL histwrite(nid_nmc,"uv",itau_w,zx_tmp_NC, 99 99 $ iim*jjmp1*nlevSTD,ndex3d) 100 100 c 101 101 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vqsumSTD(:,:,3), 102 $ zx_tmp_ 3d)103 CALL histwrite(nid_nmc,"vq",itau_w,zx_tmp_ 3d,102 $ zx_tmp_NC) 103 CALL histwrite(nid_nmc,"vq",itau_w,zx_tmp_NC, 104 104 $ iim*jjmp1*nlevSTD,ndex3d) 105 105 c 106 106 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vTsumSTD(:,:,3), 107 $ zx_tmp_ 3d)108 CALL histwrite(nid_nmc,"vT",itau_w,zx_tmp_ 3d,107 $ zx_tmp_NC) 108 CALL histwrite(nid_nmc,"vT",itau_w,zx_tmp_NC, 109 109 $ iim*jjmp1*nlevSTD,ndex3d) 110 110 c 111 111 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1, wqsumSTD(:,:,3), 112 $ zx_tmp_ 3d)113 CALL histwrite(nid_nmc,"wq",itau_w,zx_tmp_ 3d,112 $ zx_tmp_NC) 113 CALL histwrite(nid_nmc,"wq",itau_w,zx_tmp_NC, 114 114 $ iim*jjmp1*nlevSTD,ndex3d) 115 115 c 116 116 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vphisumSTD(:,:,3), 117 $ zx_tmp_ 3d)118 CALL histwrite(nid_nmc,"vphi",itau_w,zx_tmp_ 3d,117 $ zx_tmp_NC) 118 CALL histwrite(nid_nmc,"vphi",itau_w,zx_tmp_NC, 119 119 $ iim*jjmp1*nlevSTD,ndex3d) 120 120 c 121 121 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,wTsumSTD(:,:,3), 122 $ zx_tmp_ 3d)123 CALL histwrite(nid_nmc,"wT",itau_w,zx_tmp_ 3d,122 $ zx_tmp_NC) 123 CALL histwrite(nid_nmc,"wT",itau_w,zx_tmp_NC, 124 124 $ iim*jjmp1*nlevSTD,ndex3d) 125 125 c 126 126 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,u2sumSTD(:,:,3), 127 $ zx_tmp_ 3d)128 CALL histwrite(nid_nmc,"uxu",itau_w,zx_tmp_ 3d,127 $ zx_tmp_NC) 128 CALL histwrite(nid_nmc,"uxu",itau_w,zx_tmp_NC, 129 129 $ iim*jjmp1*nlevSTD,ndex3d) 130 130 c 131 131 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,v2sumSTD(:,:,3), 132 $ zx_tmp_ 3d)133 CALL histwrite(nid_nmc,"vxv",itau_w,zx_tmp_ 3d,132 $ zx_tmp_NC) 133 CALL histwrite(nid_nmc,"vxv",itau_w,zx_tmp_NC, 134 134 $ iim*jjmp1*nlevSTD,ndex3d) 135 135 c 136 136 CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,T2sumSTD(:,:,3), 137 $ zx_tmp_ 3d)138 CALL histwrite(nid_nmc,"TxT",itau_w,zx_tmp_ 3d,137 $ zx_tmp_NC) 138 CALL histwrite(nid_nmc,"TxT",itau_w,zx_tmp_NC, 139 139 $ iim*jjmp1*nlevSTD,ndex3d) 140 140 c
Note: See TracChangeset
for help on using the changeset viewer.