subroutine invlat (flg,xlamb,ier) c* c*********************************************************************** c* c* "Copyright [c] CNES 98 - tous droits reserves" c* ********************************************** c* c*PRO MAGLIB c* c*VER 99.03.31 - V 1.0 c*VER 01.06.05 - V 2.0 c*VER 03.01.06 - V 2.1 c* c*AUT spec. CNES - JC KOSIK - janvier 1991 c*AUT port. CISI c* c*ROL Theme : Calculs de geophysique c*ROL Calcul de la valeur absolue de la latitude invariante c*ROL a partir du L de Galperin. c* c*PAR flg (I) : parametre L de Galperin c* c*PAR xlamb (O) : latitude invariante (radians) c* c*PAR ier (O) : code de retour c* c*NOT ier : 1 = parametre flg = 0.0d0 c* c*INF utilise : sans objet c* c*HST version 1.0 - 99.03.31 - creation de la maglib au CDPP c*HST version 2.0 - 01.06.05 - correction de commentaires de code c*HST version 2.1 - 03.01.06 - corrections en compilation avec g77 c* c*********************************************************************** c* implicit none c c --------------------------------- c*FON Declaration identificateur rcs_id c --------------------------------- c character rcs_id*100 c c -------------------------- c*FON Declaration des parametres c -------------------------- c double precision flg double precision xlamb integer ier c c --------------------------------- c*FON Declaration des variables locales c --------------------------------- c double precision cxl c*LOC Variables de travail intermediaires c SAVE c c --------------------------------- c*FON Affectation identificateur rcs_id c --------------------------------- c data rcs_id /" >$Id$"/ c c ****************** c Debut de programme c ****************** c ier = 0 c if (flg .ne. 0.0d0) then cxl = sqrt(1.d0 / flg) xlamb = acos(cxl) else ier = 1 endif c c **************** c Fin de programme c **************** c return end