Changeset 827 for LMDZ4/trunk/libf
- Timestamp:
- Aug 28, 2007, 2:17:59 PM (17 years ago)
- Location:
- LMDZ4/trunk/libf/phylmd
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/calcul_simulISCCP.h
r766 r827 10 10 c 11 11 nbapp_isccp=30 !appel toutes les 15h 12 isccppas=NINT((itap*dtime)/3600.) !Nb. d'heures de la physique12 cIM 170107 isccppas=NINT((itap*dtime)/3600.) !Nb. d'heures de la physique 13 13 freqin_pdt(n)=ifreq_isccp(n) 14 14 c -
LMDZ4/trunk/libf/phylmd/ini_histISCCP.h
r776 r827 15 15 cIM 300505 zstophy = dtime 16 16 c appel du simulateur toutes les 3heures 17 zcals(1) = dtime *6. !toutes les 3h (en s) 17 !IM on lit la frequence d'appel dans physiq.def 18 ! zcals(1) = dtime *6. !toutes les 3h (en s) 19 zcals(1) = freq_ISCCP !toutes les freq_ISCCP secondes 18 20 DO n=1, napisccp 19 21 zcalh(n) = zcals(n)/3600. !stoutes les Xh (en heures) … … 29 31 c IF(freqout_isccp.EQ.1.) THEN 30 32 c ecriture jounaliere 31 zout_isccp(1) = ecrit_day !(en s) 33 !IM on ecrit les resultats du simulateur ISCCP toutes les 34 ! ecrit_ISCCP secondes zout_isccp(1) = ecrit_day !(en s) 35 zout_isccp(1) = ecrit_ISCCP !(en s) 32 36 c ecriture mensuelle 33 37 c zout = dtime * ecrit_mth !(en s) … … 37 41 c le nombre de sous-colonnes ncol : ncol=(100.*zcalh)/zoutd 38 42 ncol(n)=NINT((100.*zcalh(n))/zoutj(n)) 39 PRINT*,'n ncol',n,ncol(n) 43 IF(ncol(n).GT.ncolmx) THEN 44 PRINT*,'Warning: Augmenter le nombre colonnes du simulateur' 45 PRINT*,' ISCCP ncol=', ncol,' ncolmx=',ncolmx 46 c PRINT*,'n ncol',n,ncol(n) 47 CALL abort 48 ENDIF 40 49 c 41 50 DO l=1, ncol(n) … … 80 89 cym enddo 81 90 c 82 print*,'On stoke le fichier histISCCP instantannesur ',91 print*,'On stoke le fichier histISCCP sur ', 83 92 s imin_ins,imax_ins,jmin_ins,jmax_ins 84 93 cym print*,'On stoke le fichier histISCCP instantanne sur ', … … 94 103 CALL histbeg_phy("histISCCP.nc", itau_phy, zjulian, dtime, 95 104 . nhori, nid_isccp) 96 97 105 ENDIF !(1.EQ.0) THEN 98 106 c … … 106 114 CALL histbeg_phy("histISCCP.nc", itau_phy, zjulian, dtime, 107 115 . nhori, nid_isccp) 108 109 116 c 110 117 IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN … … 128 135 c 129 136 DO k=1, kmaxm1 130 CALL histdef(nid_isccp, "cldISCCP_"//taulev(k)// typinout(n),137 CALL histdef(nid_isccp, "cldISCCP_"//taulev(k)//verticaxe(n), 131 138 . "LMDZ ISCCP cld", "%", 132 139 . iim, jj_nb,nhori,lmaxm1,1,lmaxm1,nvert,32, … … 134 141 ENDDO 135 142 c 136 CALL histdef(nid_isccp, "nsunlit"// typinout(n),143 CALL histdef(nid_isccp, "nsunlit"//verticaxe(n), 137 144 . "Nb of calls with sunlit ", "%", 138 145 . iim, jj_nb,nhori,1,1,1,-99,32, … … 150 157 DO l=1, lmaxm1 151 158 c 152 CALL histdef(nid_isccp, pclev(l)//taulev(k)// typinout(n),159 CALL histdef(nid_isccp, pclev(l)//taulev(k)//verticaxe(n), 153 160 . "LMDZ ISCCP cld "//cnameisccp(l,k), "%", 154 161 . iim, jj_nb,nhori,1,1,1,-99,32, … … 159 166 c 160 167 c print*,'n=',n,' avant histdef(..Nb of calls sunlit' 161 CALL histdef(nid_isccp, "nsunlit"// typinout(n),168 CALL histdef(nid_isccp, "nsunlit"//verticaxe(n), 162 169 . "Nb of calls with sunlit ", "%", 163 170 . iim, jj_nb,nhori,1,1,1,-99,32, -
LMDZ4/trunk/libf/phylmd/write_histISCCP.h
r776 r827 20 20 cIM: champ 3d : (lon,lat,pres) pour un tau fixe 21 21 c 22 CALL histwrite_phy(nid_isccp,"cldISCCP_"//taulev(k)// typinout(n),22 CALL histwrite_phy(nid_isccp,"cldISCCP_"//taulev(k)//verticaxe(n), 23 23 . itau_w,zx_tmp_fi3d) 24 24 ENDDO !k 25 25 c 26 26 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,nbsunlit(1,:,n),zx_tmp_2d) 27 CALL histwrite_phy(nid_isccp,"nsunlit"// typinout(n),itau_w,27 CALL histwrite_phy(nid_isccp,"nsunlit"//verticaxe(n),itau_w, 28 28 . nbsunlit(1,:,n)) 29 29 c … … 55 55 cIM: champ 2d : (lon,lat) pour un tau et une pc fixes 56 56 c 57 CALL histwrite_phy(nid_isccp,pclev(l)//taulev(k)// typinout(n),58 . itau_w,fq_is_true(:, :,l,n))57 CALL histwrite_phy(nid_isccp,pclev(l)//taulev(k)//verticaxe(n), 58 . itau_w,fq_is_true(:,k,l,n)) 59 59 ENDDO !l 60 60 ENDDO !k … … 62 62 c print*,'n=',n,' write_ISCCP avant nbsunlit' 63 63 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,nbsunlit(1,:,n),zx_tmp_2d) 64 CALL histwrite_phy(nid_isccp,"nsunlit"// typinout(n),64 CALL histwrite_phy(nid_isccp,"nsunlit"//verticaxe(n), 65 65 . itau_w,nbsunlit(1,:,n)) 66 66 c
Note: See TracChangeset
for help on using the changeset viewer.