Index: LMDZ6/trunk/libf/obsolete/ini_bilKP_ave.h =================================================================== --- LMDZ6/trunk/libf/obsolete/ini_bilKP_ave.h (revision 5311) +++ LMDZ6/trunk/libf/obsolete/ini_bilKP_ave.h (revision 5311) @@ -0,0 +1,259 @@ +c +c $Id$ +c + IF (ok_journe) THEN +c + zsto = dtime + zout = ecrit_day + typeval=tave +c + idayref = day_ref + CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) +c +cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) +cym DO i = 1, iim +cym zx_lon(i,1) = rlon(i+1) +cym zx_lon(i,jjmp1) = rlon(i+1) +cym ENDDO + DO ll=1,klev + znivsig(ll)=REAL(ll) + ENDDO +cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) +cym write(*,*)'zx_lon = ',zx_lon(:,1) +cym write(*,*)'zx_lat = ',zx_lat(1,:) +cym CALL histbeg("histbilKP_ave", iim,zx_lon(:,1), jjmp1, +cym . zx_lat(1,:), +cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, +cym . nhori, nid_bilKPave) + CALL histbeg_phy("histbilKP_ave", itau_phy, zjulian, dtime, + . nhori, nid_bilKPave) + + write(*,*)'Journee ', itau_phy, zjulian + CALL histvert(nid_bilKPave, "presnivs", + . "Vertical levels","mb", + . klev, presnivs/100., nvert) +c +c +c Champs 3D: +c + CALL histdef(nid_bilKPave,"ue", + . "Zonal energy transport","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"ve", + . "Merid energy transport","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"uq", + . "Zonal humidity transport","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"vq", + . "Merid humidity transport","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c +c Champs 3D: +c + CALL histdef(nid_bilKPave,"temp", + . "Air temperature","K", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"ovap", + . "Specific humidity","Kg/Kg", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"geop", + . "Geopotential height","m", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"vitu", + . "Zonal wind","m/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"vitv", + . "Meridional wind","m/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "vitw", + . "Vertical wind", "m/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "pres", + . "Inter-Layer Air pressure", + . "Pa", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "play", + . "Mean-Layer Air pressure", + . "Pa", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "oliq", + . "Liquid water content", + . "kg/kg", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "dtdyn", + . "Dynamics dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "dqdyn", + . "Dynamics dQ", "Kg/Kg/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "dtcon", + . "Convection dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "ducon", + . "Convection du", "m/s2", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "dvcon", + . "Convection dv", "m/s2", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"dqcon", + . "Convection dQ","Kg/Kg/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "dtlsc", + . "Condensation dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"dqlsc", + . "Condensation dQ","Kg/Kg/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"dtvdf", + . "Boundary-layer dT","K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "dqvdf", + . "Boundary-layer dQ", + . "Kg/Kg/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"dtajs", + . "Ajustement sec dT","K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "dqajs", + . "Ajustement sec dQ", + . "Kg/Kg/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "dteva", + . "Reevaporation dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"dqeva", + . "Reevaporation dQ", + . "Kg/Kg/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) + +c + CALL histdef(nid_bilKPave, "dtswr", + . "SW radiation dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "dtsw0", + . "SW radiation dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "dtlwr", + . "LW radiation dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "dtlw0", + . "LW radiation dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"duvdf", + . "Boundary-layer dU","m/s2", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave,"dvvdf", + . "Boundary-layer dV","m/s2", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + IF (ok_orodr) THEN + IF (ok_orolf) THEN + CALL histdef(nid_bilKPave, "duoli", + . "Orography dU", "m/s2", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPave, "dvoli", + . "Orography dV", "m/s2", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + ENDIF + ENDIF +C + CALL histdef(nid_bilKPave, "duphy", + . "Physiq dU","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +C + CALL histdef(nid_bilKPave, "dvphy", + . "Physiq dV","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +C + CALL histdef(nid_bilKPave, "dtphy", + . "Physiq dT","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +C + CALL histdef(nid_bilKPave, "dqphy", + . "Physiq dQ","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +C + CALL histdef(nid_bilKPave, "dqlphy", + . "Physiq dQl","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +C +C + CALL histend(nid_bilKPave) +c + ndex2d = 0 + ndex3d = 0 +c + ENDIF ! fin de test sur ok_journe Index: LMDZ6/trunk/libf/obsolete/ini_bilKP_ins.h =================================================================== --- LMDZ6/trunk/libf/obsolete/ini_bilKP_ins.h (revision 5311) +++ LMDZ6/trunk/libf/obsolete/ini_bilKP_ins.h (revision 5311) @@ -0,0 +1,326 @@ +c +c $Id$ +c + IF (ok_journe) THEN +c + zsto = dtime + zout = dtime + typeval=tinst +c + idayref = day_ref + CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) +c +cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) +cym DO i = 1, iim +cym zx_lon(i,1) = rlon(i+1) +cym zx_lon(i,jjmp1) = rlon(i+1) +cym ENDDO + DO ll=1,klev + znivsig(ll)=REAL(ll) + ENDDO +cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) +cym write(*,*)'zx_lon = ',zx_lon(:,1) +cym write(*,*)'zx_lat = ',zx_lat(1,:) +c +cIM 280405 BEG +c +cIM cf. AM 081204 BEG region + imin_ins=1 + imax_ins=iim + jmin_ins=1 + jmax_ins=jjmp1 +cym do i=1,iim-1 +cym if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i +cym if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1 +cym enddo +cym do j=1,jjmp1 +cym if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j +cym if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j +cym enddo +c + print*,'On stoke le fichier bilKP instantanne sur ', + s imin_ins,imax_ins,jmin_ins,jmax_ins + print*,'On stoke le fichier bilKP instantanne sur ', + s zx_lon(imin_ins,1),zx_lon(imax_ins,1), + s zx_lat(1,jmin_ins),zx_lat(1,jmax_ins) +cIM cf. AM 081204 END region +c +cIM 280405 END +c +cym IF(1.EQ.0) THEN +cym CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), jjmp1, +cym . zx_lat(1,:), +cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, +cym . nhori, nid_bilKPins) + ENDIF +c +cIM 280405 BEG +c +cIM cf. AM 081204 BEG region +cym CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), +cym . jjmp1,zx_lat(1,:), +cym . imin_ins,imax_ins-imin_ins+1, +cym . jmin_ins,jmax_ins-jmin_ins+1, +cym . itau_phy, zjulian, dtime, +cym . nhori, nid_bilKPins) + CALL histbeg_phy("histbilKP_ins", itau_phy, zjulian, dtime, + . nhori, nid_bilKPins) +cIM 081204 END +c +cIM 280405 END +c + write(*,*)'Journee ', itau_phy, zjulian + CALL histvert(nid_bilKPins, "presnivs", + . "Vertical levels","mb", + . klev, presnivs/100., nvert) +c +c Champs 3D: +c + CALL histdef(nid_bilKPins,"ue", + . "Zonal energy transport","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"ve", + . "Merid energy transport","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"uq", + . "Zonal humidity transport","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"vq", + . "Merid humidity transport","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c +c Champs 3D: +c + CALL histdef(nid_bilKPins, "temp", + . "Air temperature", "K", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"ovap", + . "Specific humidity","Kg/Kg", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"geop", + . "Geopotential height", "m", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"vitu", + . "Zonal wind", "m/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"vitv", + . "Meridional wind", "m/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "vitw", + . "Vertical wind", "m/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "pres", + . "Inter-Layer Air pressure", + . "Pa", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "play", + . "Mean-Layer Air pressure", + . "Pa", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "oliq", + . "Liquid water content", + . "kg/kg", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "dtdyn", + . "Dynamics dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "dqdyn", + . "Dynamics dQ", "Kg/Kg/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "dtcon", + . "Convection dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "ducon", + . "Convection du", "m/s2", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "dvcon", + . "Convection dv", "m/s2", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"dqcon", + . "Convection dQ","Kg/Kg/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "dtlsc", + . "Condensation dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"dqlsc", + . "Condensation dQ","Kg/Kg/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"dtvdf", + . "Boundary-layer dT","K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "dqvdf", + . "Boundary-layer dQ", + . "Kg/Kg/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"dtajs", + . "Ajustement sec dT","K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"dqajs", + . "Ajustement sec dQ", + . "Kg/Kg/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"dteva", + . "Reevaporation dT","K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"dqeva", + . "Reevaporation dQ", + . "Kg/Kg/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) + +c + CALL histdef(nid_bilKPins, "dtswr", + . "SW radiation dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "dtsw0", + . "SW radiation dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "dtlwr", + . "LW radiation dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "dtlw0", + . "LW radiation dT", "K/s", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"duvdf", + . "Boundary-layer dU","m/s2", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins,"dvvdf", + . "Boundary-layer dV","m/s2", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + IF (ok_orodr) THEN + IF (ok_orolf) THEN + CALL histdef(nid_bilKPins, "duoli", + . "Orography dU", "m/s2", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "dvoli", + . "Orography dV", "m/s2", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +c + ENDIF + ENDIF +C + CALL histdef(nid_bilKPins, "duphy", + . "Physiq dU","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +C + CALL histdef(nid_bilKPins, "dvphy", + . "Physiq dV","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +C + CALL histdef(nid_bilKPins, "dtphy", + . "Physiq dT","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +C + CALL histdef(nid_bilKPins, "dqphy", + . "Physiq dQ","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +C + CALL histdef(nid_bilKPins, "dqlphy", + . "Physiq dQl","-", + . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, + . typeval, zsto,zout) +cIM 280405 BEG +c +c Champs 2D: +c +c u850, v850 +c DO k=1, nlevSTD + DO k=1, 12 +c + IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k) + IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k) +c + IF(bb2.EQ."850") THEN +c + CALL histdef(nid_bilKPins, "u"//bb2, + . "Zonal wind "//bb2//"mb","m/s", + . iim,jjphy_nb,nhori, 1,1,1, -99, 32, + . typeval, zsto,zout) +c + CALL histdef(nid_bilKPins, "v"//bb2, + . "Meridional wind "//bb2//"mb","m/s", + . iim,jjphy_nb,nhori, 1,1,1, -99, 32, + . typeval, zsto,zout) +c + ENDIF !(bb2.EQ."850") +c + ENDDO !k=1, 12 +c +cIM 280405 END +c + CALL histend(nid_bilKPins) +c + ndex2d = 0 + ndex3d = 0 +c + ENDIF ! fin de test sur ok_journe Index: LMDZ6/trunk/libf/obsolete/ini_coord_REGDYN.h =================================================================== --- LMDZ6/trunk/libf/obsolete/ini_coord_REGDYN.h (revision 5311) +++ LMDZ6/trunk/libf/obsolete/ini_coord_REGDYN.h (revision 5311) @@ -0,0 +1,81 @@ +c +c $Header$ +c + nsrf=3 + DO nreg=1, nbregdyn + DO i=1, klon + +c IF (debut) THEN + IF(rlon(i).LT.0.) THEN + rlonPOS(i)=rlon(i)+360. + ELSE + rlonPOS(i)=rlon(i) + ENDIF +c ENDIF + + pct_ocean(i,nreg)=0 + +c test si c'est 1 point d'ocean + IF(pctsrf(i,nsrf).EQ.1.) THEN + + IF(nreg.EQ.1) THEN + +c TROP + IF(rlat(i).GE.-30.AND.rlat(i).LE.30.) THEN + pct_ocean(i,nreg)=1 + ENDIF + +c PACIFIQUE NORD + ELSEIF(nreg.EQ.2) THEN + IF(rlat(i).GE.40.AND.rlat(i).LE.60.) THEN + IF(rlonPOS(i).GE.160..AND.rlonPOS(i).LE.235.) THEN + pct_ocean(i,nreg)=1 + ENDIF + ENDIF +c CALIFORNIE ST-CU + ELSEIF(nreg.EQ.3) THEN + IF(rlonPOS(i).GE.220..AND.rlonPOS(i).LE.250.) THEN + IF(rlat(i).GE.15.AND.rlat(i).LE.35.) THEN + pct_ocean(i,nreg)=1 + ENDIF + ENDIF +c HAWAI + ELSEIF(nreg.EQ.4) THEN + IF(rlonPOS(i).GE.180..AND.rlonPOS(i).LE.220.) THEN + IF(rlat(i).GE.15.AND.rlat(i).LE.35.) THEN + pct_ocean(i,nreg)=1 + ENDIF + ENDIF +c WARM POOL + ELSEIF(nreg.EQ.5) THEN + IF(rlonPOS(i).GE.70..AND.rlonPOS(i).LE.150.) THEN + IF(rlat(i).GE.-5.AND.rlat(i).LE.20.) THEN + pct_ocean(i,nreg)=1 + ENDIF + ENDIF + ENDIF !nbregdyn +c TROP +c IF(rlat(i).GE.-30.AND.rlat(i).LE.30.) THEN +c pct_ocean(i)=.TRUE. +c WRITE(*,*) 'pct_ocean =',i, rlon(i), rlat(i) +c ENDIF !lon +c ENDIF !lat + + ENDIF !pctsrf + ENDDO !klon + ENDDO !nbregdyn +cIM 190504 ENDIF !ok_regdyn + +cIM somme de toutes les nhistoW BEG + IF (debut) THEN + DO nreg = 1, nbregdyn + DO k = 1, kmaxm1 + DO l = 1, lmaxm1 + DO iw = 1, iwmax + nhistoWt(k,l,iw,nreg)=0. + ENDDO !iw + ENDDO !l + ENDDO !k + ENDDO !nreg + ENDIF !(debut) THEN +cIM 190504 BEG Index: LMDZ6/trunk/libf/obsolete/ini_histREGDYN.h =================================================================== --- LMDZ6/trunk/libf/obsolete/ini_histREGDYN.h (revision 5311) +++ LMDZ6/trunk/libf/obsolete/ini_histREGDYN.h (revision 5311) @@ -0,0 +1,127 @@ +! +! $Header$ +! + + IF (ok_regdyn) THEN + + if (is_sequential) then +c +cIM PRINT*, 'La frequence de sortie REGDYN est de ', ecrit_mth +c + idayref = day_ref + CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) +c +c axe vertical pour les differents niveaux des histogrammes + DO iw=1, iwmax + zx_o500(iw)=wmin+(iw-1./2.)*pas_w + ENDDO + + CALL histbeg("histREGDYN", kmaxm1,zx_tau, lmaxm1,zx_pc, + . 1,kmaxm1,1,lmaxm1, itau_phy, zjulian, dtime, + . nhoriRD, nid_regdyn) + + CALL histvert(nid_regdyn, "omeganivs", "Omega levels", + . "mb/day", + . iwmax, zx_o500, komega) +c +c pour les champs instantannes, il faut mettre la meme valeur pour +c zout et zsto. +c dtime est passe par ailleurs a histbeg +c +c zout = dtime * REAL(NINT(86400./dtime*ecrit_regdyn)) +c zsto = zout +c print*,'zout,zsto=',zout,zsto +c +c stockage a chaque pas de temps de la physique +c + zstophy = dtime +cIM 020904 zstophy = dtime * nbapp_isccp + +c ecriture mensuelle +c + zout = dtime * ecrit_mth +cIM 020904 +c zout = dtime * ecrit_day +c zout = dtime * REAL(NINT(86400./dtime*ecrit_regdyn)) + +c +c Champs 3D: +c +c TROP + CALL histdef(nid_regdyn, "hw1", "Tropics Histogram ", "%", + & kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, + & "ave(X)", zstophy,zout) + + CALL histdef(nid_regdyn, "nh1", "Nb of pixels Tropics Histo", + & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, + & 32,"ave(X)", zstophy,zout) +c + + CALL histdef(nid_regdyn, "nht1", + & "Total Nb pixels Tropics Histo", + & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, + & 32,"ave(X)", zstophy,zout) +c +c PAN + CALL histdef(nid_regdyn, "hw2", "North Pacific Histogram", "%", + & kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, + & "ave(X)", zstophy,zout) + + CALL histdef(nid_regdyn, "nh2", "Nb of pixels North Pacific", + & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, + & 32,"ave(X)", zstophy,zout) +c + + CALL histdef(nid_regdyn, "nht2", + & "Total Nb pixels North Pacific Histo", + & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, + & 32,"ave(X)", zstophy,zout) +c CAL + CALL histdef(nid_regdyn, "hw3", "California Histogram", "%", + & kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, + & "ave(X)", zstophy,zout) + + CALL histdef(nid_regdyn, "nh3", + & "Nb of pixels California Histo", + & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, + & 32,"ave(X)", zstophy,zout) +c + + CALL histdef(nid_regdyn, "nht3", + & "Total Nb pixels California Histo", + & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, + & 32,"ave(X)", zstophy,zout) +c HAW + CALL histdef(nid_regdyn, "hw4", "Hawai Histogram", "%", + & kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, + & "ave(X)", zstophy,zout) + + CALL histdef(nid_regdyn, "nh4", "Nb of pixels Hawai Histo", + & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, + & 32,"ave(X)", zstophy,zout) +c + + CALL histdef(nid_regdyn, "nht4", + & "Total Nb pixels Hawai Histo", + & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, + & 32,"ave(X)", zstophy,zout) +c WAP + CALL histdef(nid_regdyn, "hw5", "Warm Pool Histogram", "%", + & kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, + & "ave(X)", zstophy,zout) + + CALL histdef(nid_regdyn, "nh5", "Nb of pixels Warm Pool Histo", + & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, + & 32,"ave(X)", zstophy,zout) +c + + CALL histdef(nid_regdyn, "nht5", + & "Total Nb pixels Warm Pool Histo", + & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, + & 32,"ave(X)", zstophy,zout) +c + CALL histend(nid_regdyn) + + endif ! is_sequential + + endif ! ok_regdyn Index: LMDZ6/trunk/libf/obsolete/ini_histday_seri.h =================================================================== --- LMDZ6/trunk/libf/obsolete/ini_histday_seri.h (revision 5311) +++ LMDZ6/trunk/libf/obsolete/ini_histday_seri.h (revision 5311) @@ -0,0 +1,131 @@ +! +! $Id$ +! +!ym Ne fonctionnera pas en mode parallele + IF (is_sequential) THEN + + IF (type_run.EQ."AMIP") THEN +! + zstophy = pdtphys + zout = ecrit_day +! + idayref = day_ref + CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) +! + CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,longitude_deg,zx_lon) + DO i = 1, nbp_lon + zx_lon(i,1) = longitude_deg(i+1) + zx_lon(i,nbp_lat) = longitude_deg(i+1) + ENDDO + DO ll=1,klev + znivsig(ll)=REAL(ll) + ENDDO + CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,latitude_deg,zx_lat) +! + imin_debut=1 + nbpti=1 + jmin_debut=1 + nbptj=1 +! + CALL histbeg("histday_seri.nc", & + nbp_lon,zx_lon(:,1),nbp_lat,zx_lat(1,:), & + imin_debut,nbpti,jmin_debut,nbptj, & + itau_phy, zjulian, dtime, & + nhori, nid_day_seri) +! + CALL histvert(nid_day_seri, "presnivs", & + "Vertical levels","mb", & + klev, presnivs/100., nvert) +! + CALL histdef(nid_day_seri, "bilTOA", & + "Net radiation at model top", "W/m2", & + nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & + "ave(X)", zstophy,zout) +! + CALL histdef(nid_day_seri, "bils", & + "Net downward energy flux at surface","W/m2", & + nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & + "ave(X)", zstophy,zout) +! + CALL histdef(nid_day_seri, "ecin", & + "Total kinetic energy (per unit area)","J/m2", & + nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & + "ave(X)", zstophy,zout) +! +!IM 151004 BEG + IF(1.EQ.0) THEN +! + CALL histdef(nid_day_seri, "momang", & + "Total relative angular momentum (per unit area)", & + "kg/s", & + nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & + "ave(X)", zstophy,zout) +! + CALL histdef(nid_day_seri, "frictor", & + "Friction torque (per unit area)", "N/m", & + nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & + "ave(X)", zstophy,zout) +! + CALL histdef(nid_day_seri, "mountor", & + "Mountain torque (per unit area)", "N/m", & + nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & + "ave(X)", zstophy,zout) +! + ENDIF !(1.EQ.0) THEN +! + CALL histdef(nid_day_seri, "momang", & + "Axial angular momentum (per unit area)", & + "kg/s", & + nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & + "ave(X)", zstophy,zout) +! + CALL histdef(nid_day_seri, "torsfc", & + "Total surface torque (including mountain torque)", "N/m", & + nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & + "ave(X)", zstophy,zout) +! +!IM 151004 END +! + CALL histdef(nid_day_seri, "tamv", & + "Temperature (mass-weighted vert. ave)", "K", & + nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & + "ave(X)", zstophy,zout) +! + CALL histdef(nid_day_seri, "psol", & + "Surface pressure", "Pa", & + nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & + "ave(X)", zstophy,zout) +! + CALL histdef(nid_day_seri, "evap", & + "Evaporation and sublimation (per unit area)", & + "kg/(m2*s)", & + nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & + "ave(X)", zstophy,zout) +! +! call histdef(nid_day_seri, +! . "SnowFrac", +! . "Snow-covered area ", "%", +! . nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, +! . "ave(X)", zstophy,zout) +! +! CALL histdef(nid_day_seri, "snow_depth", +!IM 080904 . "Snow Depth (water equivalent)", "m", +!IM 191104 . "Snow Depth (water equivalent)", "kg/m2", +! . "Snow Mass", "kg/m2", +! . nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, +! . "ave(X)", zstophy,zout) +! + call histdef(nid_day_seri, & + "tsol_"//clnsurf(is_oce), & + "SST over open (ice-free) ocean ", "K", & + nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & + "ave(X)", zstophy,zout) +! +!================================================================= +! + CALL histend(nid_day_seri) +! +!================================================================= + ENDIF ! fin de test sur type_run.EQ.AMIP + + ENDIF ! is_sequential Index: LMDZ6/trunk/libf/phylmd/ini_bilKP_ave.h =================================================================== --- LMDZ6/trunk/libf/phylmd/ini_bilKP_ave.h (revision 5310) +++ (revision ) @@ -1,259 +1,0 @@ -c -c $Id$ -c - IF (ok_journe) THEN -c - zsto = dtime - zout = ecrit_day - typeval=tave -c - idayref = day_ref - CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) -c -cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) -cym DO i = 1, iim -cym zx_lon(i,1) = rlon(i+1) -cym zx_lon(i,jjmp1) = rlon(i+1) -cym ENDDO - DO ll=1,klev - znivsig(ll)=REAL(ll) - ENDDO -cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) -cym write(*,*)'zx_lon = ',zx_lon(:,1) -cym write(*,*)'zx_lat = ',zx_lat(1,:) -cym CALL histbeg("histbilKP_ave", iim,zx_lon(:,1), jjmp1, -cym . zx_lat(1,:), -cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, -cym . nhori, nid_bilKPave) - CALL histbeg_phy("histbilKP_ave", itau_phy, zjulian, dtime, - . nhori, nid_bilKPave) - - write(*,*)'Journee ', itau_phy, zjulian - CALL histvert(nid_bilKPave, "presnivs", - . "Vertical levels","mb", - . klev, presnivs/100., nvert) -c -c -c Champs 3D: -c - CALL histdef(nid_bilKPave,"ue", - . "Zonal energy transport","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"ve", - . "Merid energy transport","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"uq", - . "Zonal humidity transport","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"vq", - . "Merid humidity transport","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c -c Champs 3D: -c - CALL histdef(nid_bilKPave,"temp", - . "Air temperature","K", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"ovap", - . "Specific humidity","Kg/Kg", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"geop", - . "Geopotential height","m", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"vitu", - . "Zonal wind","m/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"vitv", - . "Meridional wind","m/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "vitw", - . "Vertical wind", "m/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "pres", - . "Inter-Layer Air pressure", - . "Pa", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "play", - . "Mean-Layer Air pressure", - . "Pa", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "oliq", - . "Liquid water content", - . "kg/kg", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "dtdyn", - . "Dynamics dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "dqdyn", - . "Dynamics dQ", "Kg/Kg/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "dtcon", - . "Convection dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "ducon", - . "Convection du", "m/s2", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "dvcon", - . "Convection dv", "m/s2", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"dqcon", - . "Convection dQ","Kg/Kg/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "dtlsc", - . "Condensation dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"dqlsc", - . "Condensation dQ","Kg/Kg/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"dtvdf", - . "Boundary-layer dT","K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "dqvdf", - . "Boundary-layer dQ", - . "Kg/Kg/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"dtajs", - . "Ajustement sec dT","K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "dqajs", - . "Ajustement sec dQ", - . "Kg/Kg/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "dteva", - . "Reevaporation dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"dqeva", - . "Reevaporation dQ", - . "Kg/Kg/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) - -c - CALL histdef(nid_bilKPave, "dtswr", - . "SW radiation dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "dtsw0", - . "SW radiation dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "dtlwr", - . "LW radiation dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "dtlw0", - . "LW radiation dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"duvdf", - . "Boundary-layer dU","m/s2", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave,"dvvdf", - . "Boundary-layer dV","m/s2", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - IF (ok_orodr) THEN - IF (ok_orolf) THEN - CALL histdef(nid_bilKPave, "duoli", - . "Orography dU", "m/s2", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPave, "dvoli", - . "Orography dV", "m/s2", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - ENDIF - ENDIF -C - CALL histdef(nid_bilKPave, "duphy", - . "Physiq dU","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -C - CALL histdef(nid_bilKPave, "dvphy", - . "Physiq dV","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -C - CALL histdef(nid_bilKPave, "dtphy", - . "Physiq dT","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -C - CALL histdef(nid_bilKPave, "dqphy", - . "Physiq dQ","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -C - CALL histdef(nid_bilKPave, "dqlphy", - . "Physiq dQl","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -C -C - CALL histend(nid_bilKPave) -c - ndex2d = 0 - ndex3d = 0 -c - ENDIF ! fin de test sur ok_journe Index: LMDZ6/trunk/libf/phylmd/ini_bilKP_ins.h =================================================================== --- LMDZ6/trunk/libf/phylmd/ini_bilKP_ins.h (revision 5310) +++ (revision ) @@ -1,326 +1,0 @@ -c -c $Id$ -c - IF (ok_journe) THEN -c - zsto = dtime - zout = dtime - typeval=tinst -c - idayref = day_ref - CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) -c -cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) -cym DO i = 1, iim -cym zx_lon(i,1) = rlon(i+1) -cym zx_lon(i,jjmp1) = rlon(i+1) -cym ENDDO - DO ll=1,klev - znivsig(ll)=REAL(ll) - ENDDO -cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) -cym write(*,*)'zx_lon = ',zx_lon(:,1) -cym write(*,*)'zx_lat = ',zx_lat(1,:) -c -cIM 280405 BEG -c -cIM cf. AM 081204 BEG region - imin_ins=1 - imax_ins=iim - jmin_ins=1 - jmax_ins=jjmp1 -cym do i=1,iim-1 -cym if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i -cym if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1 -cym enddo -cym do j=1,jjmp1 -cym if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j -cym if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j -cym enddo -c - print*,'On stoke le fichier bilKP instantanne sur ', - s imin_ins,imax_ins,jmin_ins,jmax_ins - print*,'On stoke le fichier bilKP instantanne sur ', - s zx_lon(imin_ins,1),zx_lon(imax_ins,1), - s zx_lat(1,jmin_ins),zx_lat(1,jmax_ins) -cIM cf. AM 081204 END region -c -cIM 280405 END -c -cym IF(1.EQ.0) THEN -cym CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), jjmp1, -cym . zx_lat(1,:), -cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, -cym . nhori, nid_bilKPins) - ENDIF -c -cIM 280405 BEG -c -cIM cf. AM 081204 BEG region -cym CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), -cym . jjmp1,zx_lat(1,:), -cym . imin_ins,imax_ins-imin_ins+1, -cym . jmin_ins,jmax_ins-jmin_ins+1, -cym . itau_phy, zjulian, dtime, -cym . nhori, nid_bilKPins) - CALL histbeg_phy("histbilKP_ins", itau_phy, zjulian, dtime, - . nhori, nid_bilKPins) -cIM 081204 END -c -cIM 280405 END -c - write(*,*)'Journee ', itau_phy, zjulian - CALL histvert(nid_bilKPins, "presnivs", - . "Vertical levels","mb", - . klev, presnivs/100., nvert) -c -c Champs 3D: -c - CALL histdef(nid_bilKPins,"ue", - . "Zonal energy transport","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"ve", - . "Merid energy transport","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"uq", - . "Zonal humidity transport","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"vq", - . "Merid humidity transport","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c -c Champs 3D: -c - CALL histdef(nid_bilKPins, "temp", - . "Air temperature", "K", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"ovap", - . "Specific humidity","Kg/Kg", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"geop", - . "Geopotential height", "m", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"vitu", - . "Zonal wind", "m/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"vitv", - . "Meridional wind", "m/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "vitw", - . "Vertical wind", "m/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "pres", - . "Inter-Layer Air pressure", - . "Pa", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "play", - . "Mean-Layer Air pressure", - . "Pa", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "oliq", - . "Liquid water content", - . "kg/kg", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "dtdyn", - . "Dynamics dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "dqdyn", - . "Dynamics dQ", "Kg/Kg/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "dtcon", - . "Convection dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "ducon", - . "Convection du", "m/s2", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "dvcon", - . "Convection dv", "m/s2", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"dqcon", - . "Convection dQ","Kg/Kg/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "dtlsc", - . "Condensation dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"dqlsc", - . "Condensation dQ","Kg/Kg/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"dtvdf", - . "Boundary-layer dT","K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "dqvdf", - . "Boundary-layer dQ", - . "Kg/Kg/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"dtajs", - . "Ajustement sec dT","K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"dqajs", - . "Ajustement sec dQ", - . "Kg/Kg/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"dteva", - . "Reevaporation dT","K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"dqeva", - . "Reevaporation dQ", - . "Kg/Kg/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) - -c - CALL histdef(nid_bilKPins, "dtswr", - . "SW radiation dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "dtsw0", - . "SW radiation dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "dtlwr", - . "LW radiation dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "dtlw0", - . "LW radiation dT", "K/s", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"duvdf", - . "Boundary-layer dU","m/s2", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins,"dvvdf", - . "Boundary-layer dV","m/s2", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - IF (ok_orodr) THEN - IF (ok_orolf) THEN - CALL histdef(nid_bilKPins, "duoli", - . "Orography dU", "m/s2", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "dvoli", - . "Orography dV", "m/s2", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -c - ENDIF - ENDIF -C - CALL histdef(nid_bilKPins, "duphy", - . "Physiq dU","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -C - CALL histdef(nid_bilKPins, "dvphy", - . "Physiq dV","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -C - CALL histdef(nid_bilKPins, "dtphy", - . "Physiq dT","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -C - CALL histdef(nid_bilKPins, "dqphy", - . "Physiq dQ","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -C - CALL histdef(nid_bilKPins, "dqlphy", - . "Physiq dQl","-", - . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, - . typeval, zsto,zout) -cIM 280405 BEG -c -c Champs 2D: -c -c u850, v850 -c DO k=1, nlevSTD - DO k=1, 12 -c - IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k) - IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k) -c - IF(bb2.EQ."850") THEN -c - CALL histdef(nid_bilKPins, "u"//bb2, - . "Zonal wind "//bb2//"mb","m/s", - . iim,jjphy_nb,nhori, 1,1,1, -99, 32, - . typeval, zsto,zout) -c - CALL histdef(nid_bilKPins, "v"//bb2, - . "Meridional wind "//bb2//"mb","m/s", - . iim,jjphy_nb,nhori, 1,1,1, -99, 32, - . typeval, zsto,zout) -c - ENDIF !(bb2.EQ."850") -c - ENDDO !k=1, 12 -c -cIM 280405 END -c - CALL histend(nid_bilKPins) -c - ndex2d = 0 - ndex3d = 0 -c - ENDIF ! fin de test sur ok_journe Index: LMDZ6/trunk/libf/phylmd/ini_coord_REGDYN.h =================================================================== --- LMDZ6/trunk/libf/phylmd/ini_coord_REGDYN.h (revision 5310) +++ (revision ) @@ -1,81 +1,0 @@ -c -c $Header$ -c - nsrf=3 - DO nreg=1, nbregdyn - DO i=1, klon - -c IF (debut) THEN - IF(rlon(i).LT.0.) THEN - rlonPOS(i)=rlon(i)+360. - ELSE - rlonPOS(i)=rlon(i) - ENDIF -c ENDIF - - pct_ocean(i,nreg)=0 - -c test si c'est 1 point d'ocean - IF(pctsrf(i,nsrf).EQ.1.) THEN - - IF(nreg.EQ.1) THEN - -c TROP - IF(rlat(i).GE.-30.AND.rlat(i).LE.30.) THEN - pct_ocean(i,nreg)=1 - ENDIF - -c PACIFIQUE NORD - ELSEIF(nreg.EQ.2) THEN - IF(rlat(i).GE.40.AND.rlat(i).LE.60.) THEN - IF(rlonPOS(i).GE.160..AND.rlonPOS(i).LE.235.) THEN - pct_ocean(i,nreg)=1 - ENDIF - ENDIF -c CALIFORNIE ST-CU - ELSEIF(nreg.EQ.3) THEN - IF(rlonPOS(i).GE.220..AND.rlonPOS(i).LE.250.) THEN - IF(rlat(i).GE.15.AND.rlat(i).LE.35.) THEN - pct_ocean(i,nreg)=1 - ENDIF - ENDIF -c HAWAI - ELSEIF(nreg.EQ.4) THEN - IF(rlonPOS(i).GE.180..AND.rlonPOS(i).LE.220.) THEN - IF(rlat(i).GE.15.AND.rlat(i).LE.35.) THEN - pct_ocean(i,nreg)=1 - ENDIF - ENDIF -c WARM POOL - ELSEIF(nreg.EQ.5) THEN - IF(rlonPOS(i).GE.70..AND.rlonPOS(i).LE.150.) THEN - IF(rlat(i).GE.-5.AND.rlat(i).LE.20.) THEN - pct_ocean(i,nreg)=1 - ENDIF - ENDIF - ENDIF !nbregdyn -c TROP -c IF(rlat(i).GE.-30.AND.rlat(i).LE.30.) THEN -c pct_ocean(i)=.TRUE. -c WRITE(*,*) 'pct_ocean =',i, rlon(i), rlat(i) -c ENDIF !lon -c ENDIF !lat - - ENDIF !pctsrf - ENDDO !klon - ENDDO !nbregdyn -cIM 190504 ENDIF !ok_regdyn - -cIM somme de toutes les nhistoW BEG - IF (debut) THEN - DO nreg = 1, nbregdyn - DO k = 1, kmaxm1 - DO l = 1, lmaxm1 - DO iw = 1, iwmax - nhistoWt(k,l,iw,nreg)=0. - ENDDO !iw - ENDDO !l - ENDDO !k - ENDDO !nreg - ENDIF !(debut) THEN -cIM 190504 BEG Index: LMDZ6/trunk/libf/phylmd/ini_histREGDYN.h =================================================================== --- LMDZ6/trunk/libf/phylmd/ini_histREGDYN.h (revision 5310) +++ (revision ) @@ -1,127 +1,0 @@ -! -! $Header$ -! - - IF (ok_regdyn) THEN - - if (is_sequential) then -c -cIM PRINT*, 'La frequence de sortie REGDYN est de ', ecrit_mth -c - idayref = day_ref - CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) -c -c axe vertical pour les differents niveaux des histogrammes - DO iw=1, iwmax - zx_o500(iw)=wmin+(iw-1./2.)*pas_w - ENDDO - - CALL histbeg("histREGDYN", kmaxm1,zx_tau, lmaxm1,zx_pc, - . 1,kmaxm1,1,lmaxm1, itau_phy, zjulian, dtime, - . nhoriRD, nid_regdyn) - - CALL histvert(nid_regdyn, "omeganivs", "Omega levels", - . "mb/day", - . iwmax, zx_o500, komega) -c -c pour les champs instantannes, il faut mettre la meme valeur pour -c zout et zsto. -c dtime est passe par ailleurs a histbeg -c -c zout = dtime * REAL(NINT(86400./dtime*ecrit_regdyn)) -c zsto = zout -c print*,'zout,zsto=',zout,zsto -c -c stockage a chaque pas de temps de la physique -c - zstophy = dtime -cIM 020904 zstophy = dtime * nbapp_isccp - -c ecriture mensuelle -c - zout = dtime * ecrit_mth -cIM 020904 -c zout = dtime * ecrit_day -c zout = dtime * REAL(NINT(86400./dtime*ecrit_regdyn)) - -c -c Champs 3D: -c -c TROP - CALL histdef(nid_regdyn, "hw1", "Tropics Histogram ", "%", - & kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, - & "ave(X)", zstophy,zout) - - CALL histdef(nid_regdyn, "nh1", "Nb of pixels Tropics Histo", - & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, - & 32,"ave(X)", zstophy,zout) -c - - CALL histdef(nid_regdyn, "nht1", - & "Total Nb pixels Tropics Histo", - & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, - & 32,"ave(X)", zstophy,zout) -c -c PAN - CALL histdef(nid_regdyn, "hw2", "North Pacific Histogram", "%", - & kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, - & "ave(X)", zstophy,zout) - - CALL histdef(nid_regdyn, "nh2", "Nb of pixels North Pacific", - & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, - & 32,"ave(X)", zstophy,zout) -c - - CALL histdef(nid_regdyn, "nht2", - & "Total Nb pixels North Pacific Histo", - & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, - & 32,"ave(X)", zstophy,zout) -c CAL - CALL histdef(nid_regdyn, "hw3", "California Histogram", "%", - & kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, - & "ave(X)", zstophy,zout) - - CALL histdef(nid_regdyn, "nh3", - & "Nb of pixels California Histo", - & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, - & 32,"ave(X)", zstophy,zout) -c - - CALL histdef(nid_regdyn, "nht3", - & "Total Nb pixels California Histo", - & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, - & 32,"ave(X)", zstophy,zout) -c HAW - CALL histdef(nid_regdyn, "hw4", "Hawai Histogram", "%", - & kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, - & "ave(X)", zstophy,zout) - - CALL histdef(nid_regdyn, "nh4", "Nb of pixels Hawai Histo", - & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, - & 32,"ave(X)", zstophy,zout) -c - - CALL histdef(nid_regdyn, "nht4", - & "Total Nb pixels Hawai Histo", - & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, - & 32,"ave(X)", zstophy,zout) -c WAP - CALL histdef(nid_regdyn, "hw5", "Warm Pool Histogram", "%", - & kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, 32, - & "ave(X)", zstophy,zout) - - CALL histdef(nid_regdyn, "nh5", "Nb of pixels Warm Pool Histo", - & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, - & 32,"ave(X)", zstophy,zout) -c - - CALL histdef(nid_regdyn, "nht5", - & "Total Nb pixels Warm Pool Histo", - & "%",kmaxm1,lmaxm1,nhoriRD, iwmax,1,iwmax, komega, - & 32,"ave(X)", zstophy,zout) -c - CALL histend(nid_regdyn) - - endif ! is_sequential - - endif ! ok_regdyn Index: LMDZ6/trunk/libf/phylmd/ini_histday_seri.h =================================================================== --- LMDZ6/trunk/libf/phylmd/ini_histday_seri.h (revision 5310) +++ (revision ) @@ -1,131 +1,0 @@ -! -! $Id$ -! -!ym Ne fonctionnera pas en mode parallele - IF (is_sequential) THEN - - IF (type_run.EQ."AMIP") THEN -! - zstophy = pdtphys - zout = ecrit_day -! - idayref = day_ref - CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) -! - CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,longitude_deg,zx_lon) - DO i = 1, nbp_lon - zx_lon(i,1) = longitude_deg(i+1) - zx_lon(i,nbp_lat) = longitude_deg(i+1) - ENDDO - DO ll=1,klev - znivsig(ll)=REAL(ll) - ENDDO - CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,latitude_deg,zx_lat) -! - imin_debut=1 - nbpti=1 - jmin_debut=1 - nbptj=1 -! - CALL histbeg("histday_seri.nc", & - nbp_lon,zx_lon(:,1),nbp_lat,zx_lat(1,:), & - imin_debut,nbpti,jmin_debut,nbptj, & - itau_phy, zjulian, dtime, & - nhori, nid_day_seri) -! - CALL histvert(nid_day_seri, "presnivs", & - "Vertical levels","mb", & - klev, presnivs/100., nvert) -! - CALL histdef(nid_day_seri, "bilTOA", & - "Net radiation at model top", "W/m2", & - nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & - "ave(X)", zstophy,zout) -! - CALL histdef(nid_day_seri, "bils", & - "Net downward energy flux at surface","W/m2", & - nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & - "ave(X)", zstophy,zout) -! - CALL histdef(nid_day_seri, "ecin", & - "Total kinetic energy (per unit area)","J/m2", & - nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & - "ave(X)", zstophy,zout) -! -!IM 151004 BEG - IF(1.EQ.0) THEN -! - CALL histdef(nid_day_seri, "momang", & - "Total relative angular momentum (per unit area)", & - "kg/s", & - nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & - "ave(X)", zstophy,zout) -! - CALL histdef(nid_day_seri, "frictor", & - "Friction torque (per unit area)", "N/m", & - nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & - "ave(X)", zstophy,zout) -! - CALL histdef(nid_day_seri, "mountor", & - "Mountain torque (per unit area)", "N/m", & - nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & - "ave(X)", zstophy,zout) -! - ENDIF !(1.EQ.0) THEN -! - CALL histdef(nid_day_seri, "momang", & - "Axial angular momentum (per unit area)", & - "kg/s", & - nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & - "ave(X)", zstophy,zout) -! - CALL histdef(nid_day_seri, "torsfc", & - "Total surface torque (including mountain torque)", "N/m", & - nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & - "ave(X)", zstophy,zout) -! -!IM 151004 END -! - CALL histdef(nid_day_seri, "tamv", & - "Temperature (mass-weighted vert. ave)", "K", & - nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & - "ave(X)", zstophy,zout) -! - CALL histdef(nid_day_seri, "psol", & - "Surface pressure", "Pa", & - nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & - "ave(X)", zstophy,zout) -! - CALL histdef(nid_day_seri, "evap", & - "Evaporation and sublimation (per unit area)", & - "kg/(m2*s)", & - nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & - "ave(X)", zstophy,zout) -! -! call histdef(nid_day_seri, -! . "SnowFrac", -! . "Snow-covered area ", "%", -! . nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, -! . "ave(X)", zstophy,zout) -! -! CALL histdef(nid_day_seri, "snow_depth", -!IM 080904 . "Snow Depth (water equivalent)", "m", -!IM 191104 . "Snow Depth (water equivalent)", "kg/m2", -! . "Snow Mass", "kg/m2", -! . nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, -! . "ave(X)", zstophy,zout) -! - call histdef(nid_day_seri, & - "tsol_"//clnsurf(is_oce), & - "SST over open (ice-free) ocean ", "K", & - nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, & - "ave(X)", zstophy,zout) -! -!================================================================= -! - CALL histend(nid_day_seri) -! -!================================================================= - ENDIF ! fin de test sur type_run.EQ.AMIP - - ENDIF ! is_sequential Index: LMDZ6/trunk/libf/phylmdiso/ini_bilKP_ave.h =================================================================== --- LMDZ6/trunk/libf/phylmdiso/ini_bilKP_ave.h (revision 5310) +++ (revision ) @@ -1,1 +1,0 @@ -link ../phylmd/ini_bilKP_ave.h Index: LMDZ6/trunk/libf/phylmdiso/ini_bilKP_ins.h =================================================================== --- LMDZ6/trunk/libf/phylmdiso/ini_bilKP_ins.h (revision 5310) +++ (revision ) @@ -1,1 +1,0 @@ -link ../phylmd/ini_bilKP_ins.h Index: LMDZ6/trunk/libf/phylmdiso/ini_coord_REGDYN.h =================================================================== --- LMDZ6/trunk/libf/phylmdiso/ini_coord_REGDYN.h (revision 5310) +++ (revision ) @@ -1,1 +1,0 @@ -link ../phylmd/ini_coord_REGDYN.h Index: LMDZ6/trunk/libf/phylmdiso/ini_histREGDYN.h =================================================================== --- LMDZ6/trunk/libf/phylmdiso/ini_histREGDYN.h (revision 5310) +++ (revision ) @@ -1,1 +1,0 @@ -link ../phylmd/ini_histREGDYN.h Index: LMDZ6/trunk/libf/phylmdiso/ini_histday_seri.h =================================================================== --- LMDZ6/trunk/libf/phylmdiso/ini_histday_seri.h (revision 5310) +++ (revision ) @@ -1,1 +1,0 @@ -link ../phylmd/ini_histday_seri.h