| 1 | !$Id $ |
|---|
| 2 | !*************************************** |
|---|
| 3 | ! ECRITURE DU FICHIER : histrac.nc |
|---|
| 4 | !*************************************** |
|---|
| 5 | IF (ecrit_tra > 0. .AND. config_inca == 'none') THEN |
|---|
| 6 | |
|---|
| 7 | itau_w = itau_phy + nstep |
|---|
| 8 | |
|---|
| 9 | CALL histwrite_phy(nid_tra,"phis",itau_w,pphis) |
|---|
| 10 | CALL histwrite_phy(nid_tra,"aire",itau_w,airephy) |
|---|
| 11 | |
|---|
| 12 | !TRACEURS |
|---|
| 13 | !---------------- |
|---|
| 14 | DO it=1,nbtr |
|---|
| 15 | iiq=niadv(it+2) |
|---|
| 16 | |
|---|
| 17 | ! CONCENTRATIONS |
|---|
| 18 | CALL histwrite_phy(nid_tra,tname(iiq),itau_w,tr_seri(:,:,it)) |
|---|
| 19 | |
|---|
| 20 | ! TD LESSIVAGE |
|---|
| 21 | IF (lessivage .AND. aerosol(it)) THEN |
|---|
| 22 | CALL histwrite_phy(nid_tra,"fl"//tname(iiq),itau_w,flestottr(:,:,it)) |
|---|
| 23 | ENDIF |
|---|
| 24 | |
|---|
| 25 | ! TD THERMIQUES |
|---|
| 26 | IF (iflag_thermals.gt.0) THEN |
|---|
| 27 | CALL histwrite_phy(nid_tra,"d_tr_th_"//tname(iiq),itau_w,d_tr_th(:,:,it)) |
|---|
| 28 | ENDIF |
|---|
| 29 | |
|---|
| 30 | ! TD CONVECTION |
|---|
| 31 | IF (iflag_con.GE.2) THEN |
|---|
| 32 | CALL histwrite_phy(nid_tra,"d_tr_cv_"//tname(iiq),itau_w,d_tr_cv(:,:,it)) |
|---|
| 33 | ENDIF |
|---|
| 34 | |
|---|
| 35 | ! TD COUCHE-LIMITE |
|---|
| 36 | CALL histwrite_phy(nid_tra,"d_tr_cl_"//tname(iiq),itau_w,d_tr_cl(:,:,it)) |
|---|
| 37 | ENDDO |
|---|
| 38 | !--------------- |
|---|
| 39 | ! |
|---|
| 40 | ! |
|---|
| 41 | ! VENT (niveau 1) |
|---|
| 42 | CALL histwrite_phy(nid_tra,"pyu1",itau_w,yu1) |
|---|
| 43 | CALL histwrite_phy(nid_tra,"pyv1",itau_w,yv1) |
|---|
| 44 | ! |
|---|
| 45 | ! TEMPERATURE DU SOL |
|---|
| 46 | zx_tmp_fi2d(:)=ftsol(:,1) |
|---|
| 47 | CALL histwrite_phy(nid_tra,"ftsol1",itau_w,zx_tmp_fi2d) |
|---|
| 48 | zx_tmp_fi2d(:)=ftsol(:,2) |
|---|
| 49 | CALL histwrite_phy(nid_tra,"ftsol2",itau_w,zx_tmp_fi2d) |
|---|
| 50 | zx_tmp_fi2d(:)=ftsol(:,3) |
|---|
| 51 | CALL histwrite_phy(nid_tra,"ftsol3",itau_w,zx_tmp_fi2d) |
|---|
| 52 | zx_tmp_fi2d(:)=ftsol(:,4) |
|---|
| 53 | CALL histwrite_phy(nid_tra,"ftsol4",itau_w,zx_tmp_fi2d) |
|---|
| 54 | ! |
|---|
| 55 | ! NATURE DU SOL |
|---|
| 56 | zx_tmp_fi2d(:)=pctsrf(:,1) |
|---|
| 57 | CALL histwrite_phy(nid_tra,"psrf1",itau_w,zx_tmp_fi2d) |
|---|
| 58 | zx_tmp_fi2d(:)=pctsrf(:,2) |
|---|
| 59 | CALL histwrite_phy(nid_tra,"psrf2",itau_w,zx_tmp_fi2d) |
|---|
| 60 | zx_tmp_fi2d(:)=pctsrf(:,3) |
|---|
| 61 | CALL histwrite_phy(nid_tra,"psrf3",itau_w,zx_tmp_fi2d) |
|---|
| 62 | zx_tmp_fi2d(:)=pctsrf(:,4) |
|---|
| 63 | CALL histwrite_phy(nid_tra,"psrf4",itau_w,zx_tmp_fi2d) |
|---|
| 64 | |
|---|
| 65 | ! DIVERS |
|---|
| 66 | CALL histwrite_phy(nid_tra,"pplay",itau_w,pplay) |
|---|
| 67 | CALL histwrite_phy(nid_tra,"t",itau_w,t_seri) |
|---|
| 68 | CALL histwrite_phy(nid_tra,"mfu",itau_w,pmfu) |
|---|
| 69 | CALL histwrite_phy(nid_tra,"mfd",itau_w,pmfd) |
|---|
| 70 | CALL histwrite_phy(nid_tra,"en_u",itau_w,pen_u) |
|---|
| 71 | CALL histwrite_phy(nid_tra,"en_d",itau_w,pen_d) |
|---|
| 72 | CALL histwrite_phy(nid_tra,"de_d",itau_w,pde_d) |
|---|
| 73 | CALL histwrite_phy(nid_tra,"de_u",itau_w,pde_u) |
|---|
| 74 | CALL histwrite_phy(nid_tra,"coefh",itau_w,coefh) |
|---|
| 75 | |
|---|
| 76 | IF (ok_sync) THEN |
|---|
| 77 | !$OMP MASTER |
|---|
| 78 | CALL histsync(nid_tra) |
|---|
| 79 | !$OMP END MASTER |
|---|
| 80 | ENDIF |
|---|
| 81 | |
|---|
| 82 | ENDIF !ecrit_tra>0. .AND. config_inca == 'none' |
|---|
| 83 | |
|---|
