! ! $Id: initfluxsto.f90 5285 2024-10-28 13:33:29Z abarral $ ! subroutine initfluxsto & (infile,tstep,t_ops,t_wrt, & fileid,filevid,filedid) USE iniprint_mod_h USE comgeom_mod_h USE IOIPSL USE comconst_mod, ONLY: pi USE comvert_mod, ONLY: nivsigs USE temps_mod, ONLY: annee_ref, day_ref, itau_dyn USE dimensions_mod, ONLY: iim, jjm, llm, ndm USE paramet_mod_h implicit none ! ! Routine d'initialisation des ecritures des fichiers histoires LMDZ ! au format IOIPSL ! ! Appels succesifs des routines: histbeg ! histhori ! histver ! histdef ! histend ! ! Entree: ! ! infile: nom du fichier histoire a creer ! day0,anne0: date de reference ! tstep: duree du pas de temps en seconde ! t_ops: frequence de l'operation pour IOIPSL ! t_wrt: frequence d'ecriture sur le fichier ! ! Sortie: ! fileid: ID du fichier netcdf cree ! filevid:ID du fichier netcdf pour la grille v ! ! L. Fairhead, LMD, 03/99 ! ! ===================================================================== ! ! Declarations ! Arguments ! character(len=*) :: infile real :: tstep, t_ops, t_wrt integer :: fileid, filevid,filedid ! This routine needs IOIPSL to work ! Variables locales ! real :: nivd(1) integer :: tau0 real :: zjulian character(len=3) :: str character(len=10) :: ctrac integer :: iq real :: rlong(iip1,jjp1), rlat(iip1,jjp1),rl(1,1) integer :: uhoriid, vhoriid, thoriid, zvertiid,dhoriid,dvertiid integer :: ii,jj integer :: zan, idayref logical :: ok_sync ! ! Initialisations ! pi = 4. * atan (1.) str='q ' ctrac = 'traceur ' ok_sync = .true. ! ! Appel a histbeg: creation du fichier netcdf et initialisations diverses ! zan = annee_ref idayref = day_ref CALL ymds2ju(zan, 1, idayref, 0.0, zjulian) tau0 = itau_dyn do jj = 1, jjp1 do ii = 1, iip1 rlong(ii,jj) = rlonu(ii) * 180. / pi rlat(ii,jj) = rlatu(jj) * 180. / pi enddo enddo call histbeg(infile, iip1, rlong(:,1), jjp1, rlat(1,:), & 1, iip1, 1, jjp1, & tau0, zjulian, tstep, uhoriid, fileid) ! ! Creation du fichier histoire pour la grille en V (oblige pour l'instant, ! IOIPSL ne permet pas de grilles avec des nombres de point differents dans ! un meme fichier) do jj = 1, jjm do ii = 1, iip1 rlong(ii,jj) = rlonv(ii) * 180. / pi rlat(ii,jj) = rlatv(jj) * 180. / pi enddo enddo call histbeg('fluxstokev.nc', iip1, rlong(:,1), jjm, rlat(1,:), & 1, iip1, 1, jjm, & tau0, zjulian, tstep, vhoriid, filevid) rl(1,1) = 1. call histbeg('defstoke.nc', 1, rl, 1, rl, & 1, 1, 1, 1, & tau0, zjulian, tstep, dhoriid, filedid) ! ! Appel a histhori pour rajouter les autres grilles horizontales ! do jj = 1, jjp1 do ii = 1, iip1 rlong(ii,jj) = rlonv(ii) * 180. / pi rlat(ii,jj) = rlatu(jj) * 180. / pi enddo enddo call histhori(fileid, iip1, rlong, jjp1, rlat, 'scalar', & 'Grille points scalaires', thoriid) ! ! Appel a histvert pour la grille verticale ! call histvert(fileid, 'sig_s', 'Niveaux sigma', & 'sigma_level', & llm, nivsigs, zvertiid) ! Pour le fichier V call histvert(filevid, 'sig_s', 'Niveaux sigma', & 'sigma_level', & llm, nivsigs, zvertiid) ! pour le fichier def nivd(1) = 1 call histvert(filedid, 'sig_s', 'Niveaux sigma', & 'sigma_level', & 1, nivd, dvertiid) ! ! Appels a histdef pour la definition des variables a sauvegarder CALL histdef(fileid, "phis", "Surface geop. height", "-", & iip1,jjp1,thoriid, 1,1,1, -99, 32, & "once", t_ops, t_wrt) CALL histdef(fileid, "aire", "Grid area", "-", & iip1,jjp1,thoriid, 1,1,1, -99, 32, & "once", t_ops, t_wrt) CALL histdef(filedid, "dtvr", "tps dyn", "s", & 1,1,dhoriid, 1,1,1, -99, 32, & "once", t_ops, t_wrt) CALL histdef(filedid, "istdyn", "tps stock", "s", & 1,1,dhoriid, 1,1,1, -99, 32, & "once", t_ops, t_wrt) CALL histdef(filedid, "istphy", "tps stock phy", "s", & 1,1,dhoriid, 1,1,1, -99, 32, & "once", t_ops, t_wrt) ! ! Masse ! call histdef(fileid, 'masse', 'Masse', 'kg', & iip1, jjp1, thoriid, llm, 1, llm, zvertiid, & 32, 'inst(X)', t_ops, t_wrt) ! ! Pbaru ! call histdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s', & iip1, jjp1, uhoriid, llm, 1, llm, zvertiid, & 32, 'inst(X)', t_ops, t_wrt) ! ! Pbarv ! call histdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s', & iip1, jjm, vhoriid, llm, 1, llm, zvertiid, & 32, 'inst(X)', t_ops, t_wrt) ! ! w ! call histdef(fileid, 'w', 'flx de masse vert', 'kg m/s', & iip1, jjp1, thoriid, llm, 1, llm, zvertiid, & 32, 'inst(X)', t_ops, t_wrt) ! ! Temperature potentielle ! call histdef(fileid, 'teta', 'temperature potentielle', '-', & iip1, jjp1, thoriid, llm, 1, llm, zvertiid, & 32, 'inst(X)', t_ops, t_wrt) ! ! ! Geopotentiel ! call histdef(fileid, 'phi', 'geopotentiel instantane', '-', & iip1, jjp1, thoriid, llm, 1, llm, zvertiid, & 32, 'inst(X)', t_ops, t_wrt) ! ! Fin ! call histend(fileid) call histend(filevid) call histend(filedid) if (ok_sync) then call histsync(fileid) call histsync(filevid) call histsync(filedid) endif return end subroutine initfluxsto