source: LMDZ4/branches/LMDZ4-dev/libf/phylmd/write_histrac.h @ 1191

Last change on this file since 1191 was 1191, checked in by jghattas, 15 years ago

Reecriture de phytrac et les routines concernes (Anthony Jamelot)

  • les suffix change de F -> F90 (nflxtr.F90,cltracrn.F90,initrrnpb.F90,cvltr.F90,minmaxqfi.F90,cltrac.F90,phytrac.F90)

Traitement d'un nouveau traceur berelium (optionel, toujours pour des
tests)(Anthony Jamelot)

  • radiornpb.F change du nom pour radio_decay.F90 car il traite maintenant tout les traceurs radioactives
  • ajoute init_be.F90

Nouveau interface dans phytrac pour serparer les calculs et appels
specifique a INCA avec les traitements des traceurs specifiques au LMDZ
(JG)

  • ajoute tracinca_mod.F90 pour les appeles a INCA
  • ajoute traclmdz_mod.F90 pour les calculs des traceurs specifiques a LMDZ
  • enleve fichier restartrac et ajoute la variable trs dans restartphy.nc

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